Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,972,881 A→T 9.0% F223I (TTT→ATT)  flgL ← flagellar hook‑associated protein FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,972,8810AT 9.0% 171.2 / 5.9 64F223I (TTT→ATT) flgLflagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base T (4/2);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GCAGCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTCGC  >  NC_002947/4972734‑4973028
                                                                                                                                                   |                                                                                                                                                    
gcAGCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAAc                                                                                                                                                     >  3:251531/1‑149 (MQ=255)
 cAGCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACt                                                                                                                                                    <  6:252660/149‑1 (MQ=255)
      agaTCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTg                                                                                                                                               <  7:192713/149‑1 (MQ=255)
            tcTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATCCGGCTGACCAGCGGCAAACTGATTGCCGTAg                                                                                                                                         >  4:226815/1‑149 (MQ=255)
                gccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAATCAGTTTGCCgctg                                                                                                                                         <  5:191661/145‑4 (MQ=255)
                gccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAggtgg                                                                                                                                     >  1:390139/1‑149 (MQ=255)
                gccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAg                                                                                                                                         <  4:251531/145‑1 (MQ=255)
                gccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAg                                                                                                                                         <  3:286495/145‑1 (MQ=255)
                 ccgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAggtgg                                                                                                                                     <  2:390139/148‑1 (MQ=255)
                  cgccTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGgtggtg                                                                                                                                   >  4:188686/1‑149 (MQ=255)
                        cGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCg                                                                                                                               <  5:95413/147‑1 (MQ=255)
                           tATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGt                                                                                                                           >  3:252279/1‑148 (MQ=255)
                             ttGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCaa                                                                                                                        >  2:57900/1‑149 (MQ=255)
                                aCGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAg                                                                                                                     >  8:87706/1‑149 (MQ=255)
                                aCGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAg                                                                                                                     >  2:362892/1‑149 (MQ=255)
                                  gTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgc                                                                                                                   >  3:267892/1‑149 (MQ=255)
                                  gTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATAAGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGgc                                                                                                                   >  3:358179/1‑149 (MQ=255)
                                             ccttgccGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGccaccaacaa                                                                                                          >  8:5473/1‑143 (MQ=255)
                                                           tCATTGGTGATGTCATTGCCGGCGGCATACCTCACGCTGTACTCGGAGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGt                                                                                           >  1:64999/1‑148 (MQ=255)
                                                            cATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGAc                                                                                                           >  3:211015/1‑131 (MQ=255)
                                                            cATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGAc                                                                                                           <  4:211015/131‑1 (MQ=255)
                                                             aTTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCAt                                                                                        >  5:202614/1‑149 (MQ=255)
                                                              ttGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGAGGCAAAATGATTGCAGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATTCCTTGCCGTCATc                                                                                       >  1:117712/1‑149 (MQ=255)
                                                                 gTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCgg                                                                                    >  3:188688/1‑149 (MQ=255)
                                                                  tGATGTCATTGCCGGCGGCATCCCTCACGCTGGACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGc                                                                                   >  3:268483/1‑149 (MQ=255)
                                                                       tCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACgg                                                                              <  6:293295/149‑1 (MQ=255)
                                                                       tCATTGCCGGCGGCATCCCTCACGCTGTACTAGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACgg                                                                              >  1:36432/1‑149 (MQ=255)
                                                                         aTTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTg                                                                            >  1:157759/1‑149 (MQ=255)
                                                                         aTTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTg                                                                            >  4:123855/1‑149 (MQ=255)
                                                                               ggcggcATCCCTCACGCTGCACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCa                                                                      <  4:268483/149‑1 (MQ=255)
                                                                                gcggcATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAg                                                                     >  3:305128/1‑149 (MQ=255)
                                                                                gcggcATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAg                                                                     >  8:150906/1‑149 (MQ=255)
                                                                                     aTCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAgg                                                                                                                            >  3:52201/1‑89 (MQ=255)
                                                                                     aTCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAgg                                                                                                                            <  4:52201/89‑1 (MQ=255)
                                                                                      tCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGa                                                              <  5:162816/149‑1 (MQ=255)
                                                                                        ccTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACt                                                            >  3:367191/1‑149 (MQ=255)
                                                                                        ccTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCCGCGACGGTGCTACCAGGCTCGACt                                                            >  2:111865/1‑149 (MQ=255)
                                                                                         cTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTg                                                           >  5:245475/1‑149 (MQ=255)
                                                                                            aCGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTggg                                                         <  3:226815/148‑1 (MQ=255)
                                                                                            aCGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTgg                                                          <  3:123855/147‑1 (MQ=255)
                                                                                            aCGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTgg                                                          <  7:87706/147‑1 (MQ=255)
                                                                                              gCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGc                                                      <  6:245475/149‑1 (MQ=255)
                                                                                                  tACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCt                                                  >  5:157349/1‑149 (MQ=255)
                                                                                                       ggTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCtgtt                                               <  6:202614/147‑1 (MQ=255)
                                                                                                           gCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTtgg                                          <  4:188688/148‑1 (MQ=255)
                                                                                                            cACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTtggtg                                        <  3:188686/149‑1 (MQ=255)
                                                                                                             aCTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGCGTGCGGCTGTAGTtcgtgg                                       >  7:116815/1‑149 (MQ=255)
                                                                                                                 gTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTgttgtt                                            <  6:301139/140‑1 (MQ=255)
                                                                                                                 gTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTgttgtt                                            >  5:301139/1‑140 (MQ=255)
                                                                                                                         gACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGaaaa                           >  5:334608/1‑149 (MQ=255)
                                                                                                                          aCAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGaaaa                           <  1:57900/148‑1 (MQ=255)
                                                                                                                            aGAATACGTCTTACCAGCGGCAATCTGATTGCCGTAGGTGGTGGCCGATGTCAACAGGCAGCCAGACATCAATACCTTGCCGTCATCCGGAGAAGGTGCTACCAGTCTCGACTGGGTGCGGCTGTTGTTGGTTGCCTGGTCGAAAATGc                        <  2:117712/149‑1 (MQ=255)
                                                                                                                            aGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGc                        >  4:282575/1‑149 (MQ=255)
                                                                                                                                tACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTa                    <  7:150906/149‑1 (MQ=255)
                                                                                                                                 aCGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAg                   >  4:269331/1‑149 (MQ=255)
                                                                                                                                   ggCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTg                                   <  4:132201/131‑1 (MQ=255)
                                                                                                                                   ggCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTg                                   >  3:132201/1‑131 (MQ=255)
                                                                                                                                       cacctacGGCAATCAGTTTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGgtg             >  5:289213/7‑149 (MQ=255)
                                                                                                                                        acctacGGCAATCAGTTTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGAAGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGgtgt            >  2:285391/6‑149 (MQ=255)
                                                                                                                                        aCCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGgtg             <  1:362892/148‑1 (MQ=255)
                                                                                                                                          ctacGGCAATCAGTTTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCg          >  6:349303/4‑149 (MQ=255)
                                                                                                                                          ctacGGCAATCAGTTTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCg          >  8:129034/4‑149 (MQ=255)
                                                                                                                                             cGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTg       <  4:252279/149‑1 (MQ=255)
                                                                                                                                                cAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTc    >  3:194348/1‑149 (MQ=255)
                                                                                                                                                 aaaCTGATTGCCGTAGGTGGGGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGTGTTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTc    <  1:111865/149‑1 (MQ=255)
                                                                                                                                                  aaCTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTCGc  >  2:3402/1‑149 (MQ=255)
                                                                                                                                                   aCTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTCGc  >  6:73094/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGCGAGATCATCTTGCCGCCTTCGGTATTGACGTCGAACTTGCCCTTGCCGGCGGTTTCATTGGTGATGTCATTGCCGGCGGCATCCCTCACGCTGTACTCGGTGGCACTGGTGAAGGTGACAGAATACGGCTGACCAGCGGCAAACTGATTGCCGTAGGTGGTGGCCGAGGTCAACAGGCCGCCAGACAGCAATACCTTGCCGTCATCCGGCGACGGTGCTACCAGGCTCGACTGGGTGCGGCTGTTGTTGGTGGCCTGGTCGAAAATGCTGTAGCCGGTGTCGTTGGTCGC  >  NC_002947/4972734‑4973028

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