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.

GTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGA  >  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)
                       ttCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCTCTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATa                                                                                                                              <  1:676118/150‑1 (MQ=255)
                          ggtggtGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTa                                                                                                                           >  1:659050/1‑150 (MQ=255)
                           gtggtgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGCCGGCGGTACTTCCAGTCCCGGGATCCCATTGC‑AATAGTAt                                                                                                                          >  2:107085/1‑150 (MQ=255)
                               tgATGCGGCTCAACACACGAACATCTTTACCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGCCGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTgcg                                                                                                                      >  1:1088843/1‑150 (MQ=255)
                               tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAATGC‑AATAGTATTgcg                                                                                                                      <  2:1359496/150‑1 (MQ=255)
                               tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTGGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTgcg                                                                                                                      <  2:1146442/150‑1 (MQ=255)
                               ggATGCGGCTCAACACACGAACATCTTTACCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTgcg                                                                                                                      <  2:1657661/149‑1 (MQ=255)
                                   gCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCAGCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGAACCGGGATCACATTGC‑AATAGTATTGCGCCAt                                                                                                                  <  1:881651/150‑1 (MQ=255)
                                     ggCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGCCGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGa                                                                                                                >  2:1631282/1‑150 (MQ=255)
                                     ggCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGa                                                                                                                <  1:164927/150‑1 (MQ=255)
                                      gCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACATCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGAc                                                                                                               <  1:1138229/150‑1 (MQ=255)
                                       cTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGAGCCATGACg                                                                                                              >  2:218610/1‑150 (MQ=255)
                                           acacacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCACACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTa                                                                                                          >  1:1918138/1‑150 (MQ=255)
                                             acacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTTC‑AATAGTATTGCGCCATGACGCTTaaa                                                                                                        >  1:175361/1‑149 (MQ=255)
                                              cacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAAATAGTATTGCGCCATGACGCTTAAt                                                                                                        >  2:972971/1‑150 (MQ=255)
                                              cacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGGGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATc                                                                                                       <  2:1335778/150‑1 (MQ=255)
                                               acGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGTACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCg                                                                                                      <  2:783841/150‑1 (MQ=255)
                                                 gAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGAc                                                                                                    >  1:600569/1‑150 (MQ=255)
                                                  aaCATCTTTATCCGGGTGAGGCATCGTTGCATGAGTAACCAGGCGCACTACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACTGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                   <  2:1285480/150‑1 (MQ=255)
                                                  aaCATCTTTAACCGGGTGAGGCATCGTTTCATTAGGAACCAGGCGCACGACGGCGGAATCCTCTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                   <  2:1458655/150‑1 (MQ=255)
                                                  aaCATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGAcc                                                                                                   <  2:462422/150‑1 (MQ=255)
                                                   aCATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCa                                                                                                  >  1:1679397/1‑150 (MQ=255)
                                                      tCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATAc                                                                                               >  2:918334/1‑150 (MQ=255)
                                                       cTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACg                                                                                              >  1:1647176/1‑150 (MQ=255)
                                                        tttAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTCCCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGc                                                                                             >  2:1936380/1‑150 (MQ=255)
                                                          tAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACATTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGAAGCTTAATCGACCATACGCtt                                                                                           >  1:556638/1‑150 (MQ=255)
                                                           aaCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGCCCATACGCttt                                                                                          >  2:1629384/1‑150 (MQ=255)
                                                            aCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACTGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                         >  2:294547/1‑150 (MQ=255)
                                                            aCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                         >  2:1075323/1‑150 (MQ=255)
                                                             ccGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                        <  1:1295531/150‑1 (MQ=255)
                                                             ccGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                        <  1:2045300/150‑1 (MQ=255)
                                                                    aGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCgg                                                                                 >  2:562872/1‑150 (MQ=255)
                                                                    aGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCgg                                                                                 >  1:1365431/1‑150 (MQ=255)
                                                                      gCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTc                                                                               >  2:377035/1‑150 (MQ=255)
                                                                       cATCGTTGCATGAGGATCCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTAGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                              <  2:457865/150‑1 (MQ=255)
                                                                       cATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                              <  2:397389/150‑1 (MQ=255)
                                                                          cGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGc                                                                           <  2:1255610/150‑1 (MQ=255)
                                                                          cGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGc                                                                           >  1:2003460/1‑150 (MQ=255)
                                                                             tGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAAc                                                                        >  1:487076/1‑150 (MQ=255)
                                                                               cATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCa                                                                      >  1:1432490/1‑150 (MQ=255)
                                                                               cATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCa                                                                      >  1:1860274/1‑150 (MQ=255)
                                                                                  gAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGAc                                                                   <  2:724536/150‑1 (MQ=255)
                                                                                          aGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacc                                                           <  2:613275/150‑1 (MQ=255)
                                                                                             cgcACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCaccacc                                                        <  2:600569/150‑1 (MQ=255)
                                                                                              gcACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACct                                                       >  1:850631/1‑150 (MQ=255)
                                                                                               cACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCCTACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACctc                                                      <  2:961118/150‑1 (MQ=255)
                                                                                               cACGACGGCGGAATCCACTTTGGGGGGTGGTGTAAAGGCGGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACctc                                                      <  2:1676729/150‑1 (MQ=255)
                                                                                                 cgacgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                                                    >  1:950779/1‑150 (MQ=255)
                                                                                                 cgacgGCGGAATACACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                                                    <  2:1982016/150‑1 (MQ=255)
                                                                                                   acgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtttt                                                  >  1:1715826/1‑150 (MQ=255)
                                                                                                     ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg                                                <  2:667468/150‑1 (MQ=255)
                                                                                                     ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg                                                <  2:1712042/150‑1 (MQ=255)
                                                                                                     ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGAACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg                                                >  2:872730/1‑150 (MQ=255)
                                                                                                      gcggATTCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                                               <  1:982268/150‑1 (MQ=255)
                                                                                                      gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTTATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTgg                                               <  2:869108/150‑2 (MQ=255)
                                                                                                      gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                                               <  1:519416/150‑1 (MQ=255)
                                                                                                      gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGGCCTGCAACCAGACGATTCACCACCTCTTTTTGc                                               <  1:1766198/150‑1 (MQ=255)
                                                                                                       cggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTACa                                              <  2:123386/150‑1 (MQ=255)
                                                                                                        ggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCTCATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCCACCAGACGATTCACCACCTCTTTTTGcta                                             >  2:419746/1‑148 (MQ=255)
                                                                                                         gAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAAc                                            <  1:2158872/150‑1 (MQ=255)
                                                                                                           aTCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCAGGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCTACAt                                          >  1:998986/1‑150 (MQ=255)
                                                                                                             ccACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATaa                                        >  2:1856423/1‑150 (MQ=255)
                                                                                                              cACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATaaa                                       >  1:1412382/1‑150 (MQ=255)
                                                                                                                cTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGt                                     >  1:2118298/1‑150 (MQ=255)
                                                                                                                   tGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCa                                  <  1:217049/150‑1 (MQ=255)
                                                                                                                     ggtggttgtGTAAAGGATGACGGCGGTACTTCCAGTACCTGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGAACATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATg                                <  2:110664/150‑1 (MQ=255)
                                                                                                                       tggtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTc                              >  2:1195016/1‑150 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  1:1698102/150‑1 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  1:1911777/150‑1 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGA‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  2:2209920/150‑1 (MQ=255)
                                                                                                                         gtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGAGCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCgg                            <  1:343070/150‑1 (MQ=255)
                                                                                                                          tggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc                           >  2:233435/1‑150 (MQ=255)
                                                                                                                            gtgtAAAGGCTGCCGGCTGTACTGCCAGTAACGGGATAACATTGC‑TATAGTATTGCGCCATGAAGCTTAATCGACCATACGCTTTGCTGTTCGGTACTGCAACCAGACGATTCACCACCTCTTTTTGCAACAAAAAGTGCATGTCGGCaa                         >  1:1073778/1‑150 (MQ=255)
                                                                                                                             tgtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCATGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAAt                        >  1:1325112/1‑150 (MQ=255)
                                                                                                                               tAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATgg                      >  2:1442946/1‑150 (MQ=255)
                                                                                                                                     cTGACGGCGGTACTGCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAg                <  2:1543352/150‑1 (MQ=255)
                                                                                                                                      tGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGt               <  2:1919906/150‑1 (MQ=255)
                                                                                                                                         tcccggTACTGGAAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata            >  1:11455/4‑150 (MQ=255)
                                                                                                                                         cggcggTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtata            >  1:564328/1‑150 (MQ=255)
                                                                                                                                          ggcggTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg           <  1:236516/150‑1 (MQ=255)
                                                                                                                                          cccggTACTGGAAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg           >  2:559795/3‑150 (MQ=255)
                                                                                                                                            cggTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCt         <  1:1440651/150‑1 (MQ=255)
                                                                                                                                               tACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTaaa      >  1:1210734/1‑150 (MQ=255)
                                                                                                                                                aCTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc     <  2:2169380/150‑1 (MQ=255)
                                                                                                                                                aCTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc     <  2:1959988/150‑1 (MQ=255)
                                                                                                                                                aCTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc     <  2:2119415/150‑1 (MQ=255)
                                                                                                                                                   tCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGa  >  1:1188874/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GTTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGC‑AATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGA  >  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: ‑