Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,017,429 (C)6→7 coding (876/1374 nt) fliI → flagellum‑specific ATP synthase FliI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,4231.C96.9% 221.1 / ‑0.4 64coding (870/1374 nt)fliIflagellum‑specific ATP synthase FliI
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (34/28);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CTGCGAATGCAGGGTGCCGCCTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCG  >  NC_000913/2017286‑2017550
                                                                                                                                           |                                                                                                                               
cTGCGAATGCAGGGTGCCGCCTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAAcc                                                                                                                               <  1:170458/140‑1 (MQ=255)
  gCGAATGCAGGGTGCCGCCTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAAcccc                                                                                                                             >  2:2807/1‑140 (MQ=255)
      aTGCAGGGTGCCACCTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAAcccccccg                                                                                                                         <  2:283694/140‑3 (MQ=255)
                ccgccTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAgg                                                                                                                >  1:1038401/1‑140 (MQ=255)
                ccgccTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAgg                                                                                                                <  2:1038401/140‑1 (MQ=255)
                   ccTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGTCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTa                                                                                                             >  1:378070/1‑140 (MQ=255)
                   ccTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTa                                                                                                             >  1:443456/1‑140 (MQ=255)
                      aTGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATcc                                                                                                          >  1:715086/1‑140 (MQ=255)
                       tGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATc                                                                                                           >  2:16370/1‑138 (MQ=255)
                       tGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATc                                                                                                           <  1:16370/138‑1 (MQ=255)
                             gcgcATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCaa                                                                                                                   >  1:341774/1‑124 (MQ=255)
                             gcgcATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCaa                                                                                                                   <  2:341774/124‑1 (MQ=255)
                             gcgcATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGt                                                                                                     <  2:491177/138‑1 (MQ=255)
                             gcgcATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGt                                                                                                     >  1:491177/1‑138 (MQ=255)
                              cgcATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCgg                                                                                                  >  1:609954/1‑140 (MQ=255)
                                        aaGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGTCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCaa                                                                                        <  2:378070/140‑1 (MQ=255)
                                        aaGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCaa                                                                                        <  2:609954/140‑1 (MQ=255)
                                        aaGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCaa                                                                                        <  2:443456/140‑1 (MQ=255)
                                        aaGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCaa                                                                                        <  2:790786/140‑1 (MQ=255)
                                        aaGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCaa                                                                                        <  2:715086/140‑1 (MQ=255)
                                        aaGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCaa                                                                                        >  2:605202/1‑140 (MQ=255)
                                                cgcgATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCgg                                                                                <  2:1047558/140‑1 (MQ=255)
                                                  cgATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCa                                                                              <  1:437641/140‑1 (MQ=255)
                                                      cGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTgg                                                                          >  1:144134/1‑140 (MQ=255)
                                                               cATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTg                                                                 >  1:125304/1‑140 (MQ=255)
                                                               cATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTg                                                                 >  1:374918/1‑140 (MQ=255)
                                                                aTGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGc                                                                <  2:917088/140‑1 (MQ=255)
                                                                 tgtgTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGc                                                                <  2:707979/139‑1 (MQ=255)
                                                                 tgtgTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGc                                                                >  1:707979/1‑139 (MQ=255)
                                                                  gtgtTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCg                                                              >  1:732473/1‑140 (MQ=255)
                                                                    gtTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGa                                                            >  1:720231/1‑140 (MQ=255)
                                                                    gtTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGa                                                            >  1:898569/1‑140 (MQ=255)
                                                                     ttGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGaa                                                           >  2:945801/1‑140 (MQ=255)
                                                                       gCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAAt                                                         >  2:245165/1‑140 (MQ=255)
                                                                       gCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAAt                                                         >  2:208916/1‑140 (MQ=255)
                                                                         tGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATgg                                                       >  2:611759/1‑140 (MQ=255)
                                                                           aTTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCa                                                     >  2:1014103/1‑140 (MQ=255)
                                                                                  aCTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAgcgg                                              >  1:53461/1‑140 (MQ=255)
                                                                                      ccTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGaa                                                           >  2:598245/1‑123 (MQ=255)
                                                                                      ccTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGaa                                                           <  1:598245/123‑1 (MQ=255)
                                                                                      ccTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAtcggcggc                                          >  1:860609/1‑140 (MQ=255)
                                                                                      ccTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAgcggcggc                                          >  2:169674/1‑140 (MQ=255)
                                                                                      ccTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAgcggcggc                                          >  2:604894/1‑140 (MQ=255)
                                                                                       cTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAgcggcggcg                                         >  1:533121/1‑140 (MQ=255)
                                                                                        tCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGcggcggcgg                                        <  1:939152/140‑1 (MQ=255)
                                                                                         cACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCggcggcggc                                       >  2:442810/1‑140 (MQ=255)
                                                                                                tACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTa                                >  1:903901/1‑140 (MQ=255)
                                                                                                  cgcgATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              <  1:945801/140‑1 (MQ=255)
                                                                                                  cgcgATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              <  2:53461/140‑1 (MQ=255)
                                                                                                  cgcgATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              <  2:533121/140‑1 (MQ=255)
                                                                                                  cgcgATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              <  1:442810/140‑1 (MQ=255)
                                                                                                  cgcgATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              <  2:420423/140‑1 (MQ=255)
                                                                                                  cgcgATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              <  2:374918/140‑1 (MQ=255)
                                                                                                      aTGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              <  2:54679/135‑1 (MQ=255)
                                                                                                      aTGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTAcc                              >  1:54679/1‑135 (MQ=255)
                                                                                                             aGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTtata                    <  1:644924/139‑1 (MQ=255)
                                                                                                             aGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTtata                    >  2:644924/1‑139 (MQ=255)
                                                                                                              gCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATAcc                  <  2:898569/140‑1 (MQ=255)
                                                                                                                 tGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTg               <  1:611759/140‑1 (MQ=255)
                                                                                                                 tGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTg               <  2:125304/140‑1 (MQ=255)
                                                                                                                 tGAAATTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTg               <  1:245165/140‑1 (MQ=255)
                                                                                                                 ggAAATGGCGCTGGCGATTGGCGAGCCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTg               <  1:47332/139‑1 (MQ=255)
                                                                                                                   aaaTTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCt             >  2:733832/1‑140 (MQ=255)
                                                                                                                   aaaTTGCGCTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGGCGGTGTTTGCCAAATTACCGTCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGAGGCTCGATTACCCCGGTTTATACCGtttt             >  2:47332/1‑137 (MQ=255)
                                                                                                                         cgcTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTg               <  1:404947/132‑1 (MQ=255)
                                                                                                                         cgcTGGCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTg               >  2:404947/1‑132 (MQ=255)
                                                                                                                              gCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCg  <  1:2807/140‑1 (MQ=255)
                                                                                                                              gCGATTGGCGAACCCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTGTATACCGTGCTCACTGAAGGCg  >  1:741556/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                               
CTGCGAATGCAGGGTGCCGCCTATGCCACGCGCATTGCCGAAGATTTTCGCGATCGTGGTCAGCATGTGTTGCTGATTATGGACTCCCTCACCCGCTACGCGATGGCCCAGCGTGAAATTGCGCTGGCGATTGGCGAACCCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCG  >  NC_000913/2017286‑2017550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: