Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,022,851 C→G 5.4% intergenic (‑502/+166) PP_3547 ← / ← emrB NAD(P)‑binding oxidoreductase/multidrug transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,022,8510CG 5.4% 208.1 / 2.3 66intergenic (‑502/+166)PP_3547/emrBNAD(P)‑binding oxidoreductase/multidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base G (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGG  >  NC_002947/4022708‑4022997
                                                                                                                                                 |                                                                                                                                                  
aaccgcgACAAAAATTAACCAAAAACCCCTGGCTCCTATCATCCGGCCCCAGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGccc                                                                                                                                              <  1:207048/147‑1 (MQ=14)
  gctcttgctttggcttttgcttctaagcgcgcggtagttcaggcgacgcagattgcgacttcaggaggccgagcgcagggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGt                                                                                                                                            <  1:1022003‑M2/65‑1 (MQ=255)
      gACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTCaag                                                                                                                                        >  2:691967/1‑148 (MQ=255)
      gACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTCaag                                                                                                                                        >  2:691352/1‑148 (MQ=255)
           aaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATGGATGCCCGTCaagcgctg                                                                                                                                   <  1:691352/150‑8 (MQ=14)
           aaaGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCAAAAGGGAGGACCGCGCTGCGCCTACCCGGCCGGGTCCCACCGGTTGACGGGCATGGATGCCCGTCaagcgctg                                                                                                                                   <  1:691967/150‑8 (MQ=9)
             tgattttgcttctaagcgcgcgatagcccaggcgccgcggattgcgacttcaggaggccgagcgcagggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                                                 <  1:1324047‑M2/76‑1 (MQ=255)
             aGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                                                 <  1:884266/150‑1 (MQ=255)
                   cAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGc                                                                                                                           >  2:154794/1‑150 (MQ=255)
                    aaaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGcc                                                                                                                          >  2:931901/1‑150 (MQ=255)
                     aaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGccc                                                                                                                         <  1:317971/150‑1 (MQ=255)
                      ttgtaagcgcgcggtagttcaggcgacgcagattgcgacttcagaaggccgagcgcagggcttgcTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGccc                                                                                                                         <  1:945986‑M2/85‑1 (MQ=255)
                      ttctaagcgcgcggtagttcaggcgatgcagattgcgacttcaggaggccgagcgcagggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCa                                                                                                                                          <  2:1635075‑M2/67‑1 (MQ=255)
                      ttctaagcgcgcggtagttcaggcgatgcagattgcgacttcaggaggccgagcgcagggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCa                                                                                                                                          >  1:1635075‑M2/66‑132 (MQ=255)
                      aaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCt                                                                                                                        <  1:1432664/150‑1 (MQ=255)
                                 tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                           <  2:1070429/120‑1 (MQ=38)
                                 tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                                           >  1:1070429/1‑120 (MQ=37)
                                    tATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGcctcct                                                                                                          <  1:123098/150‑1 (MQ=37)
                                      tcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGa                                                                                                        <  2:1097174/150‑1 (MQ=34)
                                                 cTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATcc                                                                                             >  1:242344/1‑150 (MQ=35)
                                                          gCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTc                                                                                                    <  2:743007/134‑1 (MQ=37)
                                                          gCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTc                                                                                                    >  1:743007/1‑134 (MQ=37)
                                                            tCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTgg                                                                                  <  2:1330022/150‑1 (MQ=40)
                                                                 ggTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTAt                                                                             <  2:223650/150‑1 (MQ=255)
                                                                  gTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGcc                                                                                                                          >  2:423843/1‑104 (MQ=38)
                                                                  gTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGcc                                                                                                                          <  1:423843/104‑1 (MQ=38)
                                                                       cTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATcgcgcg                                                                       >  1:874172/1‑150 (MQ=255)
                                                                         cGCTCCGGGCTTGCTCCCGGGAGGACCGCGCCGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCt                                                                     <  1:1535481/150‑1 (MQ=255)
                                                                              cGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGaa                                                                >  1:1097745/1‑150 (MQ=255)
                                                                               gggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAg                                                               >  2:31204‑M2/9‑150 (MQ=255)
                                                                               gggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAg                                                               >  1:555483‑M2/9‑150 (MQ=255)
                                                                                    tgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCaaaa                                                          <  1:931901‑M2/147‑1 (MQ=255)
                                                                                       tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCa                                                       >  2:866422/1‑150 (MQ=255)
                                                                                       tCCCGGGAGGACCGCGCTGCAGACCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCa                                                       >  2:409565/1‑150 (MQ=255)
                                                                                         ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCaaaa                                                          >  1:671293/1‑145 (MQ=255)
                                                                                         ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCaaaa                                                          <  2:671293/145‑1 (MQ=255)
                                                                                         ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaa                                                     >  1:637555/1‑150 (MQ=255)
                                                                                         ccGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCCGCACCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaa                                                     >  1:85532/1‑150 (MQ=255)
                                                                                          cGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaaa                                                    <  2:555483/150‑1 (MQ=255)
                                                                                          cGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCACGCTTAGAAGCAAAATCaaaa                                                    <  1:154794/150‑1 (MQ=255)
                                                                                            ggaggaCCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGc                                                  >  2:1270251/1‑150 (MQ=255)
                                                                                                 aCCGCGCTGCAGGCCCCATCCTGGGGCACCGAGTACGCCATGCATGCATTGCCGTCACCTCCCTTCGCAAGACCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAag                                             <  1:365268/150‑1 (MQ=18)
                                                                                                   cgcgCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcg                                           >  2:58852/1‑150 (MQ=255)
                                                                                                             gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagag                                 >  1:1387703/1‑150 (MQ=255)
                                                                                                             gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagag                                 >  1:1436763/1‑150 (MQ=255)
                                                                                                             gCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATCACCGTCATCTTGGTCTCCAAGCCCTGCGCGCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagag                                 >  1:850676/1‑150 (MQ=18)
                                                                                                              ccccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTTGTGGTCGGTGACCAAGGTGCGAACGCGCGCGTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagc                                >  1:603580/1‑150 (MQ=17)
                                                                                                              ccccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagcgcgagagc                                >  2:119392/1‑150 (MQ=255)
                                                                                                               cccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGGCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcg                               >  1:334267/1‑150 (MQ=255)
                                                                                                                ccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCcgcg                                                                                          >  2:1495764/1‑90 (MQ=40)
                                                                                                                ccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCcgcg                                                                                          <  1:1495764/90‑1 (MQ=39)
                                                                                                                         gggCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgaga                      <  2:850676/150‑1 (MQ=255)
                                                                                                                            cccAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcg                   >  1:663396/1‑150 (MQ=255)
                                                                                                                             ccAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcga                  <  2:1097745/150‑1 (MQ=255)
                                                                                                                               cAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgag                 <  2:54015/150‑1 (MQ=255)
                                                                                                                                aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgag                 >  1:1397188/1‑149 (MQ=255)
                                                                                                                                aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgag                 <  2:1397188/149‑1 (MQ=255)
                                                                                                                                  cgcTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgagagc              <  2:637555/150‑1 (MQ=255)
                                                                                                                                    cTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAagcgagagcgagagcgagagcgagagcgagagcga            <  2:242344/150‑1 (MQ=255)
                                                                                                                                      tGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAGAATCAAAGGCCGAAgcgagagcgagagcgagagcgagagcgagagcgaga          <  1:31204/150‑1 (MQ=255)
                                                                                                                                      tGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAgcgagagcgagagcgagagcgagagcgagagcgaga          <  1:1270251/150‑1 (MQ=255)
                                                                                                                                            gCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaaa                                                    >  2:1568051/1‑102 (MQ=255)
                                                                                                                                            gCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCaaaa                                                    <  1:1568051/102‑1 (MQ=255)
                                                                                                                                            gCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCg    >  2:458943/1‑150 (MQ=255)
                                                                                                                                             cATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGACGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCgg   <  2:334267/150‑1 (MQ=255)
                                                                                                                                              aTCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCggg  >  2:1410814/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCCGCGGCGCCTGGGCTATCGCGCGCTTAGAAGCAAAATCAAAAGCCGAAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCGAGAGCATCGGG  >  NC_002947/4022708‑4022997

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: