Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 51,907 T→G 6.0% E175A (GAA→GCA)  rsmA ← 16S rRNA m(6)A1518, m(6)A1519 dimethyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091351,9070TG 6.0% 466.0 / 11.1 154E175A (GAA→GCA) rsmA16S rRNA m(6)A1518, m(6)A1519 dimethyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (70/74);  new base G (6/4);  total (76/78)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

GTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGA  >  NC_000913/51760‑52056
                                                                                                                                                   |                                                                                                                                                      
gTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACTTCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTcc                                                                                                                                                      <  2:1840993/150‑1 (MQ=255)
gTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTcc                                                                                                                                                      <  1:911417/150‑1 (MQ=255)
 ttttACGACGCGGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGTGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCa                                                                                                                                                     <  2:1565793/150‑1 (MQ=255)
   ttACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCAGGGCGGCGGAATCCACTTTGGGTGGTGGTGTATAGGCTGCCGGCGGTCCGGCCAGt                                                                                                                                                   >  1:395090/1‑150 (MQ=255)
         cgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgcc                                                                                                                                             <  1:2135229/150‑3 (MQ=255)
          gCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGTCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGa                                                                                                                                            <  2:1393815/150‑1 (MQ=255)
          gCTGGCTAACGGCTTCGGTGGTGATGCGGCTCAACACACGAACAGCATTAACCGGGGGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTATGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgccg                                                                                                                                            <  2:256532/150‑4 (MQ=255)
             ggTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATca                                                                                                                                         >  1:123386/1‑150 (MQ=255)
             ggTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATca                                                                                                                                         <  2:1320116/150‑1 (MQ=255)
             ggTTAAAGGCTTCGGTGGTGAGGCGGCTCAACACACGAACATCTTTAACCGGGGGAGGCATCGGTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGGGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATca                                                                                                                                         <  1:1965965/150‑1 (MQ=255)
              gTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                        >  2:1911777/1‑150 (MQ=255)
                tAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAt                                                                                                                                      >  1:1826657/1‑150 (MQ=255)
                   aGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGGAGCATGAGGAACCAGGCGCACGCCGGCGGAAGCCACGTGGGGGGGGGGTGTTACGGCTGCCGGCGGTACTGCCAGAACAGAGATCTCAGTGc                                                                                                                                   >  2:381767/1‑150 (MQ=255)
                       ttCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCTCTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATa                                                                                                                              <  1:676118/150‑1 (MQ=255)
                          ggtggtGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTa                                                                                                                           >  1:659050/1‑150 (MQ=255)
                           gtggtgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGCCGGCGGTACTTCCAGTCCCGGGATCCCATTGCAATAGTAt                                                                                                                          >  2:107085/1‑150 (MQ=255)
                               tgATGCGGCTCAACACACGAACATCTTTACCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGCCGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTgcg                                                                                                                      >  1:1088843/1‑150 (MQ=255)
                               tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAATGCAATAGTATTgcg                                                                                                                      <  2:1359496/150‑1 (MQ=255)
                               tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTGGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTgcg                                                                                                                      <  2:1146442/150‑1 (MQ=255)
                               ggATGCGGCTCAACACACGAACATCTTTACCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTgcg                                                                                                                      <  2:1657661/149‑1 (MQ=255)
                                   gCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCAGCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGAACCGGGATCACATTGCAATAGTATTGCGCCAt                                                                                                                  <  1:881651/150‑1 (MQ=255)
                                     ggCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGCCGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGa                                                                                                                >  2:1631282/1‑150 (MQ=255)
                                     ggCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGa                                                                                                                <  1:164927/150‑1 (MQ=255)
                                      gCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACATCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGAc                                                                                                               <  1:1138229/150‑1 (MQ=255)
                                       cTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGAGCCATGACg                                                                                                              >  2:218610/1‑150 (MQ=255)
                                           acacacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCACACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTa                                                                                                          >  1:1918138/1‑150 (MQ=255)
                                             acacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTTCAATAGTATTGCGCCATGACGCTTaaa                                                                                                        >  1:175361/1‑149 (MQ=255)
                                              cacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAAATAGTATTGCGCCATGACGCTTAAt                                                                                                        >  2:972971/1‑150 (MQ=255)
                                              cacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGGGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATc                                                                                                       <  2:1335778/150‑1 (MQ=255)
                                               acGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGTACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCg                                                                                                      <  2:783841/150‑1 (MQ=255)
                                                 gAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAc                                                                                                    >  1:600569/1‑150 (MQ=255)
                                                  aaCATCTTTATCCGGGTGAGGCATCGTTGCATGAGTAACCAGGCGCACTACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACTGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                   <  2:1285480/150‑1 (MQ=255)
                                                  aaCATCTTTAACCGGGTGAGGCATCGTTTCATTAGGAACCAGGCGCACGACGGCGGAATCCTCTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                   <  2:1458655/150‑1 (MQ=255)
                                                  aaCATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                   <  2:462422/150‑1 (MQ=255)
                                                   aCATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCa                                                                                                  >  1:1679397/1‑150 (MQ=255)
                                                      tCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATAc                                                                                               >  2:918334/1‑150 (MQ=255)
                                                       cTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACg                                                                                              >  1:1647176/1‑150 (MQ=255)
                                                        tttAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTCCCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGc                                                                                             >  2:1936380/1‑150 (MQ=255)
                                                          tAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACATTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGAAGCTTAATCGACCATACGCtt                                                                                           >  1:556638/1‑150 (MQ=255)
                                                           aaCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGCCCATACGCttt                                                                                          >  2:1629384/1‑150 (MQ=255)
                                                            aCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACTGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                         >  2:294547/1‑150 (MQ=255)
                                                            aCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                         >  2:1075323/1‑150 (MQ=255)
                                                             ccGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                        <  1:1295531/150‑1 (MQ=255)
                                                             ccGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                        <  1:2045300/150‑1 (MQ=255)
                                                                    aGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCgg                                                                                 >  2:562872/1‑150 (MQ=255)
                                                                    aGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCgg                                                                                 >  1:1365431/1‑150 (MQ=255)
                                                                      gCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTc                                                                               >  2:377035/1‑150 (MQ=255)
                                                                       cATCGTTGCATGAGGATCCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTAGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                              <  2:457865/150‑1 (MQ=255)
                                                                       cATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                              <  2:397389/150‑1 (MQ=255)
                                                                          cGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGc                                                                           <  2:1255610/150‑1 (MQ=255)
                                                                          cGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGc                                                                           >  1:2003460/1‑150 (MQ=255)
                                                                             tGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAc                                                                        >  1:487076/1‑150 (MQ=255)
                                                                               cATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCa                                                                      >  1:1432490/1‑150 (MQ=255)
                                                                               cATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCa                                                                      >  1:1860274/1‑150 (MQ=255)
                                                                                  gAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGAc                                                                   <  2:724536/150‑1 (MQ=255)
                                                                                          aGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacc                                                           <  2:613275/150‑1 (MQ=255)
                                                                                             cgcACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCaccacc                                                        <  2:600569/150‑1 (MQ=255)
                                                                                              gcACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACct                                                       >  1:850631/1‑150 (MQ=255)
                                                                                               cACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCCTACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACctc                                                      <  2:961118/150‑1 (MQ=255)
                                                                                               cACGACGGCGGAATCCACTTTGGGGGGTGGTGTAAAGGCGGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACctc                                                      <  2:1676729/150‑1 (MQ=255)
                                                                                                 cgacgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                                                    >  1:950779/1‑150 (MQ=255)
                                                                                                 cgacgGCGGAATACACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                                                    <  2:1982016/150‑1 (MQ=255)
                                                                                                   acgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtttt                                                  >  1:1715826/1‑150 (MQ=255)
                                                                                                     ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg                                                <  2:667468/150‑1 (MQ=255)
                                                                                                     ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg                                                <  2:1712042/150‑1 (MQ=255)
                                                                                                     ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGAACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg                                                >  2:872730/1‑150 (MQ=255)
                                                                                                      gcggATTCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                                               <  1:982268/150‑1 (MQ=255)
                                                                                                      gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTTATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTgg                                               <  2:869108/150‑2 (MQ=255)
                                                                                                      gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                                               <  1:519416/150‑1 (MQ=255)
                                                                                                      gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGGCCTGCAACCAGACGATTCACCACCTCTTTTTGc                                               <  1:1766198/150‑1 (MQ=255)
                                                                                                       cggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTACa                                              <  2:123386/150‑1 (MQ=255)
                                                                                                        ggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCTCATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCCACCAGACGATTCACCACCTCTTTTTGcta                                             >  2:419746/1‑148 (MQ=255)
                                                                                                         gAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAAc                                            <  1:2158872/150‑1 (MQ=255)
                                                                                                           aTCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCAGGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCTACAt                                          >  1:998986/1‑150 (MQ=255)
                                                                                                             ccACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATaa                                        >  2:1856423/1‑150 (MQ=255)
                                                                                                              cACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATaaa                                       >  1:1412382/1‑150 (MQ=255)
                                                                                                                cTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGt                                     >  1:2118298/1‑150 (MQ=255)
                                                                                                                   tGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCa                                  <  1:217049/150‑1 (MQ=255)
                                                                                                                     ggtggttgtGTAAAGGATGACGGCGGTACTTCCAGTACCTGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAACATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATg                                <  2:110664/150‑1 (MQ=255)
                                                                                                                       tggtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTc                              >  2:1195016/1‑150 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  1:1698102/150‑1 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  1:1911777/150‑1 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGAAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  2:2209920/150‑1 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGAGCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  1:343070/150‑1 (MQ=255)
                                                                                                                          tggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc                           >  2:233435/1‑150 (MQ=255)
                                                                                                                            gtgtAAAGGCTGCCGGCTGTACTGCCAGTAACGGGATAACATTGCTATAGTATTGCGCCATGAAGCTTAATCGACCATACGCTTTGCTGTTCGGTACTGCAACCAGACGATTCACCACCTCTTTTTGCAACAAAAAGTGCATGTCGGCaa                         >  1:1073778/1‑150 (MQ=255)
                                                                                                                             tgtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCATGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAAt                        >  1:1325112/1‑150 (MQ=255)
                                                                                                                               tAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATgg                      >  2:1442946/1‑150 (MQ=255)
                                                                                                                                     cTGACGGCGGTACTGCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAg                <  2:1543352/150‑1 (MQ=255)
                                                                                                                                      tGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGt               <  2:1919906/150‑1 (MQ=255)
                                                                                                                                         tcccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata            >  1:11455/4‑150 (MQ=255)
                                                                                                                                         cggcggTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata            >  1:564328/1‑150 (MQ=255)
                                                                                                                                          ggcggTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg           <  1:236516/150‑1 (MQ=255)
                                                                                                                                          cccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg           >  2:559795/3‑150 (MQ=255)
                                                                                                                                            cggTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCt         <  1:1440651/150‑1 (MQ=255)
                                                                                                                                               tACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTaaa      >  1:1210734/1‑150 (MQ=255)
                                                                                                                                                aCTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc     <  2:2169380/150‑1 (MQ=255)
                                                                                                                                                aCTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc     <  2:1959988/150‑1 (MQ=255)
                                                                                                                                                aCTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc     <  2:2119415/150‑1 (MQ=255)
                                                                                                                                                   tCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGa  >  1:1188874/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGA  >  NC_000913/51760‑52056

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