Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,532,298 A→T T15T (ACA→ACT murN → FemAB family peptidoglycan biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,532,2980AT94.3% 209.0 / 3.5 70T15T (ACA→ACTmurNFemAB family peptidoglycan biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base T (25/41);  total (25/45)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

ATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532166‑2532433
                                                                                                                                    |                                                                                                                                       
aTTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCa                                                                                                                                 <  2:671003/141‑1 (MQ=255)
   tGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAttt                                                                                                                              <  2:600404/141‑1 (MQ=255)
   tGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAttt                                                                                                                              >  1:818914/1‑141 (MQ=255)
     tAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGTTGGAATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAttttt                                                                                                                            >  2:282540/1‑141 (MQ=255)
      aaTCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTc                                                                                                                           <  1:460623/141‑1 (MQ=255)
       aTCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTACAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTca                                                                                                                          <  1:916750/141‑1 (MQ=255)
        tCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTcac                                                                                                                         >  1:903973/1‑141 (MQ=255)
          aGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACAt                                                                                                                       <  1:60268/141‑1 (MQ=255)
           gTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCGCAtt                                                                                                                      <  1:286049/141‑1 (MQ=255)
             ttttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCa                                                                                                                                 <  2:547160/128‑1 (MQ=255)
             ttttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCa                                                                                                                                 >  1:547160/1‑128 (MQ=255)
             ttttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTac                                                                                                                    >  1:633947/1‑141 (MQ=255)
             ttttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTac                                                                                                                    >  2:504053/1‑141 (MQ=255)
              tttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTaca                                                                                                                   <  2:675026/141‑1 (MQ=255)
               ttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacac                                                                                                                  <  2:85840/141‑1 (MQ=255)
                tGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGAGAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacaca                                                                                                                 <  1:788198/141‑1 (MQ=255)
                 gATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacacac                                                                                                                >  2:858094/1‑141 (MQ=255)
                  ataaaataaaAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTAcacaca                                                                                                               <  1:43264/141‑1 (MQ=255)
                     aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATc                                                                                                            >  1:860854/1‑141 (MQ=255)
                     aaataaaaAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATc                                                                                                            <  2:731134/141‑1 (MQ=255)
                           aaaaaTAGCATTATTAATACATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAAt                                                                                                             >  1:997807/1‑134 (MQ=255)
                           aaaaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAAt                                                                                                             <  2:997807/134‑1 (MQ=255)
                             aaaTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTc                                                                                                    >  1:142038/1‑141 (MQ=255)
                               aTAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAt                                                                                                  >  1:685588/1‑141 (MQ=255)
                                tAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAtt                                                                                                 <  2:270897/141‑1 (MQ=255)
                                tAGCATTATTAATAAATTCAAAAATTACAGTTCAATTCAGAGAGGGATGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTa                                                                                                <  2:912152/141‑1 (MQ=255)
                                 aGCATTATTAATAAATTCAAAAATTACAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTa                                                                                                <  2:912275/141‑1 (MQ=255)
                                  gCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                               <  2:822141/141‑1 (MQ=255)
                                  gCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                               <  1:443066/141‑1 (MQ=255)
                                  gCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                               <  2:882371/141‑1 (MQ=255)
                                  gCATTACTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGCGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                               <  2:969347/141‑1 (MQ=255)
                                        ttaATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGACGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCGATTCATTACAATCAt                                                                                         <  1:749578/141‑1 (MQ=255)
                                          aataaATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATAAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGCTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATag                                                                                       >  2:437171/1‑141 (MQ=255)
                                            taaATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATagag                                                                                     <  2:595124/141‑1 (MQ=255)
                                             aaaTTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGt                                                                                    >  2:288981/1‑141 (MQ=255)
                                                    aaaaTTATAGTTCAATTCAGAGATGGATTGGATAAAATGAACTTTGTAACGGTGACTTCGGATGAGTTCAATGCGTTTACTACAAATCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTa                                                                             <  2:12278/141‑1 (MQ=255)
                                                    aaaaTTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTa                                                                             <  2:513764/141‑1 (MQ=255)
                                                        ttATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAg                                                                         >  2:678761/1‑141 (MQ=255)
                                                          ataGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGc                                                                       >  1:142748/1‑141 (MQ=255)
                                                                cAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtgtg                                                                 <  2:342234/141‑1 (MQ=255)
                                                                cAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAtgtg                                                                 <  2:727759/141‑1 (MQ=255)
                                                                   ttCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCAt                                                              <  2:581131/141‑1 (MQ=255)
                                                                   ttCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCAt                                                              <  1:837268/141‑1 (MQ=255)
                                                                   ttCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCAt                                                              >  2:37193/1‑141 (MQ=255)
                                                                   ttCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATGACAATCATAGAGTTGATTTAAAAGGCGATGTGCAt                                                              <  2:187116/141‑1 (MQ=255)
                                                                                gggatAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTa                                                 <  1:976750/141‑1 (MQ=255)
                                                                                   atAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAg                                              <  2:531329/141‑1 (MQ=255)
                                                                                    tAAACTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGa                                             <  1:817076/141‑1 (MQ=255)
                                                                                        aTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                         <  1:601252/141‑1 (MQ=255)
                                                                                        aTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                         <  2:32290/141‑1 (MQ=255)
                                                                                        aTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGGGCATCTTGTAGGGGTTAAAGATGAc                                         <  2:513009/141‑1 (MQ=255)
                                                                                                tGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCtt                                                           >  1:117146/1‑115 (MQ=255)
                                                                                                tGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCtt                                                           <  2:117146/115‑1 (MQ=255)
                                                                                                   aaCGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                         >  2:438633/1‑130 (MQ=255)
                                                                                                   aaCGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                         <  1:438633/130‑1 (MQ=255)
                                                                                                    aCGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTg                             <  1:843139/141‑1 (MQ=255)
                                                                                                    aCGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTg                             <  1:441382/141‑1 (MQ=255)
                                                                                                           cTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATg                                           <  2:795274/120‑1 (MQ=255)
                                                                                                           cTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATg                                           >  1:795274/1‑120 (MQ=255)
                                                                                                              cGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGAt                   <  1:523738/141‑1 (MQ=255)
                                                                                                               ggATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATg                  <  2:726081/141‑1 (MQ=255)
                                                                                                               ggATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATg                  <  2:58130/141‑1 (MQ=255)
                                                                                                                     gTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTAtt            >  2:693285/1‑141 (MQ=255)
                                                                                                                      ttCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTg           <  2:818914/141‑1 (MQ=255)
                                                                                                                        cAATGCGTTTACTACAAAGCATTTTTCTTTTTTCTCTCAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCGTGTAGGGGTTCAAGATTACAATCGTCACGTGATTGCaccatcctcattcac         >  2:482538/1‑127 (MQ=255)
                                                                                                                         aaTGCGTTTACTACAAAGCATTTTTTTTTTTTTTCACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACa        >  2:222381/1‑141 (MQ=255)
                                                                                                                          aTGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCttctgtcaat       <  2:216744/141‑9 (MQ=255)
                                                                                                                              gTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTAttatctttttc   >  1:874452/1‑132 (MQ=255)
                                                                                                                              gTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgc   >  2:262018/1‑141 (MQ=255)
                                                                                                                               tttACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  >  1:3814/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                       
ATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTATAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532166‑2532433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: