Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,328,184 A→T 5.8% intergenic (+101/+55) yhbY → / ← greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,328,1840AT 5.8% 748.2 / 24.5 208intergenic (+101/+55)yhbY/greARNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Reads supporting (aligned to +/- strand):  ref base A (100/96);  new base T (8/4);  total (108/100)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

Only 100 of 211 total aligned reads displayed.

GAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAGATCGCGAAACTAAAACCTTGATCGTGGAAGCTATCGTGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCACGAGCAATAGGCGAGTTTACAGAAATCAGGTTTTGTTTAAAGTCAGCTTCGTCATCGCCAACGATGCGATAAGTCTGTTCTTCGTCAGAATCCAGATTCAGCACCGTTACGGTAGCACCAAAA  >  NC_000913/3327913‑3328455
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                               
gagCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAGATCGCGAAACTAAAACCTTGATCGTGGAAGCTATCGTGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAg                                                                                                                                                                                                                                                                                <  4:250674/273‑1 (MQ=255)
               atcaAGGTGAAAATCGCCACCGAAGATCGCGAAACTAAAACCTTGATCGTGGAAGCTATCGTGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttct                                                                                                                                                                                                                                                               >  2:487641/1‑275 (MQ=255)
                                    gAAGATCGCAAAAGTAAAGCCTTTATCGTGGAAGCTAACGTGCGAAAAACCGGCGCCTGTAAGGTACAGGTCATCGGTAAAACGCTGGGGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATGATCCTAAAGTTACACACATTTGCTGTGTAAAAAGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCtt                                                                                                                                                                                                                                           <  1:189566/274‑1 (MQ=255)
                                               gagtAAAACACTGATCCTGGAAGCGATGCTGAGCGAAGCCGGCGCCAGTAATGTACAGGTTAGAGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAGAATGTCGCTGCAACGTTAAGATGATCCTAAAGTTGCACACATTTGCTGTGTAAAACGATGGGTTTTCCGCAGGCAAGAGGGAAAAGGGCCACGCTCTGTTCGTTGATAAAAGGCCACATAGCGGCCGTTTTGCTTTCTTTACGATACATCAACATCTGGAGTATTGGGTa                                                                                                                                                                                                                               <  4:534129/272‑1 (MQ=255)
                                                        cTTGATCGTGGAAGCTATCGTGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACa                                                                                                                                                                                                                      <  4:487641/275‑1 (MQ=255)
                                                                      cTATCGTGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCt                                                                                                                                                                                                                                                                          >  4:278838/1‑209 (MQ=255)
                                                                          cGTGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCttttt                                                                                                                                                                                                                                                                      <  2:278838/209‑1 (MQ=255)
                                                                            tGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCtttt                                                                                                                                                                                                                                                                       >  2:678823/1‑206 (MQ=255)
                                                                                                             tCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAGTa                                                                                                                                                                                                                                                       <  2:689189/189‑1 (MQ=255)
                                                                                                                gTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                                                                                                                                                       <  2:735724/186‑1 (MQ=255)
                                                                                                                  aaaaCGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtg                                                                                                                                                                    <  4:301524/267‑1 (MQ=255)
                                                                                                                   aaaCGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTa                                                                                                                                                                                                    >  2:635811/1‑234 (MQ=255)
                                                                                                                       gCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTc                                                                                                                                                                                                <  4:635811/234‑1 (MQ=255)
                                                                                                                        cTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCacaaca                                                                                                                                                        <  2:246843/273‑1 (MQ=255)
                                                                                                                           gTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTg                                                                                                                                                                                                                                          >  4:684335/1‑188 (MQ=255)
                                                                                                                            tGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAg                                                                                                                                                                                                                                                                                <  6:247352/149‑1 (MQ=255)
                                                                                                                               tttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATCAAAGGCCGCATAGCgg                                                                                                                                                                                                                                                                             <  8:738128/149‑1 (MQ=255)
                                                                                                                               tttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCgg                                                                                                                                                                                                                                                                             >  5:728277/1‑149 (MQ=255)
                                                                                                                                 tATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGcc                                                                                                                                                                                                                                                                           >  7:403179/1‑149 (MQ=255)
                                                                                                                                   tCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTTGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCTCCGCCCGGCGTTTTGATGACCCAACCATCAtcttcttct                                                                                                                                            >  1:93441/1‑274 (MQ=255)
                                                                                                                                     gCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttcttt                                                                                                                                                                                                                                                             >  1:115380/1‑159 (MQ=255)
                                                                                                                                      cccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGGTTTGCTGACCACAACATCATCATCTTCTTTGc                                                                                                                                        >  2:327884/1‑275 (MQ=255)
                                                                                                                                         aaCTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACaa                                                                                                                                                                                                                                                         <  3:115380/159‑1 (MQ=255)
                                                                                                                                             aaaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttct                                                                                                                                                                                                                                                               <  6:937527/149‑1 (MQ=255)
                                                                                                                                             aaaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttct                                                                                                                                                                                                                                                               <  6:330760/149‑1 (MQ=255)
                                                                                                                                              aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttctt                                                                                                                                                                                                                                                              >  6:432334/1‑149 (MQ=255)
                                                                                                                                              aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttctt                                                                                                                                                                                                                                                              >  6:180057/1‑149 (MQ=255)
                                                                                                                                              aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCcttt                                                                                                                                                                                                                                                                 <  8:381466/146‑1 (MQ=255)
                                                                                                                                                gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttc                                                                                                                                                                                                                                                                >  5:845653/1‑145 (MQ=255)
                                                                                                                                                gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttc                                                                                                                                                                                                                                                                <  6:845653/145‑1 (MQ=255)
                                                                                                                                                gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTa                                                                                                                                                                                                                                                            >  7:839255/1‑149 (MQ=255)
                                                                                                                                                gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTa                                                                                                                                                                                                                                                            <  5:948700/149‑1 (MQ=255)
                                                                                                                                                    gTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAAt                                                                                                                                                                                                                                                        >  5:131180/1‑149 (MQ=255)
                                                                                                                                                     tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                                                                                                                                                       >  6:969541/1‑149 (MQ=255)
                                                                                                                                                     tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                                                                                                                                                       <  7:815399/149‑1 (MQ=255)
                                                                                                                                                      aaaaTCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCTCCCGGCGTTTTGATGACCACAACATCATTTCTTTCTTTGCATGTCGGGCCACggac                                                                                                                        >  4:397351/1‑272 (MQ=255)
                                                                                                                                                         aTCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATc                                                                                                                                                                                                                                                   >  7:530868/1‑149 (MQ=255)
                                                                                                                                                                    aCGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                                                                                                                                                                        <  5:432334/149‑1 (MQ=255)
                                                                                                                                                                        tAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                                                                                                                                                                    >  7:48266/1‑149 (MQ=255)
                                                                                                                                                                           gATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTggg                                                                                                                                                                                                                                 <  7:356716/149‑1 (MQ=255)
                                                                                                                                                                           gATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTggg                                                                                                                                                                                                                                 >  6:466067/1‑149 (MQ=255)
                                                                                                                                                                              tATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTaa                                                                                                                                                                                                                              <  5:475242/149‑1 (MQ=255)
                                                                                                                                                                                tCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAAtt                                                                                                                                                                                                                            <  6:504841/149‑1 (MQ=255)
                                                                                                                                                                                tCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCACGAGCAATAGGCGAGTTTACAGAAATCAGGt                                                                                              >  2:204665/1‑275 (MQ=255)
                                                                                                                                                                                tCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCACGAGCAATAGGCGAGTTTACAGAAATCAGGt                                                                                              >  2:416639/1‑275 (MQ=255)
                                                                                                                                                                                  cTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                                                                                                                                                          >  8:329836/1‑149 (MQ=255)
                                                                                                                                                                                  cTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                                                                                                                                                          <  7:344279/149‑1 (MQ=255)
                                                                                                                                                                                    aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACaa                                                                                                                                                                                                                                                         <  4:142989/116‑1 (MQ=255)
                                                                                                                                                                                    aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATc                                                                                                                                                                                                                                                   <  4:164009/122‑1 (MQ=255)
                                                                                                                                                                                    aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                                                                                                                                                          <  4:204635/147‑1 (MQ=255)
                                                                                                                                                                                    aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTa                                                                                                                                                                                                                        <  7:86965/149‑1 (MQ=255)
                                                                                                                                                                                    aaaGTCACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAg                                                                                                                                                                                                                                                                                <  1:19885/93‑1 (MQ=255)
                                                                                                                                                                                       gTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAg                                                                                                                                                                                                                     <  6:450394/149‑1 (MQ=255)
                                                                                                                                                                                          acacacATTAGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGGTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTa                                                                                                                                                                                                                  <  8:568496/149‑1 (MQ=255)
                                                                                                                                                                                            acacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                                                                                                                                                                >  5:1079770/1‑149 (MQ=255)
                                                                                                                                                                                              acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATc                                                                                                                                                                                                                                             <  7:315156/118‑1 (MQ=255)
                                                                                                                                                                                              acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATc                                                                                                                                                                                                                                             >  8:315156/1‑118 (MQ=255)
                                                                                                                                                                                                aTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAc                                                                                                                                                                                                            >  7:796384/1‑149 (MQ=255)
                                                                                                                                                                                                 tttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAcc                                                                                                                                                                                                           >  7:1080470/1‑149 (MQ=255)
                                                                                                                                                                                                 tttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAcc                                                                                                                                                                                                           >  6:660061/1‑149 (MQ=255)
                                                                                                                                                                                                      tgtgTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAt                                                                                                                                                                                                      >  5:530313/1‑149 (MQ=255)
                                                                                                                                                                                                        tgcaaaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAGTa                                                                                                                                                                                                    >  7:718551/4‑149 (MQ=255)
                                                                                                                                                                                                          tAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACt                                                                                                                                                                                                  <  8:137761/149‑1 (MQ=255)
                                                                                                                                                                                                           aaaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACtt                                                                                                                                                                                                 >  6:7633/1‑149 (MQ=255)
                                                                                                                                                                                                            aaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCg                                                                                                                                                                                                                                                                              <  2:303981/71‑1 (MQ=255)
                                                                                                                                                                                                            aaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCtt                                                                                                                                                                                                                                           <  2:191806/106‑1 (MQ=255)
                                                                                                                                                                                                                gAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAAt                                                                                                                                                                                            <  5:973956/149‑1 (MQ=255)
                                                                                                                                                                                                                 gggtgTTTTCCGCAGGCAGGAGATACAAANCCCACTTTCTGTTCGTTTATAAACGNCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTANTTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCACGAGCAATAGGCGAGTTTACAGAAATCAGGTTTTGTTTAAAGTCAGCTTCGTCATCGCCAACGa                                                             <  3:132152/271‑1 (MQ=255)
                                                                                                                                                                                                                  ggggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCACGAGCAATAGGCGAGTTTACAGAAATCAGGTTTTGTTTAAAGTCAGCTTCGTCATCGCCAACGAt                                                            >  2:21934/1‑275 (MQ=255)
                                                                                                                                                                                                                    ggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTa                                                                                                                                                                                        >  5:280775/1‑149 (MQ=255)
                                                                                                                                                                                                                      ttttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACt                                                                                                                                                                                      >  5:347623/1‑149 (MQ=255)
                                                                                                                                                                                                                      ttttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGAc                                                                                                                                                               <  2:540097/172‑1 (MQ=255)
                                                                                                                                                                                                                           cGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTcgcc                                                                                                                                                                                 >  5:183507/1‑149 (MQ=255)
                                                                                                                                                                                                                           cGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTcgcc                                                                                                                                                                                 <  7:188325/149‑1 (MQ=255)
                                                                                                                                                                                                                                gcaggCGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAgg                                                                                                                                                                                                                    >  4:71274/1‑109 (MQ=255)
                                                                                                                                                                                                                                    gagaGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGt                                                                                                                                                                                                                                <  6:979312/93‑1 (MQ=255)
                                                                                                                                                                                                                                      gagCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTCCGCCGCCCGGCGttt                                                                                                                                                                      <  8:796384/149‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgat                                                                                                                                                                  <  7:293022/149‑1 (MQ=255)
                                                                                                                                                                                                                                          aaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgat                                                                                                                                                                  <  5:671722/149‑1 (MQ=255)
                                                                                                                                                                                                                                            aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatga                                                                                                                                                                <  7:103827/149‑1 (MQ=255)
                                                                                                                                                                                                                                              tGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGAcc                                                                                                                                                              <  2:143292/149‑1 (MQ=255)
                                                                                                                                                                                                                                               gCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTgg                                                                                                                                                                                                                                  >  5:8661/1‑80 (MQ=255)
                                                                                                                                                                                                                                               gCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTgg                                                                                                                                                                                                                                  <  6:8661/80‑1 (MQ=255)
                                                                                                                                                                                                                                                ccACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcac                                                                                                                                                            <  8:875455/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 cACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCaaa                                                                                                                                                                                             <  6:1080350/115‑1 (MQ=255)
                                                                                                                                                                                                                                                      tctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatc                                                                                                                                                      >  6:557829/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            cGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttc                                                                                                                                                >  6:626318/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            cGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttc                                                                                                                                                <  5:969541/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             gTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttct                                                                                                                                               >  8:426350/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aggaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttc                                                                                                                                             >  8:809473/5‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt                                                                                                                                          >  6:315104/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGtttt                                                                                                                                                                     <  5:959179/121‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt                                                                                                                                          <  7:326829/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTAc                                                                                                                                                                                                   >  5:464432/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                   aaaaGGCCGCATAGCGGCCTTTTCCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCACCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg                                                                                                                                         >  8:105204/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGc                                                                                                                                        >  7:599607/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGc                                                                                                                                        >  7:979530/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCACGAGCAATAGGCGAGTTTACAGAAATCAGGTTTTGTTTAAAGTCAGCTTCGTCATCGCCAACGAt                                                            >  4:96526/1‑225 (MQ=255)
                                                                                                                                                                                                                                                                        gCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGAt                                                                                                                                    <  8:48266/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcat                                                                                                                                                       >  2:317718/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                               
GAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAGATCGCGAAACTAAAACCTTGATCGTGGAAGCTATCGTGCGCGAAACCGGCGCCTGTAATGTACAGGTCATCGGTAAAACGCTGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCACGAGCAATAGGCGAGTTTACAGAAATCAGGTTTTGTTTAAAGTCAGCTTCGTCATCGCCAACGATGCGATAAGTCTGTTCTTCGTCAGAATCCAGATTCAGCACCGTTACGGTAGCACCAAAA  >  NC_000913/3327913‑3328455

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