Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 934,133 G→A G563R (GGG→AGG)  ftsK → cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913934,1330GA82.7% 570.1 / 85.2 243G563R (GGG→AGG) ftsKcell division protein
Reads supporting (aligned to +/- strand):  ref base G (25/17);  new base A (112/89);  total (137/106)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 247 total aligned reads displayed.

GCACCTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTCCGCAGTTG  >  NC_000913/934043‑934225
                                                                                          |                                                                                            
 cACCTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACgaggg                                                                                           <  1:2164430/93‑1 (MQ=255)
   ccTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCt                                                                                         >  2:5295843/1‑93 (MQ=255)
   ccTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCt                                                                                         <  1:4945267/93‑1 (MQ=255)
   ccTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCt                                                                                         >  1:677284/1‑93 (MQ=255)
   ccTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCt                                                                                         >  1:1758913/1‑93 (MQ=255)
   ccTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCt                                                                                         >  1:5099874/1‑93 (MQ=255)
    cTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTg                                                                                        >  1:7440051/1‑93 (MQ=255)
     ttCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGc                                                                                       >  2:3487016/1‑93 (MQ=255)
     ttCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGc                                                                                       >  2:3398143/1‑93 (MQ=255)
     ttCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGc                                                                                       >  2:3268278/1‑93 (MQ=255)
      tCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGcc                                                                                      >  2:6815243/1‑93 (MQ=255)
         gTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCa                                                                                   >  1:3387841/1‑93 (MQ=255)
         gTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCa                                                                                   >  1:5652838/1‑93 (MQ=255)
          ttGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCaa                                                                                  >  2:3817317/1‑93 (MQ=255)
               gcagTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGtt                                                                             >  2:7090340/1‑93 (MQ=255)
                   tACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCcg                                                                         >  2:2130131/1‑93 (MQ=255)
                    aCCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCcgc                                                                        <  2:7012795/93‑1 (MQ=255)
                      cTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCgcgc                                                                      >  1:5777419/1‑93 (MQ=255)
                       tCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGcc                                                                     >  2:1655377/1‑93 (MQ=255)
                       tCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGcc                                                                     >  2:4127126/1‑93 (MQ=255)
                        ccAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCa                                                                    <  2:3577564/93‑1 (MQ=255)
                        ccAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCa                                                                    <  2:6162482/93‑1 (MQ=255)
                          aGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGt                                                                  >  2:6574621/1‑93 (MQ=255)
                           gTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTc                                                                 <  2:3714560/93‑1 (MQ=255)
                           gTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTc                                                                 >  2:3246887/1‑93 (MQ=255)
                            taaaaGCCCCTGCCGCTGGTTTCCCGCCGGCCTCCGGCGTGAAAAAAGCGAAACTGGCGACGAAGGCTGCCGGAACCGGTGCCGCGCCCGTCt                                                                >  1:3679988/4‑93 (MQ=255)
                            tAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCt                                                                >  2:4936697/1‑93 (MQ=255)
                             aGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCtt                                                               >  1:5476028/1‑93 (MQ=255)
                              gAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTc                                                              >  2:1431422/1‑93 (MQ=255)
                              gAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTc                                                              <  1:6045162/93‑1 (MQ=255)
                               aaGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCa                                                             <  2:4251021/93‑1 (MQ=255)
                               aaGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCa                                                             <  1:3774557/93‑1 (MQ=255)
                                aGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAg                                                            <  2:6931079/93‑1 (MQ=255)
                                aGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAg                                                            >  1:986992/1‑93 (MQ=255)
                                aGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAg                                                            >  1:6174019/1‑93 (MQ=255)
                                  ccgctgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTc                                                          <  2:5608335/93‑1 (MQ=255)
                                  ccgctgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTc                                                          >  1:7341510/1‑93 (MQ=255)
                                   cgctgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGCCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCt                                                         <  1:1267564/93‑1 (MQ=255)
                                    gctgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTg                                                        <  1:2599746/93‑1 (MQ=255)
                                    gctgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTg                                                        <  1:1631778/93‑1 (MQ=255)
                                     ctgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTgg                                                       >  2:6000532/1‑93 (MQ=255)
                                     ctgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTgg                                                       <  2:2502137/93‑1 (MQ=255)
                                      tgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGc                                                      <  2:1324792/93‑1 (MQ=255)
                                      tgccgctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGc                                                      >  2:802765/1‑93 (MQ=255)
                                          gctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAAt                                                  >  2:1051518/1‑93 (MQ=255)
                                          gctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAAt                                                  >  2:2638671/1‑93 (MQ=255)
                                           ctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATa                                                 >  1:1034183/1‑93 (MQ=255)
                                           ctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATa                                                 >  2:1773799/1‑93 (MQ=255)
                                           ctgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATa                                                 >  2:510742/1‑93 (MQ=255)
                                            tgTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAg                                                <  2:5777867/93‑1 (MQ=255)
                                             gTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGc                                               >  1:1134812/1‑93 (MQ=255)
                                               ttCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCgg                                             >  2:6055322/1‑93 (MQ=255)
                                                tCCCCGTAGCCATCTGGGATGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGt                                            <  1:4768932/93‑1 (MQ=255)
                                                tCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGt                                            <  1:3570599/93‑1 (MQ=255)
                                                 ccccGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTg                                           <  2:92767/93‑1 (MQ=255)
                                                    cGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGAc                                        >  1:900388/1‑93 (MQ=255)
                                                    cGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGAc                                        >  1:6977144/1‑93 (MQ=255)
                                                     gCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGAcc                                       >  2:7198412/1‑93 (MQ=255)
                                                     gCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGAcc                                       <  1:6197996/93‑1 (MQ=255)
                                                       tGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACcgc                                     >  1:6841936/1‑93 (MQ=255)
                                                       tGGCATCTGGCGTGAAAAAAGCGACACGGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACcgc                                     <  1:4330618/93‑1 (MQ=255)
                                                        ggCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACcgcg                                    <  1:3555092/93‑1 (MQ=255)
                                                         gCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGt                                   >  2:5809204/1‑93 (MQ=255)
                                                         gCATCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGt                                   >  2:6336740/1‑93 (MQ=255)
                                                            tCTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCt                                >  1:5810433/1‑93 (MQ=255)
                                                             cTGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTc                               <  1:5692615/93‑1 (MQ=255)
                                                              tGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCa                              <  2:2690542/93‑1 (MQ=255)
                                                              tGGCGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCa                              <  2:1923330/93‑1 (MQ=255)
                                                                 cGTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGt                           <  1:6668099/93‑1 (MQ=255)
                                                                  gTGAAAAAAGCGACACTGGCGTCGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTc                          >  2:3565082/1‑93 (MQ=255)
                                                                  gTGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTc                          >  2:4451652/1‑93 (MQ=255)
                                                                   tGAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCa                         >  2:463891/1‑93 (MQ=255)
                                                                    gAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCaa                        >  2:3846663/1‑93 (MQ=255)
                                                                    gAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCaa                        >  1:1468655/1‑93 (MQ=255)
                                                                    gAAAAAAGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCaa                        >  1:7191849/1‑93 (MQ=255)
                                                                      aaaaaGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAag                      >  2:6616115/1‑93 (MQ=255)
                                                                       aaaaGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAaga                     >  2:4827142/1‑93 (MQ=255)
                                                                       aaaaGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAaga                     <  1:6015809/93‑1 (MQ=255)
                                                                        aaaGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAagag                    <  1:6300452/93‑1 (MQ=255)
                                                                         aaGCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAgg                   >  1:2740122/1‑93 (MQ=255)
                                                                          aGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAggg                  <  1:413825/93‑1 (MQ=255)
                                                                          aGCGACACGGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAggg                  <  1:1393425/93‑1 (MQ=255)
                                                                           gCGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAgggg                 <  1:7157985/93‑1 (MQ=255)
                                                                            cGACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGa                <  2:5985181/93‑1 (MQ=255)
                                                                             gACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGAt               <  1:461627/93‑1 (MQ=255)
                                                                             gACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGAt               >  2:120613/1‑93 (MQ=255)
                                                                             gACACTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGAt               >  1:2205785/1‑93 (MQ=255)
                                                                              acacTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGAtt              <  1:3323477/93‑1 (MQ=255)
                                                                                acTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTgg            >  2:489517/1‑93 (MQ=255)
                                                                                acTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTgg            >  2:3615045/1‑93 (MQ=255)
                                                                                acTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTgg            <  2:2360104/93‑1 (MQ=255)
                                                                                 cTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGt           <  2:5173056/93‑1 (MQ=255)
                                                                                 cTGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGt           >  1:7243803/1‑93 (MQ=255)
                                                                                  tGGCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTc          >  1:519335/1‑93 (MQ=255)
                                                                                  nggCGACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGNc          >  2:5514073/2‑93 (MQ=255)
                                                                                      gACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTCCGCa      >  2:2904733/1‑93 (MQ=255)
                                                                                      gACGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTCCGCa      >  2:6787235/1‑93 (MQ=255)
                                                                                       aCGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTCCGCAg     >  1:3508027/1‑93 (MQ=255)
                                                                                       aCGAGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTCCGCAg     >  2:3914340/1‑93 (MQ=255)
                                                                                         gggggCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTCCGCAGtt   >  2:3026315/1‑93 (MQ=255)
                                                                                          |                                                                                            
GCACCTTCTGTTGCAGCAGTACCTCCAGTAGAAGCCGCTGCCGCTGTTTCCCCGCTGGCATCTGGCGTGAAAAAAGCGACACTGGCGACGGGGGCTGCCGCAACCGTTGCCGCGCCAGTCTTCAGTCTGGCAAATAGCGGTGGACCGCGTCCTCAGGTCAAAGAGGGGATTGGTCCGCAGTTG  >  NC_000913/934043‑934225

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