Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 51,909 C→A 9.7% L174L (CTG→CTT)  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,9090CA 9.7% 339.2 / 20.8 120L174L (CTG→CTT) rsmA16S rRNA m(6)A1518, m(6)A1519 dimethyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (60/48);  new base A (6/6);  total (66/54)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 121 total aligned reads displayed.

TTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGATG  >  NC_000913/51761‑52058
                                                                                                                                                    |                                                                                                                                                     
ttttACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCa                                                                                                                                                      >  1:1512199/1‑150 (MQ=255)
 tttACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCat                                                                                                                                                     >  1:971330/1‑149 (MQ=255)
 tttACGACGCTGGTTAAAGGATTCGGTGGGGATGCGGCGCAACACACGAACAGCTTTAACCGGGTGAGGCATCGTTGCATGCGGAACCAGGCGCTCGACGGCGGAATCCACTTTGGGTGGTTGTGTAATGGCTGACGGCGGTACTTCCAg                                                                                                                                                     <  1:130391/150‑1 (MQ=255)
  ttACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGt                                                                                                                                                    >  2:1113534/1‑150 (MQ=255)
    acgacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTAc                                                                                                                                                  <  1:999202/150‑1 (MQ=255)
      gacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCg                                                                                                                                                <  2:143574/150‑1 (MQ=255)
      gacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCg                                                                                                                                                <  1:1628277/150‑1 (MQ=255)
      gacgCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCg                                                                                                                                                <  1:1624888/150‑1 (MQ=255)
         gCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgccg                                                                                                                                             <  1:494047/150‑4 (MQ=255)
         gCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgccg                                                                                                                                             <  1:825382/150‑4 (MQ=255)
          cTGGTTAACGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGAt                                                                                                                                            <  1:1570289/150‑1 (MQ=255)
          cTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGAt                                                                                                                                            >  1:1168421/1‑150 (MQ=255)
            ggTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATca                                                                                                                                          <  1:633225/150‑1 (MQ=255)
            ggTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgccgtca                                                                                                                                          <  2:263961/150‑7 (MQ=255)
             gTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACTTCTTTAACCGGGTGAGGCATCGTTGAATTAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgccgtcag                                                                                                                                         <  1:111830/150‑8 (MQ=255)
             gTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACGTTGGGTGGTGGTGTAAGGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                         >  1:343213/1‑150 (MQ=255)
             gTTAAAGGCTTCGGTGGGGATGCGGCTCAACACACGAACATCTTTACCCGGGTGAGGCATCGTTGCATGAGGTACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATcac                                                                                                                                         <  2:631340/150‑1 (MQ=255)
               tAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACAt                                                                                                                                       >  2:1539691/1‑150 (MQ=255)
                  aGGCTACGGTGGTGAGGCGCCTCAACACACGAACATCTTTAAGCGGGTGAGGCATGGTTGCAAGAGGAACCAGCCGCAGGCCGGCGGCATCCCCTTTGGGTGGTGGTGTAAAGGCTGCCGGCGGAACGACCAGAACCGGGAGAACATTga                                                                                                                                    >  2:177353/1‑149 (MQ=255)
                   ggCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGTAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGCCTGACGCCGGTACTTCAAGTACCGGGATCACATTGCa                                                                                                                                   >  2:1643101/1‑150 (MQ=255)
                     cTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAAt                                                                                                                                 >  1:1286827/1‑150 (MQ=255)
                       tCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAg                                                                                                                               <  1:1539691/150‑1 (MQ=255)
                           tggtgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTAtt                                                                                                                           >  2:1641879/1‑150 (MQ=255)
                           tggtgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTAtt                                                                                                                           >  2:1635294/1‑150 (MQ=255)
                              tgATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTGGAAGTACCgccg                                                                                                                                             >  2:61530/1‑126 (MQ=255)
                                 tGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCa                                                                                                                     <  2:377831/150‑1 (MQ=255)
                                  gCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGGAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCAt                                                                                                                    >  2:1210455/1‑150 (MQ=255)
                                   cGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGACCCAGGCTCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATg                                                                                                                   >  2:538188/1‑150 (MQ=255)
                                    ggCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGa                                                                                                                  <  2:1514034/150‑1 (MQ=255)
                                             cacGAACATCTTTAGCCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATc                                                                                                         >  1:325463/1‑150 (MQ=255)
                                             cacGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATc                                                                                                         >  1:464018/1‑150 (MQ=255)
                                               cGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGa                                                                                                       >  1:579777/1‑150 (MQ=255)
                                               cGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGa                                                                                                       >  2:766183/1‑150 (MQ=255)
                                                gAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAc                                                                                                      <  2:877034/150‑1 (MQ=255)
                                                   cATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCAt                                                                                                   <  2:1595438/150‑1 (MQ=255)
                                                    aTCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATa                                                                                                  >  2:461970/1‑150 (MQ=255)
                                                     tCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATAc                                                                                                 <  1:932416/150‑1 (MQ=255)
                                                        ttAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCt                                                                                              >  2:828118/1‑150 (MQ=255)
                                                            ccGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                          <  1:1489506/150‑1 (MQ=255)
                                                            ccGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGc                                                                                          <  1:184829/150‑1 (MQ=255)
                                                             cGGGTGAGGCATCTTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGTTGGTGTAAAGGCTGACGGCGGTACTTCCATTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCt                                                                                         <  2:374563/150‑1 (MQ=255)
                                                              gggTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTg                                                                                        >  2:596990/1‑150 (MQ=255)
                                                                  gAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCg                                                                                    <  2:403970/150‑1 (MQ=255)
                                                                  gAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCg                                                                                    <  2:340021/150‑1 (MQ=255)
                                                                      cATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                                >  2:714272/1‑150 (MQ=255)
                                                                      cATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGGGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTcc                                                                                <  1:196740/150‑1 (MQ=255)
                                                                        tCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTTTTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTg                                                                              <  1:1261245/150‑1 (MQ=255)
                                                                          gTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCa                                                                            <  1:1663991/150‑1 (MQ=255)
                                                                          gTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCa                                                                            >  1:1558483/1‑150 (MQ=255)
                                                                                 gAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGAc                                                                     <  2:1607425/150‑1 (MQ=255)
                                                                                  aGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACg                                                                    >  2:1240888/1‑150 (MQ=255)
                                                                                  aGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATAGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACg                                                                    >  1:173479/1‑150 (MQ=255)
                                                                                     aaCCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGAtt                                                                 >  2:1207317/1‑150 (MQ=255)
                                                                                      aCCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTc                                                                >  1:367993/1‑150 (MQ=255)
                                                                                      aCCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTc                                                                >  2:1435597/1‑150 (MQ=255)
                                                                                       ccAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTca                                                               <  1:1207317/150‑1 (MQ=255)
                                                                                        cAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                           >  2:296650/1‑121 (MQ=255)
                                                                                        cAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTg                                                                                           <  1:296650/121‑1 (MQ=255)
                                                                                         aGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacc                                                             >  1:470084/1‑150 (MQ=255)
                                                                                         aGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacc                                                             <  2:1512199/150‑1 (MQ=255)
                                                                                          ggCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacca                                                            <  2:1168421/150‑1 (MQ=255)
                                                                                           gcgcACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac                                                           >  1:375222/1‑150 (MQ=255)
                                                                                                cgacgGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                                                      >  1:118294/1‑150 (MQ=255)
                                                                                                    ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac                                                           >  1:1341630/1‑141 (MQ=255)
                                                                                                    ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac                                                           <  2:1341630/141‑1 (MQ=255)
                                                                                                    ggcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATCCACCTCCTCTTTTTg                                                  >  1:161227/1‑150 (MQ=255)
                                                                                                     gcggAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCACCCCGCCGAGTCCCCCCCTCTTTTTGc                                                 >  1:1215431/1‑150 (MQ=255)
                                                                                                        gAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAAc                                              <  1:1014145/150‑1 (MQ=255)
                                                                                                        gAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAAc                                              >  1:119498/1‑150 (MQ=255)
                                                                                                         aaTCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACATCTTTTTGCAACa                                             >  1:124255/1‑150 (MQ=255)
                                                                                                            ccACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCCGACGATTCACCACCTCTTTTTGCAACATaa                                          >  2:536432/1‑150 (MQ=255)
                                                                                                            ccACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATaa                                          >  2:1304076/1‑150 (MQ=255)
                                                                                                              aCTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAg                                        >  2:799543/1‑150 (MQ=255)
                                                                                                                tttGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCCCCTCTATTTGCAACATAAAGTg                                      >  1:444301/1‑150 (MQ=255)
                                                                                                                  tGGGTGGTGGTGTAAAGGCTGACGGCGTTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACTCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCa                                    <  1:225236/150‑1 (MQ=255)
                                                                                                                  tGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCa                                    <  1:548329/150‑1 (MQ=255)
                                                                                                                    ggtggtggtGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATg                                  >  2:1443846/1‑150 (MQ=255)
                                                                                                                    ggtggtggtGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATg                                  <  1:1113534/150‑1 (MQ=255)
                                                                                                                    ggtggtggtGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATg                                  <  2:470084/150‑1 (MQ=255)
                                                                                                                      gggtggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTc                                <  2:1428645/149‑1 (MQ=255)
                                                                                                                         tggttTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc                             <  2:1419823/150‑1 (MQ=255)
                                                                                                                          ggtgTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCa                            >  2:1298729/1‑150 (MQ=255)
                                                                                                                            tgtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAAt                          >  2:1519497/1‑150 (MQ=255)
                                                                                                                             gtAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATg                         >  1:493028/1‑150 (MQ=255)
                                                                                                                              tAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATgg                        <  2:367993/150‑1 (MQ=255)
                                                                                                                                    gtGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACTCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAg                  <  1:1643101/149‑1 (MQ=255)
                                                                                                                                    cTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAg                  >  1:435731/1‑150 (MQ=255)
                                                                                                                                     tGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGt                 >  1:1008311/1‑150 (MQ=255)
                                                                                                                                      gACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCGTACTCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGta                >  1:877445/1‑150 (MQ=255)
                                                                                                                                       aCGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGATTTGCGCCATGACGCTTCATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGtat               <  2:278992/150‑1 (MQ=255)
                                                                                                                                         cccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg             >  1:1098209/3‑150 (MQ=255)
                                                                                                                                         cccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTGTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAg             >  1:684112/3‑150 (MQ=255)
                                                                                                                                          gcggTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGc            <  2:1351401/150‑1 (MQ=255)
                                                                                                                                          ccggTACTGGAAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGc            >  1:724315/2‑150 (MQ=255)
                                                                                                                                            ggTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTa          >  1:944530/1‑150 (MQ=255)
                                                                                                                                            ggTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTa          <  2:1229981/150‑1 (MQ=255)
                                                                                                                                              tACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGAACATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGACGGCAATGGCATCAGTATAGCTaaa        <  2:882539/150‑1 (MQ=255)
                                                                                                                                               aCTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAAc       >  2:1522360/1‑150 (MQ=255)
                                                                                                                                                cTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACa      <  1:461970/150‑1 (MQ=255)
                                                                                                                                                    cAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACGATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGATg  <  2:325463/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTTACGACGCTGGTTAAAGGCTTCGGTGGTGATGCGGCTCAACACACGAACATCTTTAACCGGGTGAGGCATCGTTGCATGAGGAACCAGGCGCACGACGGCGGAATCCACTTTGGGTGGTGGTGTAAAGGCTGACGGCGGTACTTCCAGTACCGGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGCAATGGCATCAGTATAGCTAAACAGATG  >  NC_000913/51761‑52058

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