Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,034 G→A 6.0% intergenic (+28/‑99) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0340GA 6.0% 517.3 / 13.3 166intergenic (+28/‑99)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base G (81/73);  new base A (6/4);  total (87/79)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

Only 100 of 169 total aligned reads displayed.

GCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAA  >  NC_000913/4409885‑4410183
                                                                                                                                                     |                                                                                                                                                       
gCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAg                                                                                                                                                         >  2:1685393/1‑150 (MQ=255)
  cTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGAt                                                                                                                                                       >  1:208777/1‑150 (MQ=255)
        cGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACAGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                                                                                 >  2:1082087/1‑150 (MQ=255)
         gTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                                                                                >  1:802938/1‑150 (MQ=255)
            aaCATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATggct                                                                                                                                             <  2:1115507/150‑3 (MQ=255)
             aCATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑ttt                                                                                                                                          >  2:228741/1‑150 (MQ=255)
               aTTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGG‑‑ctttt                                                                                                                                        <  1:800812/150‑1 (MQ=255)
                 tGCGAAGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGGGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑‑cttttca                                                                                                                                      <  2:1133855/150‑3 (MQ=255)
                   cGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGCCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAGGG‑‑ttttttttt                                                                                                                                    >  2:204391/1‑150 (MQ=255)
                    taTGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGAGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCATATTTGGATGG‑‑TTTTTTTTTg                                                                                                                                   <  1:743621/149‑1 (MQ=255)
                    gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTg                                                                                                                                   >  1:368235/1‑150 (MQ=255)
                    gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTg                                                                                                                                   >  1:498274/1‑150 (MQ=255)
                    gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGGTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTg                                                                                                                                   <  2:895330/150‑1 (MQ=255)
                     aTGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑ttttttttttg                                                                                                                                  <  2:386008/150‑3 (MQ=255)
                     aTGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGACATCCAGATTTGGATGG‑‑TTTTTTTTTGt                                                                                                                                  >  2:1441982/1‑150 (MQ=255)
                      tgcgTTGATCAGCATCCCGATGGCTGGAAGCGTTGCTTCCCTGACCGTTTCGGTGGCGACCGGAATTTGCTTATTGGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTc                                                                                                                                 <  1:1221531/147‑1 (MQ=255)
                       gAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCATCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑ttttttttttgtc                                                                                                                                <  2:776613/150‑5 (MQ=255)
                       gAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTATTTTTTGTCt                                                                                                                                >  1:857053/1‑150 (MQ=255)
                       aaGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCt                                                                                                                                >  1:1362752/2‑150 (MQ=255)
                        aGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTc                                                                                                                                 >  2:756023/1‑148 (MQ=255)
                        aGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCta                                                                                                                               >  1:139030/1‑150 (MQ=255)
                          ttGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCtata                                                                                                                             >  1:1605907/1‑150 (MQ=255)
                          ttGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCtata                                                                                                                             <  1:487895/150‑1 (MQ=255)
                          ttGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCtata                                                                                                                             <  2:221093/150‑1 (MQ=255)
                           tGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAg                                                                                                                            >  1:873806/1‑150 (MQ=255)
                            gATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAg                                                                                                                            >  2:353927/1‑150 (MQ=255)
                            gATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGc                                                                                                                           >  2:343521/1‑150 (MQ=255)
                              tCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTg                                                                                                                         >  2:599637/1‑150 (MQ=255)
                               cAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTgg                                                                                                                        <  1:1281573/150‑1 (MQ=255)
                                 ccATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGCAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGt                                                                                                                       <  1:1710025/149‑1 (MQ=255)
                                     cccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGTAAGa                                                                                                                   <  1:353927/150‑1 (MQ=255)
                                       cGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATaa                                                                                                                >  2:461915/1‑150 (MQ=255)
                                          tGGCTGAAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGACGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                                                             <  2:549182/150‑1 (MQ=255)
                                             cTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                                                                          <  1:226094/150‑1 (MQ=255)
                                              tGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                                                                         >  1:299258/1‑150 (MQ=255)
                                                gAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                                                                       <  2:138996/150‑1 (MQ=255)
                                                  aGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTAAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                                                                     <  2:171126/150‑1 (MQ=255)
                                                   gCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                                                    >  1:938767/1‑150 (MQ=255)
                                                     gTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATa                                                                                                                 <  1:935178/135‑1 (MQ=255)
                                                      tttCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                                                                 >  1:66799/1‑150 (MQ=255)
                                                      tttCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAat                                                                                                 <  2:299258/150‑1 (MQ=255)
                                                       ttcttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                                                                <  1:978103/150‑1 (MQ=255)
                                                         cttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAatat                                                                                               >  1:1364700/1‑150 (MQ=255)
                                                          ttcCCTGAATGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCAATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                                                             <  2:984511/150‑1 (MQ=255)
                                                            cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                                                           >  2:614199/1‑150 (MQ=255)
                                                            cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                                                           >  1:1135700/1‑150 (MQ=255)
                                                             ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                          >  2:666013/1‑150 (MQ=255)
                                                             ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                          >  1:1393870/1‑150 (MQ=255)
                                                              cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                                                           >  2:1036977/1‑150 (MQ=255)
                                                              cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                                                         <  1:599637/150‑1 (MQ=255)
                                                               tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                                                                        >  1:167341/1‑150 (MQ=255)
                                                               tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtttt                                                                                        >  1:1378293/1‑150 (MQ=255)
                                                                      ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAata                                                                                 <  1:1082087/150‑1 (MQ=255)
                                                                         ggTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAttt                                                                              >  1:417311/1‑150 (MQ=255)
                                                                              cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                                                         >  1:1108957/1‑150 (MQ=255)
                                                                                    gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                                                                   <  2:368235/150‑1 (MQ=255)
                                                                                    gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                                                                   <  2:498274/150‑1 (MQ=255)
                                                                                         tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATCATGAAAAGCCATCCAGATTTGGATGG‑‑GTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatatt                                                              >  2:1066001/1‑149 (MQ=255)
                                                                                          gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                                                             >  1:224099/1‑150 (MQ=255)
                                                                                          gCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                                                             <  1:925280/150‑1 (MQ=255)
                                                                                             tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                                                          >  2:451842/1‑150 (MQ=255)
                                                                                             tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                                                          <  2:246734/150‑1 (MQ=255)
                                                                                             tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat                                                          <  2:1506762/150‑1 (MQ=255)
                                                                                                  gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat                                                     <  2:1393870/150‑1 (MQ=255)
                                                                                                  gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat                                                     >  2:873014/1‑150 (MQ=255)
                                                                                                   aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata                                                    >  2:317222/1‑150 (MQ=255)
                                                                                                   aaGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTATTTTTTGTCTATAGCTGGAAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata                                                    <  2:857053/150‑1 (MQ=255)
                                                                                                    aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGAAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat                                                   >  2:1360032/1‑150 (MQ=255)
                                                                                                       ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTATTCTATATGGAATAATATATGAc                                                <  2:333579/150‑1 (MQ=255)
                                                                                                         tGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAt                                              <  2:1179657/150‑1 (MQ=255)
                                                                                                          gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg                                             <  1:1720800/150‑1 (MQ=255)
                                                                                                           cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATgg                                            <  1:1049317/150‑1 (MQ=255)
                                                                                                               aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATc                                        <  2:886516/150‑1 (MQ=255)
                                                                                                               aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATc                                        <  2:1459588/150‑1 (MQ=255)
                                                                                                                gcgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATcc                                       >  2:1228607/1‑150 (MQ=255)
                                                                                                                       tAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                  <  1:1166670/150‑1 (MQ=255)
                                                                                                                        aaTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGgcgc                               <  2:1319407/150‑1 (MQ=255)
                                                                                                                         aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGATTAATATATGACATGGAATCCGTTGGCGCt                              >  1:1052500/1‑150 (MQ=255)
                                                                                                                           ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAg                            <  2:618958/150‑1 (MQ=255)
                                                                                                                           ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAg                            <  1:666013/150‑1 (MQ=255)
                                                                                                                            tctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGc                           >  2:1671489/1‑150 (MQ=255)
                                                                                                                            tctcAGAATTATGAAAAGCCATACAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGc                           <  2:802938/150‑1 (MQ=255)
                                                                                                                             ctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCg                          >  2:1016071/1‑150 (MQ=255)
                                                                                                                               cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAAAATATTACATTGAATCCGTTGGCGCTAGCGa                         >  2:1454098/1‑150 (MQ=255)
                                                                                                                               cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGAc                        >  1:1478191/1‑150 (MQ=255)
                                                                                                                               cAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATAAGTTGGCTATAGCTAc                        >  1:240107/1‑150 (MQ=255)
                                                                                                                                 gAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACgg                      >  2:965102/1‑150 (MQ=255)
                                                                                                                                   acaaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgc                     >  2:595298/7‑150 (MQ=255)
                                                                                                                                    ttATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcgc                   >  1:451565/1‑150 (MQ=255)
                                                                                                                                        gAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCa               <  1:873014/150‑1 (MQ=255)
                                                                                                                                        aaaaaACCATCCAAATCTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCa               >  2:531499/2‑150 (MQ=255)
                                                                                                                                         aaaaACCATCCAAATCTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATATATATGGAATAATATATGACATGGAATCCGTTGTCGCTAGCGACGGCGCTGCaa              >  1:178032/1‑150 (MQ=255)
                                                                                                                                           aaGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAAc            <  1:1671489/150‑1 (MQ=255)
                                                                                                                                           aaGACATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAAc            <  1:888975/150‑1 (MQ=255)
                                                                                                                                            aGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACt           >  2:926392/1‑150 (MQ=255)
                                                                                                                                                aTCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTTAAATTTAAAGTTAATATATATGGAATAATATATGACATGGAATCCGTTGGAGCTAGCGACGGCGCTGCAAACTGTAc       >  1:1306082/1‑150 (MQ=255)
                                                                                                                                                 tCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAcc      >  2:481443/1‑150 (MQ=255)
                                                                                                                                                  ccAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCt     <  1:461915/150‑1 (MQ=255)
                                                                                                                                                   cAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGCTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATTGAATCAGTAGTCGATATCGACGGCGCTTCAAACTGTACCTg    >  1:293670/1‑150 (MQ=255)
                                                                                                                                                     gATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGaa  >  2:1613342/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAA  >  NC_000913/4409885‑4410183

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: