Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,786,064 T→A 5.9% intergenic (+38/‑148) rpsB → / → tsf 30S ribosomal protein S2/elongation factor Ts

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,786,0640TA 5.9% 534.6 / 13.6 162intergenic (+38/‑148)rpsB/tsf30S ribosomal protein S2/elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base T (80/70);  new base A (2/7);  total (85/77)
Fisher's exact test for biased strand distribution p-value = 9.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 174 total aligned reads displayed.

TATCGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785915‑1786213
                                                                                                                                                     |                                                                                                                                                      
tATCGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCt                                                                                                                                                        >  2:5100050/1‑150 (MQ=255)
 aTCGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCtt                                                                                                                                                       >  2:5000529/1‑150 (MQ=255)
  tCGAGCTGTACATGACTTCGAGGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCttt                                                                                                                                                      <  2:87635/150‑1 (MQ=255)
     aGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGcc                                                                                                                                                   <  1:2311483/150‑1 (MQ=255)
        tGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGccccc                                                                                                                                                >  2:3703477/1‑150 (MQ=255)
          tACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCtt                                                                                                                                             >  1:1770745/1‑150 (MQ=255)
          tACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCtt                                                                                                                                             <  2:1642247/150‑1 (MQ=255)
           aCATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGGTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttt                                                                                                                                            <  1:1958548/150‑1 (MQ=255)
            cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCtttt                                                                                                                                           <  2:1511489/150‑1 (MQ=255)
            cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCtttt                                                                                                                                           <  1:425938/150‑1 (MQ=255)
            cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTAGGCACGCAAAAAGGGGGCTTTGCCCCCtttt                                                                                                                                           >  2:3008586/1‑150 (MQ=255)
             aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttttt                                                                                                                                          <  1:2416008/150‑1 (MQ=255)
             aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                          <  2:2082373/150‑1 (MQ=255)
             aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGGTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                          <  2:120260/150‑1 (MQ=255)
               gACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGGTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGc                                                                                                                                        <  2:43610/150‑1 (MQ=255)
                 cTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGGTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCTTTTTGcgt                                                                                                                                      <  1:3004657/150‑3 (MQ=255)
                  ttCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCAc                                                                                                                                     >  2:400821/1‑150 (MQ=255)
                  ttCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCAc                                                                                                                                     >  2:4792620/1‑150 (MQ=255)
                  ttCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCAc                                                                                                                                     >  1:1860442/1‑150 (MQ=255)
                     gATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCtt                                                                                                                                  >  2:1267419/1‑150 (MQ=255)
                     gATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCtt                                                                                                                                  <  1:694960/150‑1 (MQ=255)
                         gCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAATGGGGGCTTTGCCCCCTTTTTGCCACCTTGaaa                                                                                                                              >  1:239650/1‑150 (MQ=255)
                            gACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATcc                                                                                                                           >  2:4951025/1‑150 (MQ=255)
                             aCGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCct                                                                                                                          <  1:2485169/150‑1 (MQ=255)
                              cGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctg                                                                                                                         <  1:291235/150‑1 (MQ=255)
                               gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgc                                                                                                                        >  1:2280134/1‑150 (MQ=255)
                               gCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCCTTGCCCCCATTTTTCCCACCTGGAAAGCctc                                                                                                                         >  1:2534553/1‑149 (MQ=255)
                                cAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                       <  2:1217765/150‑1 (MQ=255)
                                    cATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCa                                                                                                                   >  1:728375/1‑150 (MQ=255)
                                       ccGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCa                                                                                                                >  2:4260934/1‑150 (MQ=255)
                                         gcgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATc                                                                                                              >  2:1118270/1‑150 (MQ=255)
                                              cGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                         <  2:666350/150‑1 (MQ=255)
                                              cGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                         >  1:1914240/1‑150 (MQ=255)
                                              cGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                         <  2:3218506/150‑1 (MQ=255)
                                               gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGcccc                                                                                                        <  2:569329/150‑1 (MQ=255)
                                               gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGcccc                                                                                                        >  1:4655250/1‑150 (MQ=255)
                                                caacaaCGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCg                                                                                                       >  1:1479545/1‑150 (MQ=255)
                                                caacaaCGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATACTGCGGTCAGCATCGGCCCCg                                                                                                       >  1:1479821/1‑150 (MQ=255)
                                                 aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCCCCTTGAAATCCTGCTGTCAGCATCGGCCCCGc                                                                                                      >  1:223886/1‑150 (MQ=255)
                                                    aacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGcatc                                                                                                   >  1:4853572/1‑150 (MQ=255)
                                                      cGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGcatcat                                                                                                 <  1:3603100/150‑1 (MQ=255)
                                                       gTTGCCGGCGGCACCGAAGTTTACGCTGCAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATg                                                                                                <  1:1810968/150‑1 (MQ=255)
                                                        ttGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATgg                                                                                               <  1:3315384/150‑1 (MQ=255)
                                                        gtGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATgg                                                                                               <  2:3937282/149‑1 (MQ=255)
                                                            cggcggCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCct                                                                                           >  2:3477202/1‑150 (MQ=255)
                                                               cggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                        <  1:5021594/150‑1 (MQ=255)
                                                                  cACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaa                                                                                     >  1:2167205/1‑150 (MQ=255)
                                                                         gTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                              <  1:3329230/150‑1 (MQ=255)
                                                                         gTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                              <  1:376368/150‑1 (MQ=255)
                                                                             aCGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGc                                                                          <  2:4137189/150‑1 (MQ=255)
                                                                              cGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCg                                                                         <  1:4544459/150‑1 (MQ=255)
                                                                                 tGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAt                                                                      >  2:203373/1‑150 (MQ=255)
                                                                                 tGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAt                                                                      >  2:3984305/1‑150 (MQ=255)
                                                                                 tGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCGGTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCAGGGGCCTGCTGACATGAAAACGGCGGag                                                                      >  2:541603/1‑149 (MQ=255)
                                                                                   aagaagCGGCTGCACCTGCTGCTGAGTAATAGGCTCTAGCGTCTACTTGGCACGCAAATAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGTAATGAAAACAGCGGAttt                                                                    <  1:4066045/150‑1 (MQ=255)
                                                                                   aagaagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAttt                                                                    >  2:414219/1‑150 (MQ=255)
                                                                                     gaagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCGGTGCCCCCTTTTTGCCACCGTGAAATCCTGCGGTCAGCCTCGGCCCCGCATCATGGGCCTGCGGAAATGCACACAGCGGAGTTGc                                                                  >  2:3246175/1‑150 (MQ=255)
                                                                                              gCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAACAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGaa                                                         <  2:3108587/150‑1 (MQ=255)
                                                                                               cACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAAc                                                        >  1:3876647/1‑150 (MQ=255)
                                                                                               cACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAAc                                                        <  2:1033269/150‑1 (MQ=255)
                                                                                                  ctgctgctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGcc                                                     >  1:548408/1‑150 (MQ=255)
                                                                                                       gctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGa                                                >  1:4422089/1‑150 (MQ=255)
                                                                                                         tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAATCAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                              >  1:12204/1‑150 (MQ=255)
                                                                                                         tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                              >  2:1607857/1‑150 (MQ=255)
                                                                                                         tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                              >  1:3162528/1‑150 (MQ=255)
                                                                                                              aaTAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACgg                                         <  1:1387799/150‑1 (MQ=255)
                                                                                                               aTAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGc                                                                          <  2:1136085/116‑1 (MQ=255)
                                                                                                               aTAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGc                                                                          >  1:1136085/1‑116 (MQ=255)
                                                                                                               aTAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACggg                                        <  1:4481829/150‑1 (MQ=255)
                                                                                                                 aGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTg                                      <  1:2674652/150‑1 (MQ=255)
                                                                                                                        aGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGt                               <  1:1185719/150‑1 (MQ=255)
                                                                                                                         gCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGtt                              <  1:4644607/150‑1 (MQ=255)
                                                                                                                         gCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAAGCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGtt                              >  1:182223/1‑150 (MQ=255)
                                                                                                                          cGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTg                             <  1:89595/150‑1 (MQ=255)
                                                                                                                             cTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGaaa                          <  2:1248125/150‑1 (MQ=255)
                                                                                                                                 ttGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACt                      <  1:62277/150‑1 (MQ=255)
                                                                                                                                   ggCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACttt                    >  1:662083/1‑150 (MQ=255)
                                                                                                                                   ggCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACttt                    <  2:4695108/150‑1 (MQ=255)
                                                                                                                                       cGCAAAACGGGGGCTGTGCCCCCTTTGTGCCACCTTGACATCCTGCTGTCAGCATCGGCCCCGCAGCAGGGGCCGGCTGAAATGACAACAGCGGATTTGCAGACTTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTGTCCaa                <  2:4729601/150‑1 (MQ=255)
                                                                                                                                        gCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAag               >  1:2000740/1‑150 (MQ=255)
                                                                                                                                         cAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga              <  1:4805326/150‑1 (MQ=255)
                                                                                                                                         cAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga              <  2:1958131/150‑1 (MQ=255)
                                                                                                                                          aataaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCTGTATTGAACGCCCGTGATTGAGGGGTGGTATGGTTGAAAATCTTTCCAagag             >  2:221546/4‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAACAACTTTCCAagag             >  2:3725200/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag             >  1:4935637/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag             >  1:968770/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag             >  1:3004507/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag             >  2:3294308/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag             >  2:1307059/1‑150 (MQ=255)
                                                                                                                                            aaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGa           <  1:4304058/150‑1 (MQ=255)
                                                                                                                                                ggggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAtttt       <  2:2450130/150‑1 (MQ=255)
                                                                                                                                                ggggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAtttt       <  2:3622767/150‑1 (MQ=255)
                                                                                                                                                  ggCTTTGCCCCCTTTTTGCCACCTTGAAATTCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGa     >  1:4747582/1‑150 (MQ=255)
                                                                                                                                                    cTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGaaa   <  2:4853572/150‑1 (MQ=255)
                                                                                                                                                    cTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGaaa   <  2:2045796/150‑1 (MQ=255)
                                                                                                                                                     tttGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGACTGGGTGAAAAACTTTCCAAGAGGAGTTTGAAAt  >  2:4108254/1‑150 (MQ=255)
                                                                                                                                                     tttGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAt  >  2:1424301/1‑150 (MQ=255)
                                                                                                                                                     tttGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAt  <  1:400814/150‑1 (MQ=255)
                                                                                                                                                     tttGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAt  >  2:4693031/1‑150 (MQ=255)
                                                                                                                                                     tttGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAt  <  2:1860442/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TATCGAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785915‑1786213

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: