Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,034 G→A 5.1% 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 5.1% 416.5 / 7.5 139intergenic (+28/‑99)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base G (62/67);  new base A (4/3);  total (66/73)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.51e-01

Only 100 of 139 total aligned reads displayed.

ACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTG  >  NC_000913/4409890‑4410181
                                                                                                                                                |                                                                                                                                                    
aCTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTg                                                                                                                                                 <  2:15449/150‑1 (MQ=255)
   cGTTAACTTTGCGATGAGTTGATCAGCATCCCGATGGCTGGTAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                                                                              <  1:437691/150‑1 (MQ=255)
   cGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGAt                                                                                                                                              >  2:1102558/1‑150 (MQ=255)
    gTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATg                                                                                                                                             <  2:1952210/150‑1 (MQ=255)
          aTTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGG‑ctttt                                                                                                                                      <  1:394040/150‑1 (MQ=255)
           ttGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAAATCTGGATGG‑cttttc                                                                                                                                     <  2:1886916/150‑2 (MQ=255)
             gCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑tttttttt                                                                                                                                   >  2:1737739/1‑150 (MQ=255)
             gCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑tttttttt                                                                                                                                   >  1:1526786/1‑150 (MQ=255)
             gCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGGTTCTTCCCTGAACGTTTCGGGTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑tttttttt                                                                                                                                   <  2:1206151/150‑1 (MQ=255)
              ccaTGATTTGATCAGCATCCCGATGGCTTGAAGCGTTTCTTCCCTGCACGTTTCGGTGCCGACCGGAATTTGCTTATTTGTAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑ttttttttt                                                                                                                                  <  2:92017/148‑1 (MQ=255)
               gATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTg                                                                                                                                 <  1:1882040/150‑1 (MQ=255)
                 tGAGTTGATCAGCATCCCGATGGCTGTAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTc                                                                                                                               <  1:880393/150‑1 (MQ=255)
                      tGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAg                                                                                                                          >  1:1056465/1‑150 (MQ=255)
                      tGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAg                                                                                                                          >  1:1063384/1‑150 (MQ=255)
                       gATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGc                                                                                                                         >  1:720396/1‑150 (MQ=255)
                       gATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGc                                                                                                                         >  1:662216/1‑150 (MQ=255)
                           aGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGt                                                                                                                     >  1:1066549/1‑150 (MQ=255)
                             cATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTaa                                                                                                                   >  2:463397/1‑150 (MQ=255)
                              aTCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAg                                                                                                                  <  2:1288248/150‑1 (MQ=255)
                              aTCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAg                                                                                                                  >  1:314105/1‑150 (MQ=255)
                                 ccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATa                                                                                                               <  2:1578117/150‑1 (MQ=255)
                                 ccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATa                                                                                                               >  2:688285/1‑150 (MQ=255)
                                   gATGGCAGGAAGCGTTTATACCCTCAACGTTACGGATGCGCCCGGAATATGCTTATTTGAAGCGCTGCGCCAGAGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTTCTAAGATaaa                                                                                                             >  2:1451447/1‑149 (MQ=255)
                                     tGGCTTGAATCGTTTCTTTCCTGAACGTTTCGGTTTCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGGAGCTAATTTCTCAGAATTATGTAAAGCCATCGAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                                                           <  2:1905915/150‑1 (MQ=255)
                                          ggAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTTAGATAATTACGTAt                                                                                                      >  2:1497457/1‑150 (MQ=255)
                                           gAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                                                                     <  2:13700/150‑1 (MQ=255)
                                           gAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                                                                     <  2:908506/150‑1 (MQ=255)
                                             aGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGc                                                                                                   <  1:1418255/150‑1 (MQ=255)
                                              gCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                                                  <  2:490216/150‑1 (MQ=255)
                                               gggttcttcCCAGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGGGCGCCAGCGCATCTAATTTCTCAGAATTATGAAAACCCATCCAGATTTGGATGGTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                                                                  <  2:798033/147‑1 (MQ=255)
                                                gTTTCTTCCCTGAATGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                                                                <  1:509454/150‑1 (MQ=255)
                                                    cttcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                                                            <  1:1674089/150‑1 (MQ=255)
                                                      tcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                                                          <  2:67137/150‑1 (MQ=255)
                                                      tcCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                                                          >  2:1636983/1‑150 (MQ=255)
                                                       cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                                                         >  2:1898554/1‑150 (MQ=255)
                                                       cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAt                                                                                         >  2:83247/1‑150 (MQ=255)
                                                        ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                        >  2:1350294/1‑150 (MQ=255)
                                                        ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                        <  2:435466/150‑1 (MQ=255)
                                                        ccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATCG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                                                                        >  2:1582622/1‑150 (MQ=255)
                                                         cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                                                                       <  2:73258/150‑1 (MQ=255)
                                                             aCGTTTCGGTTGCGACCGGAATTTGCTTATTTTAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                                                                   <  2:1027008/150‑1 (MQ=255)
                                                              cGTTTGGGTTGCGAACGGAATTTGCTTATTTGAAGCGGTGAGCCAGCGCAGCTAATTTCTCATAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                                                                  <  2:781900/150‑1 (MQ=255)
                                                              cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                                                                  <  2:675589/150‑1 (MQ=255)
                                                              cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                                                                  >  1:48559/1‑150 (MQ=255)
                                                              cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTttatta                                                                                  >  2:501348/1‑150 (MQ=255)
                                                                   gggtgccGACCGTAATTTGCTTATTTTACGCGGTGCGCAAGGGATGCTAATTACTCAGAATTATTAAATGCCATCCTCATTTGCATGG‑‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATAAAATGAATAttt                                                                            <  2:1321185/144‑1 (MQ=255)
                                                                   cGGTTTCTACCGGAATTTGCTTATTTTAAGCGGTGAGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                                                             <  1:1924558/150‑1 (MQ=255)
                                                                   cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                                                             >  2:1182264/1‑150 (MQ=255)
                                                                     gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                                                           >  2:1515115/1‑150 (MQ=255)
                                                                       tGCGACCGGAAATTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaaa                                                                         <  2:1783739/150‑1 (MQ=255)
                                                                        gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                                                                        >  1:320928/1‑150 (MQ=255)
                                                                        gCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                                                       >  2:1441360/1‑150 (MQ=255)
                                                                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                                                       <  2:1391841/150‑1 (MQ=255)
                                                                         cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAAATATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTACGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                                                                       <  1:302188/150‑1 (MQ=255)
                                                                          gcccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTGATTTTATTAATATTTAAACtt                                                                      >  1:1351669/3‑150 (MQ=255)
                                                                          gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                                                                      <  2:33024/150‑1 (MQ=255)
                                                                          gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACtt                                                                      >  2:1542708/1‑150 (MQ=255)
                                                                            ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGAAAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTaa                                                                    >  2:572013/1‑150 (MQ=255)
                                                                              ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                                                                  <  2:1833241/150‑1 (MQ=255)
                                                                                    tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                                                            >  1:1734354/1‑150 (MQ=255)
                                                                                    tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                                                            >  1:1726309/1‑150 (MQ=255)
                                                                                         aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaata                                                       <  2:885791/150‑1 (MQ=255)
                                                                                         aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaata                                                       <  1:1528421/150‑1 (MQ=255)
                                                                                            tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata                                                    >  1:802951/1‑150 (MQ=255)
                                                                                            tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata                                                    >  1:1393005/1‑150 (MQ=255)
                                                                                            tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata                                                    >  2:1837121/1‑150 (MQ=255)
                                                                                            tGAAGAGGTGCGCCAGCTCAGCTAATTTCTCAGAATTATTAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata                                                    <  2:72466/150‑1 (MQ=255)
                                                                                             gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataataa                                                   >  2:885257/1‑149 (MQ=255)
                                                                                             gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat                                                   <  1:1542708/150‑1 (MQ=255)
                                                                                                  ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATATCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc                                              <  2:153890/150‑1 (MQ=255)
                                                                                                   gTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACa                                             <  2:1021184/150‑1 (MQ=255)
                                                                                                     gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTTTATGGAATAATATATGACATg                                           >  2:427581/1‑150 (MQ=255)
                                                                                                     gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg                                           >  2:1040460/1‑150 (MQ=255)
                                                                                                        ccAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGaa                                        >  1:942769/1‑150 (MQ=255)
                                                                                                          aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATc                                      <  1:1577085/150‑1 (MQ=255)
                                                                                                          aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATc                                      >  2:1621726/1‑150 (MQ=255)
                                                                                                           gcgcAGCTAATTTCTCATAATTATGAAAATCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATcc                                     <  1:529048/150‑1 (MQ=255)
                                                                                                            cgcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCg                                    <  2:1105951/150‑1 (MQ=255)
                                                                                                             gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGt                                   >  1:38360/1‑150 (MQ=255)
                                                                                                             gcAGCTAATTTCTCAGAATTATGAAAAGCAATCCTGATTTGTATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGt                                   <  1:1803067/150‑1 (MQ=255)
                                                                                                                gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTgg                                <  2:1593554/150‑1 (MQ=255)
                                                                                                                      ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTAAAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAg                          <  1:449845/150‑1 (MQ=255)
                                                                                                                        ctcAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCg                        <  1:929915/150‑1 (MQ=255)
                                                                                                                        ctcAGAATTATGAAAAGCCATACAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCg                        >  1:1491078/1‑150 (MQ=255)
                                                                                                                            gAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATATATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACgg                    >  1:445462/1‑150 (MQ=255)
                                                                                                                            gAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAAGATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACgg                    >  1:1647244/1‑150 (MQ=255)
                                                                                                                             aaTTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgc                   >  1:623893/1‑150 (MQ=255)
                                                                                                                               ttATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGgcgc                 <  2:1953601/150‑1 (MQ=255)
                                                                                                                               caaaaaaaaACCATCCAAATCTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATTGAATAATATATGACATGTAATCCTTTGGCGCTAGCGACGgcga                 >  1:123415/6‑149 (MQ=255)
                                                                                                                                 aaaaaaaACCATCCAAATCTGGATGG‑CTTTTTTTTGTCTATAGATGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCTACGGCTCTg               >  2:755303/4‑150 (MQ=255)
                                                                                                                                  tGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGc              >  2:1761740/1‑150 (MQ=255)
                                                                                                                                   taaaaGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCa             <  1:1308729/149‑1 (MQ=255)
                                                                                                                                      aaaccATCCAAATCTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTAATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATTACATGGAATCCGTTGGATCTAGCGACGGCGCTGCaaaa          >  1:90466/4‑149 (MQ=255)
                                                                                                                                          cATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTa      <  1:74007/150‑1 (MQ=255)
                                                                                                                                           aTCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTAc     <  2:810298/150‑1 (MQ=255)
                                                                                                                                             ccAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCt   <  2:600392/150‑1 (MQ=255)
                                                                                                                                             ccAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCt   >  2:1165052/1‑150 (MQ=255)
                                                                                                                                             ccAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCt   <  1:1544346/150‑1 (MQ=255)
                                                                                                                                              cAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTg  <  1:771099/150‑1 (MQ=255)
                                                                                                                                              cAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTg  >  2:927181/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTG  >  NC_000913/4409890‑4410181

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