Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,268,081 G→C 8.6% intergenic (+33/+26) PP_2863 → / ← PP_2864 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,268,0810GC 8.6% 407.7 / 23.1 128intergenic (+33/+26)PP_2863/PP_2864membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (68/49);  new base C (3/8);  total (71/57)
Fisher's exact test for biased strand distribution p-value = 6.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 130 total aligned reads displayed.

GCACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCAC  >  NC_002947/3267934‑3268230
                                                                                                                                                    |                                                                                                                                                     
gCACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGcc                                                                                                                                                     >  1:3621256/1‑150 (MQ=255)
gCACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGcc                                                                                                                                                     >  2:3983621/1‑150 (MQ=255)
gCACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGcc                                                                                                                                                     >  1:3045387/1‑150 (MQ=255)
gCACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGcc                                                                                                                                                     >  1:3044973/1‑150 (MQ=255)
gCACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGcc                                                                                                                                                     >  1:2013204/1‑150 (MQ=255)
  acaAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCcgc                                                                                                                                                   >  1:5114350/1‑150 (MQ=255)
  acaAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCAGGTAAAGCcgc                                                                                                                                                   <  1:2560132/150‑1 (MQ=255)
    taGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGa                                                                                                                                                 <  2:1926194/149‑1 (MQ=255)
     aGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGCTTTACCCGCGaa                                                                                                                                                <  1:4033740/150‑1 (MQ=255)
     aGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGCTTTACCCGCGaa                                                                                                                                                <  1:188671/150‑1 (MQ=255)
       aCCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAg                                                                                                                                              <  1:4819935/150‑1 (MQ=255)
        cccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAgg                                                                                                                                             <  2:238721/150‑1 (MQ=255)
        cccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGCTTTACCCGCGaagag                                                                                                                                             <  1:1722385/150‑4 (MQ=255)
         ccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAggg                                                                                                                                            <  2:3403702/150‑1 (MQ=255)
          cGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGc                                                                                                                                           >  1:4534680/1‑150 (MQ=255)
           gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGc                                                                                                                                           >  1:1566672/1‑150 (MQ=255)
           gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTTCCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGcc                                                                                                                                          <  2:3603696/150‑1 (MQ=38)
           gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGcc                                                                                                                                          <  2:2233696/150‑1 (MQ=255)
           gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGcc                                                                                                                                          <  2:3596241/150‑1 (MQ=255)
           gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGcc                                                                                                                                          >  2:721507/1‑150 (MQ=255)
           gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGAAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGcc                                                                                                                                          <  2:3086413/150‑1 (MQ=255)
             aaGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCgg                                                                                                                                        <  1:3122506/150‑1 (MQ=255)
                  ttCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTTCCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGTACaa                                                                                                                                   <  2:59277/150‑2 (MQ=37)
                  ttCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAg                                                                                                                                   <  1:230769/150‑1 (MQ=255)
                  ttCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGTACaa                                                                                                                                   <  2:124504/150‑2 (MQ=38)
                      ttCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCg                                                                                                                               >  2:3869792/1‑150 (MQ=255)
                        ccAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGAt                                                                                                                             <  1:116537/150‑1 (MQ=255)
                          aGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATca                                                                                                                           <  1:4770419/150‑1 (MQ=255)
                            gtggtgTGCTTGGGTTTTGTGCTGGTGCTCCTGGTTTTCTGGGCTGACTGGCTTTTTCTTGACGGTCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATcaca                                                                                                                         <  2:217261/150‑1 (MQ=255)
                            gtggtgTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATcaca                                                                                                                         <  2:3478277/150‑1 (MQ=255)
                             tggtgTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACAt                                                                                                                        <  1:4387986/150‑1 (MQ=255)
                               gtgtGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGa                                                                                                                      <  2:2661683/150‑1 (MQ=255)
                                tgtgCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAc                                                                                                                     <  1:3778738/150‑1 (MQ=255)
                                   gCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCaa                                                                                                                  <  1:4817950/150‑1 (MQ=255)
                                    cTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAAc                                                                                                                 <  2:702084/150‑1 (MQ=255)
                                     tGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAAcc                                                                                                                >  2:4586641/1‑150 (MQ=255)
                                                  ctgctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAAcc                                                                                                    >  2:3721520/1‑150 (MQ=255)
                                                   tgctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCg                                                                                                   >  2:2299675/1‑150 (MQ=255)
                                                    gctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTAGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGc                                                                                                  <  2:34852/150‑1 (MQ=255)
                                                      tgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCtt                                                                                                <  1:665100/150‑1 (MQ=255)
                                                         aCCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTggg                                                                                             <  1:2106351/150‑1 (MQ=255)
                                                              gtgtTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAAcc                                                                                        >  2:2230831/1‑150 (MQ=255)
                                                                gtTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGa                                                                                      <  1:3439247/150‑1 (MQ=255)
                                                                gtTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGa                                                                                      <  2:2654307/150‑1 (MQ=255)
                                                                 ttCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGAc                                                                                     >  1:4281138/1‑150 (MQ=255)
                                                                 ttCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGAc                                                                                     >  1:893561/1‑150 (MQ=255)
                                                                       gCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTa                                                                               >  1:3811795/1‑150 (MQ=255)
                                                                         aGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTAcc                                                                             >  2:1484017/1‑150 (MQ=255)
                                                                           aCTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAgg                                                                                                                                             >  2:3270116/1‑84 (MQ=255)
                                                                           aCTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAgg                                                                                                                                             <  1:3270116/84‑1 (MQ=255)
                                                                           aCTGGCTGTTACTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAg                                                                           >  2:405785/1‑150 (MQ=255)
                                                                            cTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTTGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGc                                                                          >  2:4691581/1‑150 (MQ=255)
                                                                             tGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCt                                                                         >  2:2890985/1‑150 (MQ=255)
                                                                             tGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCt                                                                         >  1:6845/1‑150 (MQ=255)
                                                                               gCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTgg                                                                       <  1:721507/150‑1 (MQ=255)
                                                                                     tCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGAc                                                                 >  2:4844049/1‑150 (MQ=255)
                                                                                       ttGACGGGCGTTATCTGGCGCTTGGTTGATGTGCAGTTGTTCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGTTCACATGACCAACCCCACCTGCAACGGCTTGGGCAACCGACGCACTACCAGCTGGTGCGAGCg                                                               >  2:2885503/1‑150 (MQ=255)
                                                                                         gACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCt                                                             <  2:49472/150‑1 (MQ=255)
                                                                                          aCGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTg                                                            >  2:1029414/1‑150 (MQ=255)
                                                                                                   aTCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTc                                                   >  1:3335939/1‑150 (MQ=255)
                                                                                                      tGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTCCCAGCTGGTGCGACCGCTGCAACACCTcaca                                                >  2:2573906/1‑150 (MQ=255)
                                                                                                        ccccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTcacaca                                              >  2:1346165/1‑150 (MQ=255)
                                                                                                         cccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAg                                             >  1:1227023/1‑150 (MQ=255)
                                                                                                          ccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAAc                                                                                                     >  2:2554530/1‑93 (MQ=255)
                                                                                                          ccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAAc                                                                                                     <  1:2554530/93‑1 (MQ=255)
                                                                                                          ccTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGt                                            <  2:2013204/150‑1 (MQ=255)
                                                                                                            ttGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCCACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTtc                                          >  2:4131343/1‑150 (MQ=255)
                                                                                                                ttGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTtcatca                                      >  1:3562062/1‑150 (MQ=255)
                                                                                                                  gATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTCCCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACg                                    >  2:4584230/1‑150 (MQ=255)
                                                                                                                    tgtgCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGcc                                  >  1:3877542/1‑150 (MQ=255)
                                                                                                                      tgCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGcccc                                >  2:2547263/1‑150 (MQ=255)
                                                                                                                      tgCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGcccc                                >  1:1903497/1‑150 (MQ=255)
                                                                                                                       gCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCa                               >  2:3527696/1‑150 (MQ=255)
                                                                                                                       gCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCa                               >  1:3175851/1‑150 (MQ=255)
                                                                                                                        cAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCAt                              >  2:2646059/1‑150 (MQ=255)
                                                                                                                         aGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATg                             >  1:1764038/1‑150 (MQ=255)
                                                                                                                         aGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATg                             <  2:3027291/150‑1 (MQ=255)
                                                                                                                          gTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                            <  1:3112313/150‑1 (MQ=255)
                                                                                                                            tGTTCCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCGTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgggg                          >  2:4241461/1‑150 (MQ=37)
                                                                                                                               aCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAg                       <  1:335808/150‑1 (MQ=255)
                                                                                                                               aCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAg                       <  2:468510/150‑1 (MQ=255)
                                                                                                                                ccGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCCCTACCAGCTGGTGCGACCGCTGCAACACCGCACACAGTTCATCACGCCCCAAGGGGTAgg                      >  1:176553/1‑150 (MQ=255)
                                                                                                                                ccGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAgg                      >  2:2941860/1‑150 (MQ=255)
                                                                                                                                ccGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCAGGGGGTAgg                      >  1:176527/1‑150 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAggg                     >  2:5018401/1‑150 (MQ=255)
                                                                                                                                 cGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAggg                     >  1:5003933/1‑150 (MQ=255)
                                                                                                                                   gccctttCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcg                   >  1:498342/6‑150 (MQ=255)
                                                                                                                                   gCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcg                   >  1:1506858/1‑150 (MQ=255)
                                                                                                                                   gCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcg                   >  1:2959572/1‑150 (MQ=255)
                                                                                                                                    ccctttCGCGGCTTTACCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcgg                  >  1:1235869/5‑150 (MQ=255)
                                                                                                                                     ctctTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcggc                 <  1:3869792/150‑1 (MQ=255)
                                                                                                                                      tctTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGgcggcg                >  2:3589720/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGcc              >  1:1149015/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGcc              >  2:4350510/1‑150 (MQ=255)
                                                                                                                                          cgcgGGTACAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCAt            <  1:3204746/150‑1 (MQ=255)
                                                                                                                                          cgcgGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCAt            <  1:3614698/150‑1 (MQ=255)
                                                                                                                                          cgcgGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCAt            <  1:2210421/150‑1 (MQ=255)
                                                                                                                                            cgGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGc          >  1:2266930/1‑150 (MQ=255)
                                                                                                                                               gTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGa       <  2:1741267/150‑1 (MQ=255)
                                                                                                                                                 aaaGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGAcc     <  2:1566672/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCACAAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCAC  >  NC_002947/3267934‑3268230

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