Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,341,183 A→G 7.1% intergenic (‑220/+23) PP_2935 ← / ← ybiT membrane protein/ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,341,1830AG 7.1% 481.1 / 18.0 155intergenic (‑220/+23)PP_2935/ybiTmembrane protein/ABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base A (85/58);  new base G (5/6);  total (90/65)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 160 total aligned reads displayed.

CTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCTGA  >  NC_002947/3341034‑3341330
                                                                                                                                                     |                                                                                                                                                   
cTAACGAATTATTCGTCGGGCCGTGCCGGTCGATGGGATGGGGATGTACTCTTCGCGGGCTGGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCga                                                                                                                                                     <  2:2634451/150‑1 (MQ=255)
   aCGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGAt                                                                                                                                                  <  2:2841933/150‑1 (MQ=255)
    cGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATg                                                                                                                                                 >  2:373268/1‑150 (MQ=255)
      aaTTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGc                                                                                                                                               >  2:1107363/1‑150 (MQ=255)
        ttattCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTTCCCCCGCTCCCAAGGCCCCCCAACCGTCCAAAACCAGTGCCGGCCCAGTGGGGGCGGGCCACCCCCCCAAGGGGCCCCGCCCACCCAGATTGccc                                                                                                                                             >  1:3485698/1‑149 (MQ=255)
         tattCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAAAACCAGTGCGGTCCCTGTGGGAGCGGGCCAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAg                                                                                                                                            >  2:2401230/1‑150 (MQ=255)
          attCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAgg                                                                                                                                           <  2:717384/150‑1 (MQ=255)
           ttCGTCGGGCCGTGCCGGTCGGTGGGTTGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACAAGCGTCCAATACGTGTGCGGTCCCTGTGGGAGGGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                                                                          <  1:4363407/150‑1 (MQ=255)
           ttCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                                                                          <  2:1482957/150‑1 (MQ=39)
           tgcgtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                                                                          <  1:1816400/148‑1 (MQ=255)
               tcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTTCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATACGACATGGCCAGGCCCTc                                                                                                                                      >  2:1720048/1‑150 (MQ=255)
                  ggCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGGTTTCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTTTGCGGTACCTGTTGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGTTGGCCAGGCCGTCCCa                                                                                                                                   >  2:3224022/1‑150 (MQ=39)
                   gCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCGGCTCTCCCAAGGCCCCACAAGGGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCat                                                                                                                                  <  1:2931164/150‑1 (MQ=255)
                      gTGCCGGTCGGTGGGCTGGGGAGGTACTCTGCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGGCCAAGACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatcaa                                                                                                                               <  2:3035596/150‑1 (MQ=39)
                      gTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCGCACTCCCAAGGCCCCCCAAGCGTCCAAAACCAGTGCGGTCCCGGTGGGGGCGGGCAAGCCCCCGAACAGGCCCGGGCATACGAGATGGCCCGGCCGGCCCagcaa                                                                                                                               >  1:375488/1‑150 (MQ=255)
                       tGCCGGTCGGTGGGTTGGGGATGTACTCTTCGCGGGCTGGCCCGCTCTCCCAAGCCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatcaat                                                                                                                              <  1:1386449/150‑1 (MQ=255)
                           ggtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACCAGCGTCCAAGACCGGTGCGGGCCCTGTGGGCGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                                                                          <  2:3099537/150‑1 (MQ=255)
                           ggggggggggttggggatgTACTGTTCGCGGGCTTGCCCGCACTCCCACGGCCCCACACGCGTCCACTACCTGTGCGGGCCGGGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                                                                          <  1:3243261/143‑1 (MQ=255)
                             tcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCCTCAATCAGAcc                                                                                                                        >  2:3209857/1‑150 (MQ=255)
                                         ggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGCAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTgg                                                                                                            >  2:2166826/1‑150 (MQ=255)
                                          gatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATAAGACCAGCACACCCTGGc                                                                                                           >  2:1760556/1‑150 (MQ=255)
                                          gatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCCTCCCATCAATCAGACCAGCACACCCAGGc                                                                                                           >  2:2393307/1‑150 (MQ=255)
                                             gTACTCTTCGCGGGCGTGCCCCCACTCCCAACGCCCCCCAACCGTCCAAAACCTCTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCCTCCGAGATGGCCCGGCCGTCCCCTCAAACAGACCCGCCCCCCCGGGCTAc                                                                                                        >  2:2346978/1‑150 (MQ=25)
                                              tACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCCTCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCAGGCTACg                                                                                                       >  2:2937927/1‑150 (MQ=255)
                                               aCTCTTCGCGGGCTTGCCCGCTCCCCCAAGGCCCCACAAGCGTCCAATACGTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGc                                                                                                      <  1:3659181/150‑1 (MQ=255)
                                                 gctggGCGGGCTGGCCCGCACCCCCACGGCCCCACACGCGCCCATTCCCTGGGCGGTCCGGGTGGGAGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg                                                                                                    >  1:4198306‑M1/6‑150 (MQ=255)
                                                 cCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCGGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg                                                                                                    >  2:882787‑M1/2‑150 (MQ=255)
                                                 TCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACATGCGTCCATTACCTTTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg                                                                                                    >  2:1518543‑M1/1‑150 (MQ=255)
                                                 TCTGCGGGGGCTGGCCCGCTCTCCCAAGGCCCCACACGCGTCCAATACCTGTGCGGTCCCGGGGGGCGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg                                                                                                    >  1:278639‑M1/1‑150 (MQ=255)
                                                   agaGCGGGCAAGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGCGCTGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGt                                                                                                  <  2:4165795‑M1/147‑1 (MQ=255)
                                                    TCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                                                                                 >  2:4520998‑M1/1‑150 (MQ=255)
                                                     CGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAg                                                                                                >  2:3844774‑M1/1‑150 (MQ=255)
                                                          GCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATc                                                                                           <  2:3591613‑M1/150‑1 (MQ=255)
                                                            ttgcccgctcccacagggaccgcacaggtattggacgcttgtggggccttgggagagcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGt                                                                                         >  1:3024314‑M2/77‑150 (MQ=255)
                                                            TTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGt                                                                                         <  1:3502920‑M1/150‑1 (MQ=255)
                                                            TTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGt                                                                                         <  2:1961573‑M1/150‑1 (MQ=255)
                                                             tgcccgctcccacagggaccgcacaggtattggacgcttgtggggccttgggagagcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa                                                                                        >  1:2498630‑M2/76‑150 (MQ=255)
                                                             gGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa                                                                                        >  2:1339988‑M1/2‑150 (MQ=255)
                                                             gGCCCGCTCCCCCACGGCCCCACACGCGTCCAACCCGTGTGCGGTCCCTGGGGGCGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa                                                                                        >  1:4205979‑M1/2‑150 (MQ=255)
                                                             TGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa                                                                                        >  2:1339601‑M1/1‑150 (MQ=255)
                                                              gcccgctcccacagggaccgcacaggtattggacgcttgtggggccttgggagagcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAg                                                                                       >  2:1822933‑M2/75‑150 (MQ=255)
                                                              GCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAg                                                                                       >  2:2962822‑M1/1‑150 (MQ=255)
                                                                CCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGt                                                                                     <  2:2595175‑M1/150‑1 (MQ=255)
                                                                CCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGt                                                                                     <  1:307565‑M1/150‑1 (MQ=255)
                                                                  GCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTAc                                                                                   <  2:633133‑M1/150‑1 (MQ=255)
                                                                   CTCTCCCAAGGCCCCACAAGCGTCCAATACGTGGGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTAcc                                                                                  >  1:2478806‑M1/1‑150 (MQ=255)
                                                                      tcCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCt                                                                               >  2:1055388/1‑150 (MQ=255)
                                                                        ccACGGCCCCACATGCGTCCAATCCCTGTGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGa                                                                             <  1:843556/150‑1 (MQ=255)
                                                                         catggCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                                                                            <  2:428437/147‑1 (MQ=255)
                                                                          ccgctcccacagggaccgcacaggtattggacgcttgtgggagcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAg                                                                           >  1:690250‑M2/63‑150 (MQ=255)
                                                                          aaGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAg                                                                           >  2:4408274/1‑150 (MQ=255)
                                                                             gCCCCACAAGCGGCCAATACCTCTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCCTCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCAACGCAGGTAGTCATCGTAGGTACCGCTGAAGTcc                                                                        >  1:3161623/1‑150 (MQ=255)
                                                                              ccccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCCGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTcca                                                                       >  1:4163546/1‑150 (MQ=255)
                                                                              ccccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTcca                                                                       <  2:57829/150‑1 (MQ=255)
                                                                                ccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacc                                                                     >  2:41577/1‑150 (MQ=255)
                                                                                 cacaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacca                                                                    >  1:79280/1‑150 (MQ=255)
                                                                                 cacaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacca                                                                    >  2:2447211/1‑150 (MQ=255)
                                                                                  acaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccac                                                                   >  1:3925101/1‑150 (MQ=255)
                                                                                    aaGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccacac                                                                 >  1:2263107/1‑150 (MQ=255)
                                                                                    aaGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccacac                                                                 >  2:3113901/1‑150 (MQ=255)
                                                                                        gTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATc                                                             <  2:3956835/150‑1 (MQ=255)
                                                                                        gTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATc                                                             <  1:1073188/150‑1 (MQ=255)
                                                                                        gTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATc                                                             <  2:3068690/150‑1 (MQ=255)
                                                                                           cAATACGTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGc                                                          <  1:2242778/150‑1 (MQ=255)
                                                                                            aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGcc                                                         <  1:4408274/150‑1 (MQ=255)
                                                                                            aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCAGCGGcc                                                         >  1:473924/1‑150 (MQ=255)
                                                                                                      cGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAg                                                                                       >  1:2601706/1‑110 (MQ=255)
                                                                                                      cGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAg                                                                                       <  2:2601706/110‑1 (MQ=255)
                                                                                                       ggTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCgat                                              <  1:2198174/150‑1 (MQ=255)
                                                                                                          cccTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCgatgat                                           <  2:4232638/150‑1 (MQ=255)
                                                                                                            cTGTGGGTGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATgc                                         <  1:645754/150‑1 (MQ=255)
                                                                                                            cTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGTCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATgc                                         >  1:2600363/1‑150 (MQ=255)
                                                                                                             tgtgGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATgcg                                        >  1:3543183/1‑150 (MQ=255)
                                                                                                             tgtgGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATgcg                                        >  2:1607306/1‑150 (MQ=255)
                                                                                                               tgGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGcgcg                                      >  2:850560/1‑150 (MQ=255)
                                                                                                               tgGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGcgcg                                      >  2:850574/1‑150 (MQ=255)
                                                                                                                 ttcgcgggcttgcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTg                                    >  1:1595212‑M2/24‑150 (MQ=255)
                                                                                                                 ggAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTg                                    >  2:3498389/1‑150 (MQ=255)
                                                                                                                   agcgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGc                                  <  2:2263107‑M2/129‑1 (MQ=255)
                                                                                                                     cgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCa                                <  2:2793199‑M2/131‑1 (MQ=255)
                                                                                                                     cgggcaagcccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCa                                >  2:334923‑M2/20‑150 (MQ=255)
                                                                                                                              ccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGaca                       >  1:323974‑M2/11‑150 (MQ=255)
                                                                                                                              ccgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGaca                       >  1:1479075‑M2/11‑150 (MQ=255)
                                                                                                                               cgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGacac                      <  1:2401244‑M2/141‑1 (MQ=255)
                                                                                                                               cgcgaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGacac                      >  2:296206‑M2/10‑150 (MQ=255)
                                                                                                                                  gaagagGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGaa                   >  2:2165470‑M2/7‑150 (MQ=255)
                                                                                                                                       gGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTcgc              <  1:933431‑M2/149‑1 (MQ=255)
                                                                                                                                       gGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTcgc              <  2:4133116‑M2/149‑1 (MQ=255)
                                                                                                                                        gCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTcgcg             <  1:1811360/150‑1 (MQ=255)
                                                                                                                                         cccGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCgg            <  1:1107363/150‑1 (MQ=255)
                                                                                                                                          ctggCCATCTCGGATGGCCGGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAg                               >  1:2958415/3‑130 (MQ=255)
                                                                                                                                          ctggCCATCTCGGATGGCCGGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGt           >  2:4512341/3‑150 (MQ=255)
                                                                                                                                          ctggCCATCTCGGATGGCCGGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGt           >  2:2281651/3‑150 (MQ=255)
                                                                                                                                           tggCCATCTCGGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTc          >  1:1392809/2‑150 (MQ=255)
                                                                                                                                             gCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGt        >  1:4059568/1‑150 (MQ=255)
                                                                                                                                               cATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTgg      >  1:1733167/1‑150 (MQ=255)
                                                                                                                                               cATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTgg      <  2:2677990/150‑1 (MQ=255)
                                                                                                                                                 tCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCt    <  2:1733167/150‑1 (MQ=255)
                                                                                                                                                 tCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCt    >  1:2260861/1‑150 (MQ=255)
                                                                                                                                                   cGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCTGa  <  1:1954482/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGACAGCTCGATGATGCGCGTGGCCAGCGACGACACGAACTCGCGGTCGTGGCTGA  >  NC_002947/3341034‑3341330

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