Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 51,908 C→G 5.7% E175Q (GAA→CAA) ‡ rsmA ← 16S rRNA m(6)2A1518,m(6)2A1519 dimethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091351,9080CG 5.7% 527.2 / 14.7 156E175Q (GAA→CAA) ‡rsmA16S rRNA m(6)2A1518,m(6)2A1519 dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (69/77);  new base G (7/2);  total (77/79)
Fisher's exact test for biased strand distribution p-value = 9.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 158 total aligned reads displayed.

GGTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGA  >  NC_000913/51759‑52056
                                                                                                                                                     |                                                                                                                                                    
ggTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTc                                                                                                                                                      >  1:1496115/1‑150 (MQ=255)
  ttttACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCa                                                                                                                                                    <  1:2485787/150‑1 (MQ=255)
  ttttACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCa                                                                                                                                                    >  2:2326258/1‑150 (MQ=255)
       cgacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTACCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTAcc                                                                                                                                               >  1:1620216/1‑150 (MQ=255)
       cgacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTAcc                                                                                                                                               >  1:1876480/1‑150 (MQ=255)
       cgacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGAACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGCCTGACGGCGGTACTTCAAGTAcc                                                                                                                                               >  2:882222/1‑150 (MQ=255)
        gacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCg                                                                                                                                              <  1:1474129/150‑1 (MQ=255)
         acgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgc                                                                                                                                             <  1:1411098/150‑2 (MQ=255)
          cgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCggg                                                                                                                                            >  2:723003/1‑150 (MQ=255)
           gCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGa                                                                                                                                           <  2:587163/150‑1 (MQ=255)
            cTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGAt                                                                                                                                          <  2:2223443/150‑1 (MQ=255)
            atGGTTAAAGGCTTCGGTGGTGATGCGGCGCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCTTGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTATAGGCTGACGGCGGTACTTCCAGTACCGGGAt                                                                                                                                          <  2:856850/149‑1 (MQ=255)
              ggTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATca                                                                                                                                        <  2:1630566/150‑1 (MQ=255)
               gTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                       >  2:1315434/1‑150 (MQ=255)
               gTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                       >  1:2440505/1‑150 (MQ=255)
                 tAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAt                                                                                                                                     >  2:836982/1‑150 (MQ=255)
                  aaaGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAtt                                                                                                                                    >  2:2041627/1‑150 (MQ=255)
                   aaGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTg                                                                                                                                   <  1:1278366/150‑1 (MQ=255)
                   aaGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTg                                                                                                                                   <  1:533076/150‑1 (MQ=255)
                     ggCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCa                                                                                                                                 >  1:746744/1‑150 (MQ=255)
                     ggCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCa                                                                                                                                 >  2:452559/1‑150 (MQ=255)
                       cTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAAt                                                                                                                               >  1:1151600/1‑150 (MQ=255)
                            gtggtgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTAt                                                                                                                          >  1:661473/1‑150 (MQ=255)
                                tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTg                                                                                                                        <  1:1741095/148‑1 (MQ=255)
                                 gATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTgcgc                                                                                                                     <  1:2081255/150‑1 (MQ=255)
                                 gATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTgcgc                                                                                                                     >  1:2327201/1‑150 (MQ=255)
                                        cTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACg                                                                                                              <  2:1946068/150‑1 (MQ=255)
                                           aaCACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCtt                                                                                                           <  1:1432196/150‑1 (MQ=255)
                                           aaCACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGCGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCtt                                                                                                           >  2:1226507/1‑150 (MQ=255)
                                             cacacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGTTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTaa                                                                                                         >  1:2368984/1‑150 (MQ=255)
                                             cacacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTaa                                                                                                         >  1:1157401/1‑150 (MQ=255)
                                               cacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATc                                                                                                       <  2:148772/150‑1 (MQ=255)
                                                acGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTTACTTAGGTACTTCCAGTACCGGTATCATATTGCAATATTATTGCGTGATGACGATTAAtag                                                                                                      <  2:112803/150‑3 (MQ=255)
                                                 cGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGa                                                                                                     <  1:1798574/150‑1 (MQ=255)
                                                  gAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAc                                                                                                    <  1:1632637/150‑1 (MQ=255)
                                                  gAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAc                                                                                                    <  2:1424317/150‑1 (MQ=255)
                                                   aaCATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                   <  2:1496115/150‑1 (MQ=255)
                                                    aCATCTTTAACCGGGTGAGGCATCGTTTCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTTGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCa                                                                                                  <  1:692703/150‑1 (MQ=255)
                                                      aTCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATa                                                                                                >  2:884729/1‑150 (MQ=255)
                                                        cTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACg                                                                                              >  1:1209400/1‑150 (MQ=255)
                                                           tAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCtt                                                                                           >  2:1457695/1‑150 (MQ=255)
                                                            aaCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCttt                                                                                          >  2:2460380/1‑150 (MQ=255)
                                                             aCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                         >  1:978949/1‑150 (MQ=255)
                                                             aCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                         >  2:2363334/1‑150 (MQ=255)
                                                                gggTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTg                                                                                      >  1:2259026/1‑150 (MQ=255)
                                                                   tGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTc                                                                                   <  2:14275/150‑1 (MQ=255)
                                                                    gAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCg                                                                                  >  2:1904004/1‑150 (MQ=255)
                                                                        cATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                              <  1:535133/150‑1 (MQ=255)
                                                                        cATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                              <  1:762494/150‑1 (MQ=255)
                                                                            gTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCa                                                                          >  2:1765616/1‑150 (MQ=255)
                                                                              tGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAc                                                                        >  2:2039845/1‑150 (MQ=255)
                                                                              tGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAc                                                                        >  1:214102/1‑150 (MQ=255)
                                                                               gCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAcc                                                                       <  1:271313/150‑1 (MQ=255)
                                                                               gCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTACAGTACCGGGATCACATTGCAATAGTATTGCGACATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAca                                                                       >  2:138037/1‑149 (MQ=255)
                                                                                  tGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGa                                                                    <  2:696123/150‑1 (MQ=255)
                                                                                   gAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTTGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGAc                                                                   <  1:559766/150‑1 (MQ=255)
                                                                                   gAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGAc                                                                   <  1:2480612/150‑1 (MQ=255)
                                                                                    aGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACg                                                                  >  1:2252040/1‑150 (MQ=255)
                                                                                        aCCAGGCGCACGCAGGAGGACTCCAGTTTGGGAGGTGGTGTAAAGGAAGACGGCGGTCGGGCCAGTCCCGGGATCACATTGCCATAGTATTGCGCCATGACGCTTATTCGCCCATACGGTTTGCTGTTCGGTCCTGCAACCAGACTAtta                                                              >  2:2218314/1‑149 (MQ=255)
                                                                                         ccAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTca                                                             <  1:52459/150‑1 (MQ=255)
                                                                                         ccAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTca                                                             <  1:415433/150‑1 (MQ=255)
                                                                                          cAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcac                                                            <  1:881846/150‑1 (MQ=255)
                                                                                           aGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacc                                                           <  1:1467287/150‑1 (MQ=255)
                                                                                                  cgacgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                                                    >  2:2399439/1‑150 (MQ=255)
                                                                                                  cgacgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                                                    >  1:2039710/1‑150 (MQ=255)
                                                                                                    acgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtttt                                                  <  2:1130679/150‑1 (MQ=255)
                                                                                                       gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                                               >  1:82875/1‑150 (MQ=255)
                                                                                                        gttaaTCCACTTTGTGTGGTGGTGTAAAGGCTGAAGGCGGTACTTCCAGTAACGGGATCACATTGCAATAGTATTGCTCCATGACGCTTAATCGACCATACGCTTTGCTTTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCa                                              <  2:1620216/147‑1 (MQ=255)
                                                                                                        cggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGGTTCACCACCTCTTTTTGCa                                              <  2:303367/150‑1 (MQ=255)
                                                                                                        cggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCa                                              <  1:2326258/150‑1 (MQ=255)
                                                                                                        cggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCa                                              <  2:612686/150‑1 (MQ=255)
                                                                                                          gAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAAc                                            >  2:1742472/1‑150 (MQ=255)
                                                                                                                  tttGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAAAAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg                                    <  2:2456184/150‑1 (MQ=255)
                                                                                                                    tGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCa                                  >  1:1082711/1‑150 (MQ=255)
                                                                                                                            ggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCa                          >  1:2463214/1‑150 (MQ=255)
                                                                                                                              tgtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAAt                        <  1:1042066/150‑1 (MQ=255)
                                                                                                                              tgtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAAt                        >  1:520132/1‑150 (MQ=255)
                                                                                                                               gtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTGTTTGCAACATAAAGTGCATGTCGGCAATg                       >  2:1019473/1‑150 (MQ=255)
                                                                                                                                 aaaGGGTGACTGCGGTACTTCCAGTACCGGTATAACATTGCAATAGTATTGCGCAATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGc                     <  2:253849/150‑1 (MQ=255)
                                                                                                                                    ggCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATc                  >  1:114218/1‑150 (MQ=255)
                                                                                                                                    ggCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATc                  <  1:1969240/150‑1 (MQ=255)
                                                                                                                                    ggCTGACCGCCGTACTTCCAATACCGGGATCACATTGCAATAGTATTGCGCCAAGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATc                  >  2:2311691/1‑150 (MQ=255)
                                                                                                                                     gCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGACATGACGCTTAATCGACCATACGCTTTGATGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCa                 <  2:2440505/150‑1 (MQ=255)
                                                                                                                                      cttactgcggTACTTCCAGTAGGGGGATCAGATTGCAATAGTTTTGCGCCATTACGCTTAATCGTCGATACGGTTTGGTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAg                <  2:2152445/144‑1 (MQ=255)
                                                                                                                                      cTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAg                <  2:1736424/150‑1 (MQ=255)
                                                                                                                                      cTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAg                >  1:2350289/1‑150 (MQ=255)
                                                                                                                                       tGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGt               <  1:1315434/150‑1 (MQ=255)
                                                                                                                                       tGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGt               <  2:776227/150‑1 (MQ=255)
                                                                                                                                          tcccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata            >  2:972912/4‑150 (MQ=255)
                                                                                                                                          tcccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata            >  1:1649873/4‑150 (MQ=255)
                                                                                                                                          tcccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata            >  2:1190899/4‑150 (MQ=255)
                                                                                                                                           cccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg           >  1:2077864/3‑150 (MQ=255)
                                                                                                                                               gTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTaa       >  2:661159/1‑150 (MQ=255)
                                                                                                                                                 aCTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCa                                  >  1:1138505/1‑121 (MQ=255)
                                                                                                                                                 aCTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCa                                  <  2:1138505/121‑1 (MQ=255)
                                                                                                                                                 aCTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc     <  1:390760/150‑1 (MQ=255)
                                                                                                                                                  cTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACa    <  2:381955/150‑1 (MQ=255)
                                                                                                                                                  cTTCAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACTCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACa    <  1:882222/150‑1 (MQ=255)
                                                                                                                                                    tCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGa  >  1:1921922/1‑150 (MQ=255)
                                                                                                                                                    tCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGa  <  1:2399439/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGA  >  NC_000913/51759‑52056

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