Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,482,901 C→T V313M (GTG→ATG)  valS ← valine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,482,9010CT96.0% 227.5 / NA 75V313M (GTG→ATG) valSvaline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/36);  minor base G (0/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTAACGGCTACGCCAGT  >  NC_000913/4482655‑4483147
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
cAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTg                                                                                                                                                                                                                                                      <  1:422286/250‑1 (MQ=255)
      ccgGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAGACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCtt                                                                                                                                                                                                                                               <  1:46196/251‑1 (MQ=255)
        ggTTTCGCCAGGACTTAGGCACGCACGTACCACTGATCGGTCAGCATTGGTTCGATAACTACGCCTCCCCGGTCGCCGTAAGGCACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTTCAGCAAAACGATACTGTTTCTGGAACTCTGACGGGATTTCGCTTGCATAAACGTCCGATTCGTTACCTTTGGTATAGACCATCTGGGCGCTTTc                                                                                                                                                                                                                                             >  2:187385/1‑251 (MQ=255)
         gTTTCGCCAGGACATCGGCACGCACGTACCAATGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTACTTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCCACACGCTCCAGTTTCTGGAACACTGACGGGATTTCGCTGGACTAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTc                                                                                                                                                                                                                                             >  2:107027/1‑250 (MQ=255)
              gCCAGGACATCGTCACGCCCGTACCACTGGTCGGTTAGCATTGGTTTGATACATACGCCGCCACGGTCGCCTTAAGGAACGGTCAGGTCTGGCGTTTTAATTTCTTCCAGCAGTCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTTCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGat                                                                                                                                                                                                                                       <  2:398424/251‑1 (MQ=255)
                 aGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGCCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGCATAAAACGCAGATTAGTTACCTTTGGTATCGAACATCTGGGCGCATTCACggctaac                                                                                                                                                                                                                                    >  2:197970/1‑246 (MQ=255)
                        cGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTc                                                                                                                                                                                                                              >  2:186597/1‑250 (MQ=255)
                           cactcacgTACCACTGGTCGGTTAGCATTGGTTCGATCACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAg                                                                                                                                                                                                                          <  2:156050/247‑1 (MQ=255)
                                      cACTGGTCGGTCAGCATTGGTTCGATAACTACGCCTCCACGGTCGCCGTAAGGCACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGCACTCTGCAGGGATTTCGCTGGAATAAACGTAAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCACAGGTCAGGa                                                                                                                                                                                                                   >  2:363938/1‑247 (MQ=255)
                                          ggtcggtcAGCTTTGTTTCGATAACGACGCCTCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCATCTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGTGCGTTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGatc                                                                                                                                                                                                           <  2:370971/251‑1 (MQ=255)
                                          ggtcagtcctcattgtctCGAGAACTACTCCGCCACGGTCGACGTAAGGAACGGTCAGGTCGTGCGGTTTAATTGCTTCCAGCAGGCCAAGCGAGTCAACTGCGGCAACGACTGCGTTACGTGCATCAATACGCTCCAGTTTCTGGAACTCTGCATGGATTTCGCTGGAATAAACGGCAGATTCGTTACCTTTGGTATCGATCAGCTTTGCGCTTTCCCGGCCTTCGTCGTCCAAGGTCAGGATATTGatc                                                                                                                                                                                                           <  2:104535/234‑1 (MQ=255)
                                                 cAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGc                                                                                                                                                                                                     >  1:375169/1‑250 (MQ=255)
                                                 cAGCATAGGTTCGATAACTACGGCGCCACGGGCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGCCTGCTTTACGTGCAGCAAAACGCACAAGTTTCGGGAACTCTGCAGGGCTTTCGCTGGAATGAACGTAAGATTCGTTACCTTTGGAATCGAACATCTGGGCGCTTTCACGGATATAGCCGTCTACGGTCAGGATGGTGatcttccgc                                                                                                                                                                                                     >  2:257710/1‑244 (MQ=255)
                                                       tGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGt                                                                                                                                                                                              >  1:22205/1‑251 (MQ=255)
                                                        ggTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGt                                                                                                                                                                                              <  1:259687/250‑1 (MQ=255)
                                                           tcgttAACTACTCCGCCCCGGTGACCGTAAGGATCGGTCAGGTCGATCGGTTTAATGTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGTCTGCTTTACGTGCAGCAAAACGCGCCAGTTTCTGGAACTCTGCATGGATTTCTCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTTGGCGCTTTCTCGGATATCGCCGTCAAATGTCAGGATGTTGATCAGCGGCAGGGCGTGAc                                                                                                                                                                                           <  2:332221/246‑1 (MQ=255)
                                                               tAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCACAACACTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGCTATCGCCGTAAAAGGGCAGGATGTTGATCATCGGCAGGGCGAGCCGtttt                                                                                                                                                                                      >  2:534815/1‑250 (MQ=255)
                                                                   tACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTAccc                                                                                                                                                                                   >  1:85344/1‑250 (MQ=255)
                                                                     cccccccACGGTCGCCGTAAGGGACGGTCATGTCGTGCGGTTTAATTTCTTCCAGCAGGCCACGCGCGTAAACTGCGGCAACGACTGCGTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTTGGCGCTTTTACGGCTATCGCCTTCAAAGGGCAGGATGTCTATCATCGGCAGGGCGTGACGTTTAccctc                                                                                                                                                                                 <  2:510118/245‑3 (MQ=255)
                                                                      gactccaaggTCGCCGTAAAGAACGGTCAGGACGTGCGGTTTACTTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAAAGCCTGCTTTACGTGCAGCAACACGCTCCAGTTTCTGGACCTCTGCAGGGCTTTCTCTGGAATAACCTTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGTGCATGGCGCTTACCCAcc                                                                                                                                                                                >  2:452239/9‑249 (MQ=255)
                                                                           cACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGAATGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGCTCATCGGCAGGGCGTGACGTTTACCCACTTCAt                                                                                                                                                                            >  1:370298/1‑249 (MQ=255)
                                                                                   ccGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTaaa                                                                                                                                                                 <  1:465425/251‑1 (MQ=255)
                                                                                       caGGAACGGCCAGGTCGTGCGGTTTAATTTCTGCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGCCTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGCACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGCACATCTGGGCGCTTTCACGGCTAGCGCCGTCAAAGATCAGGATGTTGCTCATCGGCAGGGAGTGACGTTCACCCCATTCATAGTctgtaaagacg                                                                                                                                                              >  2:204989/2‑241 (MQ=255)
                                                                                       aaGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTTGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAggc                                                                                                                                                              <  1:138437/250‑3 (MQ=255)
                                                                                       aaGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGTATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGATCATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCg                                                                                                                                                             <  1:529042/251‑1 (MQ=255)
                                                                                       aaGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCg                                                                                                                                                             <  1:201373/251‑1 (MQ=255)
                                                                                                   ggTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTg                                                                                                                                                 <  1:142323/251‑1 (MQ=255)
                                                                                                   ggTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGGGACGTTTACCCACATCATAGGCGTTAAAGCCAGGCGCCGGAGTg                                                                                                                                                 >  2:563994/1‑251 (MQ=255)
                                                                                                           ggTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCAc                                                                                                                                         <  1:368236/251‑1 (MQ=255)
                                                                                                            gTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACg                                                                                                                                        >  1:313770/1‑251 (MQ=255)
                                                                                                                aacttcttcCAGCAGGCCAAGCGCGGCAACTGCGGCAACTACTGCTTTACGGGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTGCACGGATATCGCCGTCAAAGGGCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGGCTTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAga                                                                                                                                    <  2:312807/248‑2 (MQ=255)
                                                                                                                       tcCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGcc                                                                                                                             >  1:547165/1‑251 (MQ=255)
                                                                                                                       tcCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGAAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGt                                                                                                                                <  1:477322/248‑1 (MQ=255)
                                                                                                                       cccAGCAGGCGAGGCGCGTCAAGTGCGGCAACGACTGCTGTACGTGCGACAAAACGCTCCAGGTTCGGGAACTCTGCAGGGATTTCGCTGGAATTAACGTCAGATTCGTTACCTTTGGTAACGATCAGCTGGGCGCTTGCACCGATAACGCCGGCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCGCTTCATAGGCGGTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGcc                                                                                                                             <  2:46791/250‑1 (MQ=255)
                                                                                                                        ccAGTAGGACAAGCGCGTCCACTGCGGCCACGACTGCTTTACGTGCCGCAACACGCTCCAGTTTCTGCAACTATGCAGGGATTTCGCTGGAATAAACGTCCGATTCGTTACCTTTGGCATCGAACATCTGGGAGCTTTCACGGATATCGCAGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGGCGTTTACCCCCTTCCTAGTCGTTACAGTCGTGCGCCGGCGAGAGTTTCACGCAGCCGGGGCCt                                                                                                                            >  2:31772/1‑251 (MQ=255)
                                                                                                                        ccAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCt                                                                                                                            <  1:186597/251‑1 (MQ=255)
                                                                                                                          agcagGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTAGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGAGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCttt                                                                                                                          >  1:401840/1‑251 (MQ=255)
                                                                                                                             agGCCAAGCGTGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACGCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTc                                                                                                                       <  2:85344/251‑1 (MQ=255)
                                                                                                                              ggCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCCACCGGTGCCTTTTTcc                                                                                                                      >  1:297593/1‑251 (MQ=255)
                                                                                                                                    gcgcGTCAACTGCGGCAACGACTGCTTGACGTGCAGTAAAACGCTCCAGTTTCTTGAACTCTGCAGGTATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTGCATAGTCGTTAAAGTCGGGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTca                                                                                                                <  2:375169/251‑2 (MQ=255)
                                                                                                                                        gTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTAATAGTCGTTAAAGGCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCgt                                                                                                            >  1:337962/1‑251 (MQ=255)
                                                                                                                                         tCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCgtg                                                                                                           >  1:293329/1‑251 (MQ=255)
                                                                                                                                              tGCGGCAACGACTGCTTTACGTGAAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGCTTTTCCCGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTtcgt                                                                                                      >  2:352074/1‑251 (MQ=255)
                                                                                                                                                 ggCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTACATGTCGGCGAGTTCGTCGc                                                                                                   >  1:374755/1‑251 (MQ=255)
                                                                                                                                                     acgtcTGCTTTACGTTCAGCAAAACGCGCCAGGTTCTGGAACGCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCAGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCTTTAAATGCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAAc                                                                                               <  2:525393/247‑1 (MQ=255)
                                                                                                                                                     aCGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCAAAGTCGTTAAAGTCGTGCGCCGGAGTGCTTTTCACGCAACCGGTGCCTTTTTCCATGTCGGCGGGTTCGGCGcccac                                                                                               >  2:157651/1‑248 (MQ=255)
                                                                                                                                                         cTGCTTTACGTGCAGAAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAACCGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGCCGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGGGATTTTCACGCCGCCGGTGACTTTTTTCATGTCgggggtttcggagccaacgagc                                                                                           >  2:374197/1‑231 (MQ=255)
                                                                                                                                                          tGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGACTAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGACGTCAAAGGTCAGGATGATGATCTTCGGCAGGGCGGGACGTTTCCACACTTCGTAGTAACATAAATAGGGCGCCGGAGGGATTTTCACGACGCCGGTGCCTTTTTCCATGTCGGTGTGTTCGTCCCAAcagac                                                                                           >  2:526796/1‑245 (MQ=255)
                                                                                                                                                          tGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTTCCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCttttt                                                                                                                        <  2:529002/221‑1 (MQ=255)
                                                                                                                                                          tGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCttttt                                                                                                                        >  1:529002/1‑221 (MQ=255)
                                                                                                                                                           gCTTTACGGGCAGCAAAACGTTTACGTTGCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGGCAGATTCGTTATCATTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCTGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGagc                                                                                           <  2:309536/249‑3 (MQ=255)
                                                                                                                                                                      agcaAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGt                                                                              <  1:388548/251‑1 (MQ=255)
                                                                                                                                                                        caAAACGCTCCAGTTTCTGTAACTCTGCAGGGATTTCGCTGGAATAAGCGTCATATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATACCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCCCTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTa                                                                            <  1:187385/251‑1 (MQ=255)
                                                                                                                                                                        caAAACGCTCCAGTTTCTGGTACTCTTCAGGGAGTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGTGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTa                                                                            <  1:521112/251‑1 (MQ=255)
                                                                                                                                                                            aCGCTCCAGTTTCGGGAACTCTTCTGGGATTTCGCTTGCATAAACGTCAGATTCGTTACCTTTGGTATCTAACATCTGGGCGCTTTCACTGATATTGCCGTCAAAGGTCATGATGTTTACCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGGCGTGCGTTGGAGTGATTTTCACGCAGCCGGTGCCTCTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAAGACGACGTTTAACCa                                                                        <  2:303515/251‑1 (MQ=255)
                                                                                                                                                                                 ccAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGGGTTCGTCTCGAACGATCGGAATACGACGGTTAACCAGCGGc                                                                   <  1:452239/251‑1 (MQ=255)
                                                                                                                                                                                    gttgcttgaaatctGCAGGGATGTCCCTGGAATAAACGTCAGATGCGTTACCTTTGGTATCGATCCGCTGGGCGCTTTCACGGATAGCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTTCCCACGTCATAGTCGTTAAAGTCGTGCGCCTGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCCTGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGa                                                                <  2:102972/240‑1 (MQ=255)
                                                                                                                                                                                      gtCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTGTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGGTAGGAGGTTGATCATCGGCAGGGCGGGACGTTGACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGGTCGTCGCCAACGATCGGAATTCGACGGTTAACCAGCGGCAGaa                                                               <  2:22205/249‑1 (MQ=255)
                                                                                                                                                                                           gAACTCTGCAGGGATTTCGCTGGAATAAACGTCATAGTCGTTACCTTTGGTAGCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCa                                                                                                                                      <  2:527323/174‑1 (MQ=255)
                                                                                                                                                                                           gAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCa                                                                                                                                      >  1:527323/1‑174 (MQ=255)
                                                                                                                                                                                            aaCTCTGCAGGGATTTCGCTTGAAGAAACGTCATATTCTTTACCTGTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGGCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACata                                                        <  2:525168/251‑1 (MQ=255)
                                                                                                                                                                                                        aTTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTTACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGGTTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGGCGCCAACGATCGGGATACCACGGTTTACCAGCGGCAGAATGGCATAttatgcaatca                                             >  2:466803/1‑241 (MQ=255)
                                                                                                                                                                                                        aTTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAg                                            <  1:197970/251‑1 (MQ=255)
                                                                                                                                                                                                          ttCGCTGGAATAAACGTCAGATTCGTTACCTTTGTTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGAGCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGa                                           <  2:302723/250‑1 (MQ=255)
                                                                                                                                                                                                            cGCTGGAATAAACGGCAGATTCGTTACCTTTTGTATCGAACATCTGGGCGCTTTCACGGACATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCt                                        <  1:439197/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaaCGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAAcgc                               >  1:497272/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCGTCAGATTCGTTACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCgg                             >  1:454889/1‑251 (MQ=255)
                                                                                                                                                                                                                           tCAGATTCGTTACCTTTGGTATAGCACATCTGGGCGCTCTCACGGATATCGCCGTCATCGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTCCATAGTCGTTACAGTCGTGCGACGGAGAGCTTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCCTGTTCGTCGCAAACGGTAGGAATCCGCCGGTTAAACAGCGGAAGAATGAAATAATTGCCAATCAGCTCTTTGTAACTCgggtc                          >  2:261788/1‑248 (MQ=255)
                                                                                                                                                                                                                               aTTCGTTACCTGTGGTATCTAACATCTGGGCGCTTTCTCGGATATCGTCGTCAAAGGTCAGGATGGGGATCCTCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGCTTTTCACGCAGCCGGTGCCTTTTTCCATGGCGGCGTGTTCGTCGCCATCGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTcc                      <  2:144235/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ttACCTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGtt                 >  1:130117/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ccTTTGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGGATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGAGGATCTTCCGGGTTAAc              >  1:88414/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tttGGTATCGAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTAACgg            >  1:172733/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ggTATCGAACATCTGGGCGCTTTCACGGATAGCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTAACGGCt          <  1:374197/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  gAACATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAAt                                               >  1:423840/1‑207 (MQ=255)
                                                                                                                                                                                                                                                   aaCATCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTAACGGCTACGCCAGt  <  1:157651/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
CAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTATCGAACACCTGGGCGCTTTCACGGATATCGCCGTCAAAGGTCAGGATGTTGATCATCGGCAGGGCGTGACGTTTACCCACTTCATAGTCGTTAAAGTCGTGCGCCGGAGTGATTTTCACGCAGCCGGTGCCTTTTTCCATGTCGGCGTGTTCGTCGCCAACGATCGGAATACGACGGTTAACCAGCGGCAGAATGACATATTTGCCAATCAGATCTTTGTAACGCGGATCTTCCGGGTTAACGGCTACGCCAGT  >  NC_000913/4482655‑4483147

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