Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,648,290 G→A 6.5% intergenic (+75/‑48) dbpA → / → yhiN 23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,648,2900GA 6.5% 671.8 / 23.2 217intergenic (+75/‑48)dbpA/yhiN23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (125/76);  new base A (12/2);  total (139/78)
Fisher's exact test for biased strand distribution p-value = 9.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

TGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACC  >  NC_002947/5648142‑5648437
                                                                                                                                                    |                                                                                                                                                   
tGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGa                                                                                                                                                    <  1:3190032/150‑1 (MQ=255)
  gCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAc                                                                                                                                                  <  1:4715661/150‑1 (MQ=255)
  gCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAc                                                                                                                                                  >  1:2269867/1‑150 (MQ=255)
    cAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAccc                                                                                                                                                >  1:768021/1‑150 (MQ=255)
    cAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAccc                                                                                                                                                <  1:1502888/150‑1 (MQ=255)
    cAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAccc                                                                                                                                                >  1:332739/1‑150 (MQ=255)
    cAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAccc                                                                                                                                                >  1:2694304/1‑150 (MQ=255)
     aaGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCg                                                                                                                                               >  2:226921/1‑150 (MQ=255)
     aaGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCg                                                                                                                                               >  2:5600450/1‑150 (MQ=255)
      aGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGc                                                                                                                                              >  2:2161160/1‑150 (MQ=255)
      aGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGc                                                                                                                                              >  2:2160949/1‑150 (MQ=255)
       gCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCt                                                                                                                                             >  2:1285279/1‑150 (MQ=255)
                gCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAggg                                                                                                                                    <  2:3294253/150‑1 (MQ=255)
                gCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAggg                                                                                                                                    >  2:3175811/1‑150 (MQ=255)
                  aGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGtg                                                                                                                                  >  1:1053835/1‑149 (MQ=255)
                    cGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTACAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCa                                                                                                                                <  2:2946234/150‑1 (MQ=255)
                     ggCTGAACAGCGGCAAGATCACGGGTCGCGCCCTGAAAGTCCGCATCATCTAAAGCCAGGCAACGCCCCCTGTCGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCAGCGCATTCGGGGGTGAACCCGCTCCCACAGGGTTCAg                                                                                                                               <  1:5345311/150‑1 (MQ=38)
                        tGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAt                                                                                                                            >  2:3028241/1‑150 (MQ=255)
                             aGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTCGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaa                                                                                                                       <  2:4797742/150‑1 (MQ=255)
                              gCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                      <  2:3276089/150‑1 (MQ=255)
                              gCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCAGCTCCCACAGGGTTCAGCATTTCaaa                                                                                                                      >  2:4287641/1‑150 (MQ=255)
                                 gaaaGATCACGGGTCGCGCCCTTAAAGTCAGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCt                                                                                                                   <  1:1119721/148‑1 (MQ=255)
                                  cAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTg                                                                                                                  >  2:5398442/1‑150 (MQ=255)
                                    aGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGaa                                                                                                                <  2:336277/150‑1 (MQ=255)
                                                  gcCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGc                                                                                                  >  2:3322854/1‑150 (MQ=255)
                                                   cccTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCa                                                                                                 >  1:1036730/1‑150 (MQ=255)
                                                     cTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACt                                                                                               >  1:3888502/1‑150 (MQ=255)
                                                     cTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACt                                                                                               >  2:3387985/1‑150 (MQ=255)
                                                               gCATTATCTACCCCCAGGCAAAGCCCCCTGAAGGTGCGGGTTTACCCGCGAATGCGGCCGTGAATGCACAATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGt                                                                                     <  1:525804/150‑1 (MQ=38)
                                                               gCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGt                                                                                     >  1:3056888/1‑150 (MQ=255)
                                                                     tCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCGTTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGatcat                                                                               >  1:287169/1‑150 (MQ=255)
                                                                      cTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGCAGGTTCATACCGTGCACTCCACCGACGTGatcatc                                                                              >  2:3447117/1‑150 (MQ=255)
                                                                      cTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGatcatc                                                                              >  2:3003354/1‑150 (MQ=255)
                                                                      cTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGatcatc                                                                              >  2:1239340/1‑150 (MQ=255)
                                                                      cTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGatcatc                                                                              >  2:2352396/1‑150 (MQ=255)
                                                                       tAAAGCCAGGCAACGCCCCCGGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAAGGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATcc                                                                             <  2:5191302/150‑1 (MQ=255)
                                                                       tAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCACATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATcc                                                                             >  2:3764804/1‑150 (MQ=255)
                                                                        aaaGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCt                                                                            >  1:1672655/1‑150 (MQ=255)
                                                                          aGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTg                                                                          <  1:1243722/150‑1 (MQ=255)
                                                                             cAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcg                                                                       >  1:5556612/1‑150 (MQ=255)
                                                                               ggCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgcc                                                                     >  1:937346/1‑150 (MQ=255)
                                                                                gCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccg                                                                    <  1:5102860/150‑1 (MQ=255)
                                                                                  aaaGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCTGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccggc                                                                  <  1:1644562/150‑1 (MQ=255)
                                                                                    tgaacccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcgccggcgc                                                                >  1:5051617/5‑150 (MQ=17)
                                                                                             tAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAttt                                                                                                                          <  2:5078949/83‑1 (MQ=255)
                                                                                             tAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCt                                                       <  1:4689880/150‑1 (MQ=255)
                                                                                               ggAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGa                                                     <  2:2392269/150‑1 (MQ=255)
                                                                                               ggAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGa                                                     <  1:4546502/150‑1 (MQ=255)
                                                                                                gAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGACCCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTTATCATCCTTGGCGCCGGCGCCGCCGGCCTGAt                                                    <  2:2202939/150‑1 (MQ=255)
                                                                                                 aGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtg                                                   >  2:2565571/1‑150 (MQ=255)
                                                                                                 aGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtg                                                   >  2:4081935/1‑150 (MQ=255)
                                                                                                    gggTTTACCCGCGAATGCGGCAGTGAATGCACCAGCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTgc                                                <  1:1910483/150‑1 (MQ=255)
                                                                                                    gggTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTgc                                                >  1:1335597/1‑150 (MQ=21)
                                                                                                        aaaCCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGccc                                            >  2:3594797/3‑150 (MQ=255)
                                                                                                          aCCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGcc                                             >  1:5589993/1‑147 (MQ=255)
                                                                                                            ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  2:1690939/1‑150 (MQ=255)
                                                                                                            ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        <  1:1691203/150‑1 (MQ=255)
                                                                                                            ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  1:606361/1‑150 (MQ=255)
                                                                                                            ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  2:5505576/1‑150 (MQ=255)
                                                                                                            ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  1:4073719/1‑150 (MQ=255)
                                                                                                            ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCACCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  1:1947180/1‑150 (MQ=255)
                                                                                                            ccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTAAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  2:4276527/1‑150 (MQ=255)
                                                                                                            ccGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  1:3995037/1‑150 (MQ=37)
                                                                                                            ccGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  1:1891002/1‑150 (MQ=25)
                                                                                                            ccGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                                        >  1:3467234/1‑150 (MQ=25)
                                                                                                             cgcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTc                                       >  2:803312/1‑150 (MQ=255)
                                                                                                              gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                                      >  1:4122278/1‑150 (MQ=255)
                                                                                                              gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                                      <  2:745477/150‑1 (MQ=255)
                                                                                                              gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                                      <  2:3996407/150‑1 (MQ=255)
                                                                                                              gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                                      >  2:1199106/1‑150 (MQ=255)
                                                                                                              gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                                      >  1:5599412/1‑150 (MQ=255)
                                                                                                              gcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                                      >  1:2926875/1‑150 (MQ=25)
                                                                                                                gAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGt                                    >  1:4218876/1‑150 (MQ=255)
                                                                                                                gAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGt                                    >  1:4219020/1‑150 (MQ=255)
                                                                                                                   tGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCCAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGcc                                 <  1:823490/150‑1 (MQ=255)
                                                                                                                     cGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCg                               >  2:1792830/1‑150 (MQ=255)
                                                                                                                     cGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCg                               >  1:2079077/1‑150 (MQ=255)
                                                                                                                        cAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAg                            >  1:1916024/1‑150 (MQ=255)
                                                                                                                          gTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCg                          >  2:2219693/1‑150 (MQ=255)
                                                                                                                          gTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCg                          >  1:3707509/1‑150 (MQ=255)
                                                                                                                          gTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCg                          >  1:2576803/1‑150 (MQ=255)
                                                                                                                              aTGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGgccg                      >  2:415898/1‑150 (MQ=255)
                                                                                                                               tGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGgccgc                     <  2:4308063/150‑1 (MQ=255)
                                                                                                                                gCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGgccgcc                    >  2:1779295/1‑150 (MQ=255)
                                                                                                                                  aCCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCgg                  >  2:1696492/1‑150 (MQ=255)
                                                                                                                                    cATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGt                >  1:4447560/1‑150 (MQ=255)
                                                                                                                                     gccGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtg               <  1:419284/148‑1 (MQ=255)
                                                                                                                                     aTCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACACCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtg               >  2:5198703/1‑150 (MQ=255)
                                                                                                                                      tCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtgc              >  2:2830445/1‑150 (MQ=255)
                                                                                                                                      tCGCAGTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTGGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtgc              <  1:3653132/150‑1 (MQ=255)
                                                                                                                                        gCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtgctg            <  1:3028241/150‑1 (MQ=255)
                                                                                                                                         cATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtgctgc           >  2:1494489/1‑150 (MQ=255)
                                                                                                                                          gggagcgGGTTTACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGtgctgct          >  2:297970/5‑150 (MQ=255)
                                                                                                                                              gcgGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTc      <  1:1704868/150‑1 (MQ=255)
                                                                                                                                              gcgGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTc      >  2:4806705/1‑150 (MQ=255)
                                                                                                                                              gcgGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTc      >  2:2151642/1‑150 (MQ=255)
                                                                                                                                              gcgGGTAAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTc      >  2:553980/1‑150 (MQ=255)
                                                                                                                                                gggTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGa    >  1:592273/1‑150 (MQ=255)
                                                                                                                                                gggTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGa    >  2:2919849/1‑150 (MQ=255)
                                                                                                                                                 ggTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGAc   <  2:1672654/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCGGGTGCTGCTGCTCGACC  >  NC_002947/5648142‑5648437

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