Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,026 G→A 7.4% intergenic (+20/‑107) 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,0260GA 7.4% 155.3 / 2.3 54intergenic (+20/‑107)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base G (30/20);  new base A (2/2);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGTGAGTAGTATCAGCAAT  >  NC_000913/4409738‑4410324
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                            
cTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGCAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTTGACGCGTTTCTTCCATGAACGTTTCGGTTGCGCCCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAACTAATTTATCAGAATTATGAAACGCCATCCAGAgtt                                                                                                                                                                                                                                                                                                  >  1:217803/1‑301 (MQ=255)
      tttggcgcACCATGCGTTTGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATgg                                                                                                                                                                                                                                                                                            <  1:147064/297‑1 (MQ=255)
       tGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑t                                                                                                                                                                                                                                                                                         >  2:519157/1‑301 (MQ=255)
                         gTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCt                                                                                                                                                                                                                                                                       >  2:191401/1‑301 (MQ=255)
                            gCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGCTCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGCCCGGACTTTGCTTATTTTAAGCGGTGCGCCAGCGCAGCTAATTTATCAGAATTCTGAAAAGCCATCCAGATTTGGCTGG‑‑TTTTTTTTTGTCTATCGCTgtt                                                                                                                                                                                                                                                                    >  1:602491/1‑299 (MQ=255)
                            gCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTTCGCCAGCGCAGCTAATTTCTCAGCATTATGAAAAGCCATCCAGCTTTGGCAGG‑‑TTTTTTTTTTTCTAAAGCTGGt                                                                                                                                                                                                                                                                    >  1:94704/1‑301 (MQ=255)
                               ggAAGGGAATATCGGGGTCGTCGGTACGGCAGGCGAGGCGGTTCATACACTCTATCAGAGCAAAGTGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAg                                                                                                                                                                                                                                                                 <  1:359929/301‑1 (MQ=255)
                                gAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTNNNNAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTNNCGGTTGCGACCGGAANNNNCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAAAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCtaacactggaaagt                                                                                                                                                                                                                                                                >  1:321132/1‑289 (MQ=255)
                                 aaGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTGTCTATAGCTGGTAAGATa                                                                                                                                                                                                                                                              <  2:170917/301‑1 (MQ=255)
                                  agagAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCTCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATa                                                                                                                                                                                                                                                              <  1:266082/301‑1 (MQ=255)
                                    tcagacgtgtgctcttccgatctgCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                                                                                                                                                                                                            <  2:166485/278‑1 (MQ=255)
                                     gattatTTGGATGGTCGGTACGGCAGGCGAGGGGGATCATACACTCGATCAGGGCAAAATGACCGGACGTCTGGCGCTGGGGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                                                                                                                                                                                                           <  1:238419/298‑1 (MQ=255)
                                            tGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGACCATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCTGACTTTGCTTCTTTGAAGCGGTGCGCCACCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGG‑‑TTTTTTTTTGTCTAAAGCTGGTAAGATAATTACATAtt                                                                                                                                                                                                                                                    >  1:522646/1‑301 (MQ=255)
                                                 acgacgctgttccgttctggcgGTTCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaa                                                                                                                                                                                                                                               <  1:504930/283‑1 (MQ=255)
                                                    cgGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAataa                                                                                                                                                                                                                                            >  2:228531/1‑300 (MQ=255)
                                                      gTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGATAAGATAATTACGTATTGAAAATAttt                                                                                                                                                                                                                                          >  1:135846/1‑300 (MQ=255)
                                                           gCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAAtagaacggaagagcgacgagcagg                                                                                                                                                                                                                                     >  1:166288/1‑278 (MQ=255)
                                                               ccgatctggATCATACACTCTTTCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAGCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatta                                                                                                                                                                                                                                 <  1:276012/294‑1 (MQ=255)
                                                                   ggcgGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCaaaagatcggaagagcacaca                                                                                                                                                                                                                             >  2:505602/1‑283 (MQ=255)
                                                                      ggATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAGCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAagatcgg                                                                                                                                                                                                                          >  2:276315/1‑298 (MQ=255)
                                                                       taTCATACACTCTATCAGAGCAACATGACCGGACGTCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTaa                                                                                                                                                                                                                         <  2:8433/300‑1 (MQ=255)
                                                                                    aTCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAACAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                                                                                                                                                                                                            >  2:378111/1‑301 (MQ=255)
                                                                                      caggggaaagtggcgggtcGCCGGGGGCTGGTGATGGGTGCGGAAGTGGAAGGTATGCTTCGCCTGACTCGGGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                                                                                                                                                                                                          <  1:518485/283‑1 (MQ=255)
                                                                                              aaTGACCGGACGCCTGGCGTTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGACTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatat                                                                                                                                                                                                  <  2:135971/301‑1 (MQ=255)
                                                                                                           ttgtcgttGGTTATGGGTGCGGAAGGTGAAGGTTTGCGTCGCCTGACTCGTGAACATTGCGATGAGTTTATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATc                                                                                                                                                                                     <  2:218046/294‑1 (MQ=255)
                                                                                                           cTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCGGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGAAATAATATATGAAATGGaaat                                                                                                                                                                                     >  1:460318/1‑299 (MQ=255)
                                                                                                                 cTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATAGAATCAGTTga                                                                                                                                                                               >  2:16458/1‑300 (MQ=255)
                                                                                                                                               cGTTGCCTGACTTGTGGTCATTGTGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTg                                                                                                                                                 <  1:228288/301‑1 (MQ=255)
                                                                                                                                                           cGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCAAGCGACGGAGCTGAAAACTGTACCTGAAAAAAAAAttt                                                                                                                                     >  1:248185/1‑300 (MQ=255)
                                                                                                                                                              ttaCATTGTGATTAGTTTATCAGCGTCCTGGTGGTTGGAAGCGTTTCTTCGTTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGGGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATg                                                                                                                                  <  2:321536/299‑1 (MQ=255)
                                                                                                                                                              gAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATg                                                                                                                                  <  2:73125/301‑1 (MQ=255)
                                                                                                                                                                  tttGGGATGGGTTGATAAGCATCCTGATGGGTGGAAGCGTTTTTTTACTGAACGTGTCGGTTGCGACTGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAAc                                                                                                                              <  2:94828/300‑1 (MQ=255)
                                                                                                                                                                        tttaGTTGTTCTGCTTTCCTATTGGTGGAAGGGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCGGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAg                                                                                                                        <  2:460882/298‑1 (MQ=255)
                                                                                                                                                                              tGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGCTTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGaa                                                                                                                  >  2:537483/1‑301 (MQ=255)
                                                                                                                                                                              tGATCAGCATCCCGATGGCTGGAAGCGTTTATTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTCAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGAAACAAATaacaaaaa                                                                                                                  >  2:255468/1‑294 (MQ=255)
                                                                                                                                                                                        cccGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCattaata                                                                                                        >  2:643393/1‑300 (MQ=255)
                                                                                                                                                                                               gCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAAGTTCTCAGAATTATGAGAAGCCATCCAGATTCGGATGG‑‑TTTTTTTTNNNNNATAGCAGGTAAGATNNNNACGTATNNNNNNTATTCATTTNNNNAATATTTAAANNNNNNCTATATGGAATAANNNNTAACATGGAANCCGTTGGCGCTNNNNNNNGCGATGCAAACAGTACCTGAAAACAAACTTAATGTAACAAAAAACAAGAACACAATAAATATTaaaa                                                                                                 >  1:250987/1‑301 (MQ=255)
                                                                                                                                                                                                               cccTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAACAATTGATGTAACAAATCGCAAAAACGCAATAATTATTAAAATAAATAATAATGGCGa                                                                                 >  1:61344/1‑301 (MQ=255)
                                                                                                                                                                                                                           tggtttctattgtaaTTTGCTTTTTTTAAGCGGTGCGGCAGCGCAGCTAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCaa                                                                     <  2:248452/288‑1 (MQ=255)
                                                                                                                                                                                                                                 cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTGACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCt                                                               >  2:533828/1‑301 (MQ=255)
                                                                                                                                                                                                                                                   ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGACCGCATTAATTATTCAAATGAATGATTATGGAGATTTGCAAATAAATATTCTTTTTACTTCCCGCCAAAt                                             >  1:408621/1‑301 (MQ=255)
                                                                                                                                                                                                                                                    tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAAta                                            >  2:251702/1‑300 (MQ=255)
                                                                                                                                                                                                                                                          ggTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATc                                      >  2:41943/1‑301 (MQ=255)
                                                                                                                                                                                                                                                             gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCTATTTGCAAATAAATATTATTTTTACTTCCCGCCCAATGATTATCGaa                                   >  1:376214/1‑301 (MQ=255)
                                                                                                                                                                                                                                                              ttttttatttttttttttcctatcacaaaaaaaaACCATCCAAATCTGGATGGTTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGa                                    <  2:552179/271‑1 (MQ=255)
                                                                                                                                                                                                                                                              cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTACCAAATAACGAGAACGCATTAATTATTAAAATGAATGATTAAGGCGATTTGCAAATAAAAATTCTTTTTACTTCCCGCAAAATAAATATCaaaa                                  >  1:501310/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                tcagacgtgtgctcttccgatctcaaaaaaaaACCATCCAAATCTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGTTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAAcc                                <  2:23649/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCTAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTAtttt                        >  2:553639/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                              tttttaGAATTATTAAAAGCCATCCAGATTTTGATTG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTGACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGt                  <  1:533155/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     aaTTATGAAAAGTCTTTCTTATTTGGATGG‑‑TTTTTTTTTGTCTATAGTTGGTAAGATAATTACGTATTGCAAATATTTATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGTGagtagt           <  1:251435/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       caaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTAACaaaaatcggaaaagagtcgggtagggg           >  1:551506/6‑278 (MQ=255)
                                                                                                                                                                                                                                                                                       caaaaaaaaACCATCCAAATCTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGTTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAAccagatcggaagagcgtcgagtagg         >  1:23598/6‑278 (MQ=255)
                                                                                                                                                                                                                                                                                           gAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTTCGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTTTACCTGAACAAAATATTGCTGTAACAAATACCACGACCGCATTAATTCTTAAACTGAATGATTATGGCGATTTGCAAAACACTCTTCTTTTTACTTCCCGCCAAATAATTCACGAAACCTTTCTTTCTCCGGTGcgtagcacactc     >  1:456662/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                                              aaGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCCAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATATGTCCGGTGATTAGTATCAGCaaa  >  1:461233/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                            
CTAACCTTGCGCGCACCATGCGTATGTTGCAGGAAGAGAATATCTGGATCGTCGGTACGGCAGGCGAGGCGGATCATACACTCTATCAGAGCAAAATGACCGGACGCCTGGCGCTGGTGATGGGTGCGGAAGGTGAAGGTATGCGTCGCCTGACTCGTGAACATTGCGATGAGTTGATCAGCATCCCGATGGCTGGAAGCGTTTCTTCCCTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGG‑‑TTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATGGAATCCGTTGGCGCTAGCGACGGCGCTGCAAACTGTACCTGAACAAAATATTGATGTAACAAATAGCGAGAACGCATTAATTATTAAAATGAATGATTATGGCGATTTGCAAATCAATATTCTTTTTACTTCCCGCCAAATGATTATCGAAACCTTTATTTGTCCGGTGAGTAGTATCAGCAAT  >  NC_000913/4409738‑4410324

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: