Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,326,782 G→C 5.7% intergenic (+141/‑70) yciQ → / → rluB DUF2207 domain‑containing protein YciQ/23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,326,7820GC 5.7% 286.3 / 4.0 84intergenic (+141/‑70)yciQ/rluBDUF2207 domain‑containing protein YciQ/23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base G (48/31);  new base C (3/2);  total (51/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGGTAATTAAGCTGATGTTAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAA  >  NC_000913/1326633‑1326926
                                                                                                                                                     |                                                                                                                                                
tgCTGGTAATTAAGCTGATGTTAATCGAAACAAGCTTTAACAATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATg                                                                                                                                                  <  1:1750266/149‑1 (MQ=255)
  cTGGTAATTAAGCTGATGTTAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTa                                                                                                                                                >  1:1686428/1‑150 (MQ=255)
      tAATTAAGCTGATGTTAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGtt                                                                                                                                            >  2:468583/1‑150 (MQ=255)
        aTTAAGCTGATGTTAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCt                                                                                                                                          <  1:679291/150‑1 (MQ=255)
           aaGCTGATGTTAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGACTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGAc                                                                                                                                       >  1:641291/1‑150 (MQ=255)
                  tGTTAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGt                                                                                                                                >  1:1722388/1‑150 (MQ=255)
                    ttAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTg                                                                                                                              <  1:1252328/150‑1 (MQ=255)
                    ttAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTg                                                                                                                              >  2:246383/1‑150 (MQ=255)
                    ttAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGAACTTACATCTGACTCACTgata                                                                                                                              <  1:220387/150‑4 (MQ=255)
                     tAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGc                                                                                                                             >  1:1843299/1‑150 (MQ=255)
                          gAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACaaa                                                                                                                        >  2:1351968/1‑150 (MQ=255)
                            aaCAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAg                                                                                                                      <  1:1351968/150‑1 (MQ=255)
                                aGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCata                                                                                                                  >  2:814407/1‑150 (MQ=255)
                                aGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCata                                                                                                                  <  1:1545531/150‑1 (MQ=255)
                                    ttAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAg                                                                                                              <  1:306454/150‑1 (MQ=255)
                                     tAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGa                                                                                                             <  1:145313/150‑1 (MQ=255)
                                       aCCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACg                                                                                                           >  2:533596/1‑150 (MQ=255)
                                        ccATAAAGATAGTCTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTAAACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATCTAAGTTCTGACTCACTGTTTCACAAAAGCATATCAGACGc                                                                                                          >  1:1776123/1‑150 (MQ=255)
                                             aaGATAGTGTAGTATGTTGCGCCTCAACGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGGCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTggg                                                                                                     <  1:1001101/150‑1 (MQ=255)
                                              aGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGa                                                                                                    >  2:283047/1‑150 (MQ=255)
                                              aGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGAACTTACATCTGACTCACt                                                                                                                                  <  1:1865023/120‑1 (MQ=39)
                                              aGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGAACTTACATCTGACTCACt                                                                                                                                  >  2:1865023/1‑120 (MQ=39)
                                                     gtagtaTGTTGCGCCTCAAAGCCAGGCCGATAAACGGCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTCCTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGac                                                                                             <  1:1875892/150‑1 (MQ=255)
                                                         taTGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCa                                                                                         <  2:1843299/150‑1 (MQ=255)
                                                           tGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCaaa                                                                                       >  1:344494/1‑150 (MQ=255)
                                                             ttGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAg                                                                                     <  1:504930/150‑1 (MQ=255)
                                                                cgcCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGCAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGaa                                                                                  >  1:1292853/1‑150 (MQ=255)
                                                                cgcCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGaa                                                                                  >  1:1962424/1‑150 (MQ=255)
                                                                  ccTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGc                                                                                <  2:1686428/150‑1 (MQ=255)
                                                                  ccTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGc                                                                                >  1:343894/1‑150 (MQ=255)
                                                                     cAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGACTGGGAAGGCGACACCAAAAGGAAGctct                                                                             >  1:1463003/1‑150 (MQ=255)
                                                                      aaaGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTa                                                                            <  2:1594121/150‑1 (MQ=255)
                                                                       aaGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTAt                                                                           >  1:625569/1‑150 (MQ=255)
                                                                          ccAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAg                                                                        >  2:355159/1‑150 (MQ=255)
                                                                          ccAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAg                                                                        >  2:1712736/1‑150 (MQ=255)
                                                                           cAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGc                                                                       <  1:1793314/150‑1 (MQ=255)
                                                                           cAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGc                                                                       >  1:576128/1‑150 (MQ=255)
                                                                               ccGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGaaa                                                                   >  2:1205355/1‑150 (MQ=255)
                                                                                cGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGaaaa                                                                  <  1:533596/150‑1 (MQ=255)
                                                                                   tAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCt                                                               >  1:32300/1‑150 (MQ=255)
                                                                                   tAAACGTCGAGTCGTTTACTTAAGGCCGGAAGAGTTCAAACAGACTTTACTATCAGTGAGTAAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACCCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCt                                                               >  2:494593/1‑150 (MQ=255)
                                                                                     aaCGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTAc                                                             >  2:962980/1‑150 (MQ=255)
                                                                                     aaCGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTAc                                                             <  1:136175/150‑1 (MQ=255)
                                                                                        gTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGAc                                                                                                            <  1:796655/100‑1 (MQ=255)
                                                                                        gTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGAc                                                                                                            >  2:796655/1‑100 (MQ=255)
                                                                                             gTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTg                                                     >  2:729339/1‑150 (MQ=255)
                                                                                               cGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCt                                                   >  2:1020232/1‑150 (MQ=255)
                                                                                                 tttACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTgg                                                 <  2:625569/150‑1 (MQ=255)
                                                                                                          ggCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCg                                        >  1:576036/1‑150 (MQ=255)
                                                                                                            ccTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGc                                                                <  1:1699671/124‑1 (MQ=255)
                                                                                                            ccTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGc                                                                >  2:1699671/1‑124 (MQ=255)
                                                                                                            ccTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGc                                      >  2:1656260/1‑150 (MQ=255)
                                                                                                             cTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGcc                                     <  1:917512/150‑1 (MQ=255)
                                                                                                               gAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCAc                                   >  1:1926188/1‑150 (MQ=255)
                                                                                                               gAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCAc                                   >  2:735858/1‑150 (MQ=255)
                                                                                                                 agagTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACgg                                 <  1:551939/150‑1 (MQ=255)
                                                                                                                    gTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTtc                              <  2:47538/150‑1 (MQ=255)
                                                                                                                     ttCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTtct                             <  1:1446838/150‑1 (MQ=255)
                                                                                                                       cAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGTTGGCGCGTGCCGGCCACGGTTCTCg                           >  1:1461025/1‑150 (MQ=255)
                                                                                                                       cAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCg                           <  1:729339/150‑1 (MQ=255)
                                                                                                                        aaaCAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGc                          >  2:300800/1‑150 (MQ=255)
                                                                                                                         aaCAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGcc                         <  2:1722388/150‑1 (MQ=255)
                                                                                                                           cAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGt                       >  1:140338/1‑150 (MQ=255)
                                                                                                                            aGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTg                      <  2:576036/150‑1 (MQ=255)
                                                                                                                             gACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGa                     >  1:258998/1‑150 (MQ=255)
                                                                                                                               cTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGaaa                   >  2:538924/1‑150 (MQ=255)
                                                                                                                                  tACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCg                >  2:78732/1‑150 (MQ=255)
                                                                                                                                     tATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAAt             >  1:462221/1‑150 (MQ=255)
                                                                                                                                      aTCAGTGAGTCAGATGTAAGTTCTGACTCACTGATGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAAta            >  2:249863/1‑149 (MQ=255)
                                                                                                                                       tCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCt           >  2:1984169/1‑150 (MQ=255)
                                                                                                                                          gTGAGTCAGAACTTACATCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCtatt        >  2:1480000/1‑150 (MQ=255)
                                                                                                                                            gAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCtattat      <  1:1916391/150‑1 (MQ=255)
                                                                                                                                            gAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCtattat      >  1:297691/1‑150 (MQ=255)
                                                                                                                                             aGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTattatt     <  1:1003704/150‑1 (MQ=255)
                                                                                                                                              gTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTg    >  1:500711/1‑150 (MQ=255)
                                                                                                                                              gTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTg    >  2:1037154/1‑150 (MQ=255)
                                                                                                                                              gTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTg    >  1:901257/1‑150 (MQ=255)
                                                                                                                                              gTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTg    >  1:841157/1‑150 (MQ=255)
                                                                                                                                               tCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGa   <  2:344494/150‑1 (MQ=255)
                                                                                                                                               tCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGa   <  2:483300/150‑1 (MQ=255)
                                                                                                                                               tCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGa   <  1:283047/150‑1 (MQ=255)
                                                                                                                                                cAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGc                                                                       <  2:1316568/81‑1 (MQ=255)
                                                                                                                                                cAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGc                                                                       >  1:1316568/1‑81 (MQ=255)
                                                                                                                                                cAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGaa  >  1:1375393/1‑150 (MQ=255)
                                                                                                                                                cAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGaa  >  2:1381036/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GGCTGGTAATTAAGCTGATGTTAATCGAAACAAGCTTTAACCATAAAGATAGTGTAGTATGTTGCGCCTCAAAGCCAGGCCGATAAACGTCGAGTCGTTTACTTAAGGCCTGAAGAGTTCAAACAGACTTTACTATCAGTGAGTCAGATGTAAGTTCTGACTCACTGTTGCACAAAAGCATATCAGACGCCTGGGAAGGCGACACCAAAAGGAAGCTCTATGAGCGAAAAGCTACAGAAAGTGCTGGCGCGTGCCGGCCACGGTTCTCGCCGTGAAATCGAATCTATTATTGAA  >  NC_000913/1326633‑1326926

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