Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 6.5% intergenic (+96/‑276) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT 6.5% 512.4 / 19.0 149intergenic (+96/‑276)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (65/71);  new base T (6/4);  total (73/76)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-02

Only 100 of 149 total aligned reads displayed.

ATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGG  >  NC_000913/1425909‑1426442
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                        
 tCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACaga                                                                                                                                                                                                                                                                      <  2:311954/273‑1 (MQ=255)
         gAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGaa                                                                                                                                                                                                                                                            <  2:424427/275‑1 (MQ=255)
                        gCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCAt                                                                                                                                                                                                                                             >  2:59609/1‑275 (MQ=255)
                          ccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATcc                                                                                                                                                                                                                                           <  4:394158/275‑1 (MQ=255)
                                       cagcagAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACGGAGATAACTGACTATTTAGCCACTCATCCGGAATTGTTTGAg                                                                                                                                                                                                                              >  2:55231/1‑275 (MQ=255)
                                                      aagagaaaaaaGCCATAGCGCCCGTGNGNGGGTTTTGGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGNAGCTAGAACACCTCAGGTTGCCACAGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTNCGCCTATTCGGTAATATCCCGTGCCCGACGGTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTg                                                                                                                                                                                                                 <  3:48501/268‑1 (MQ=255)
                                                            aaaaaGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGt                                                                                                                                                                                                          <  1:45808/274‑1 (MQ=255)
                                                             aaaaGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGa                                                                                                                                                                                                        >  2:470045/1‑275 (MQ=255)
                                                                    aTAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCCATGACAGAGATAAATGAATATTTAGCCAATCATCCGGTATTGATTGAGCTCGATTAATTGGAAGCAGGGATTTTTGc                                                                                                                                                                                                 >  4:46024/1‑275 (MQ=255)
                                                                    aTAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCCTACGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTGTGCTGGACTGGCTCGGGTGCCTTTGCCTTTATCTCTGTCAGAGATAAAGTAATATTTATCCACTTGTCCGTTATTTATTGACCGCGATGAATTTGAATCAGTGATCTTTCCa                                                                                                                                                                                                >  4:7695/1‑275 (MQ=255)
                                                                             cAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTAAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGa                                                                                                                                                                                        <  2:129985/275‑1 (MQ=255)
                                                                                cGGTGGTTATCGGAGTTCGGGGGAGAAAGTGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa                                                                                                                                                                                     <  4:146515/275‑1 (MQ=255)
                                                                                     gtTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTc                                                                                                                                                                                                                                               >  4:52473/1‑212 (MQ=255)
                                                                                               ttCGGGGGAGCAAGGGGGGAAAAGGAAAGATGGAAACGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCa                                                                                                                                                                      <  1:118318/275‑1 (MQ=255)
                                                                                                tCGGGGGAGAAAGTGGGGAAAAGGCAAGATGGGAGCCAGAAAAACTCAGGTTGCCACCTGTACCGGAACCAGAAATAGGCCCGATGCGTAAGGGGTGGTTGTACGCCTGTGCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCa                                                                                                                                                                      <  4:359354/274‑1 (MQ=255)
                                                                                                       aggaaGGGGGTAAAAAGCAAGATGGAAGCAAGAAAAACTCAGGTTGCCACCCATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGGTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtgtg                                                                                                                                                               <  4:59608/271‑1 (MQ=255)
                                                                                                                       aaaagaagGAAGCTAGAAAAACTCAGGTTNCCACAGATACCGGAACCAGAAATTGNCCCGGTGGTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTTGGGTGCCTTNGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTa                                                                                                                                              <  3:91045/273‑1 (MQ=255)
                                                                                                                         aaGAGGGAAGCGGGTAAAACTCAGGTGGCCACCGAGACCGGAACCAGAAAGAGACCCGGTGCTTAAGGATGTTTTGTAGGCCTATTCGGTAATATCCCGTGCCCGACGTGATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTaca                                                                                                                                            <  4:183305/275‑1 (MQ=255)
                                                                                                                             tGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACag                                                                                                                                                                                                                                                                       >  8:209254/1‑148 (MQ=255)
                                                                                                                                      ggaaaaCTCCGGTTGCCACCGATACCGGAACCAGAAATAGACCCGTTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACag                                                                                                                                 <  4:55230/271‑1 (MQ=255)
                                                                                                                                       aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                                                                                                                                                                                                     <  6:878481/140‑1 (MQ=255)
                                                                                                                                       aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                                                                                                                                                                                                     >  5:878481/1‑140 (MQ=255)
                                                                                                                                       aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGa                                                                                                                                                                                                                                                             >  8:690991/1‑148 (MQ=255)
                                                                                                                                       aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTa                                                                                                                                                  >  1:91045/1‑255 (MQ=255)
                                                                                                                                             tCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAttt                                                                                                                                                                                                                                                       >  8:295574/1‑148 (MQ=255)
                                                                                                                                             tCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTc                                                                                                                                                                                                                                               >  3:106937/1‑156 (MQ=255)
                                                                                                                                                ggTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCAc                                                                                                                                                                                                                                                 >  2:569974/1‑151 (MQ=255)
                                                                                                                                                ggTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGgt                                                                                                                     >  4:476365/1‑275 (MQ=255)
                                                                                                                                                 gTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATcc                                                                                                                                                                                                                                           <  1:106937/156‑1 (MQ=255)
                                                                                                                                                      cACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCttt                                                                                                               >  2:393467/1‑275 (MQ=255)
                                                                                                                                                              ccGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                                                                                                                                                <  7:932020/122‑1 (MQ=255)
                                                                                                                                                                 gAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgn                                                                                                                                                                                                                                   <  7:209254/148‑2 (MQ=255)
                                                                                                                                                                 gAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGACATGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTAAGTACACTCTGTTTNNNAGATGTGATGGTGTTTTTATTTAAAtagt                                                                                                    >  3:152329/1‑272 (MQ=255)
                                                                                                                                                                       gAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACag                                                                                                                                                                                                                                                                       >  7:712618/1‑106 (MQ=255)
                                                                                                                                                                           tAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                                                                                                                                                >  6:633219/1‑109 (MQ=255)
                                                                                                                                                                           tAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                                                                                                                                                <  5:633219/109‑1 (MQ=255)
                                                                                                                                                                                 cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGGCAGAGATAAATGAATATTTAGCCACt                                                                                                                                                                                                                                                >  6:847426/1‑119 (MQ=255)
                                                                                                                                                                                 cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGGCAGAGATAAATGAATATTTAGCCACt                                                                                                                                                                                                                                                <  5:847426/119‑1 (MQ=255)
                                                                                                                                                                                 cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAtt                                                                                                                                                                                                                   >  5:202240/1‑148 (MQ=255)
                                                                                                                                                                                           ggAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTg                                                                                                                                                                                                         <  7:690991/148‑1 (MQ=255)
                                                                                                                                                                                              gttgttGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAgaga                                                                                                                              >  4:294184/1‑220 (MQ=255)
                                                                                                                                                                                                tgttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                                                                                                                                                                                                    <  8:531493/148‑1 (MQ=255)
                                                                                                                                                                                                tgttgtACGCCTATTCGGTAATATCACGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTATGTTACAGAGATGTGATGGGGTGTTTAATTAAATCGNNNATGCTAATGGAGAAAAGAATTGAGTGGc                                                                      >  3:153494/1‑274 (MQ=255)
                                                                                                                                                                                                 gttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCttt                                                                                                                                                                                                   <  7:190824/148‑1 (MQ=255)
                                                                                                                                                                                                  ttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgtg                                                                                                                          <  2:294185/220‑1 (MQ=255)
                                                                                                                                                                                                   tgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGc                                                                                                                                                                                                 >  7:970231/1‑148 (MQ=255)
                                                                                                                                                                                                     tACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAgaga                                                                                                                                                                                                                                                                    >  7:106458/1‑79 (MQ=255)
                                                                                                                                                                                                         ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTgg                                                                                                                                                                                            <  6:461482/147‑1 (MQ=255)
                                                                                                                                                                                                         ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTgg                                                                                                                                                                                            >  5:461482/1‑147 (MQ=255)
                                                                                                                                                                                                            aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                                                                                                                                                >  7:716012/1‑140 (MQ=255)
                                                                                                                                                                                                            aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                                                                                                                                                <  8:716012/140‑1 (MQ=255)
                                                                                                                                                                                                            aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgtt                                                                                                                                                          >  2:549430/1‑178 (MQ=255)
                                                                                                                                                                                                            aTTCGGTAAGATCCCGTGCCCGACGTGATGCTGGAATGGCTGGGGTGCCTTTGCCGTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGAGTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCt                                                          <  2:7695/274‑1 (MQ=255)
                                                                                                                                                                                                             ttCGGTAATATCCCGTGCCCGACGTTATGCTTGAATGGCTGGGGTGCCTTTTCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTCTTGATTGAGCGCGATGAATTTTAAGCAGGGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAACTCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCta                                                         >  2:202081/1‑274 (MQ=255)
                                                                                                                                                                                                              tCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGa                                                                                                                                                                                                                >  6:9048/1‑122 (MQ=255)
                                                                                                                                                                                                                  tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaa                                                                                                                                                                                                                                                                >  6:752257/1‑70 (MQ=255)
                                                                                                                                                                                                                  tAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaa                                                                                                                                                                                                                                                                <  5:752257/70‑1 (MQ=255)
                                                                                                                                                                                                                        cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                                                                                                                                                                                                >  6:661086/1‑96 (MQ=255)
                                                                                                                                                                                                                        cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                                                                                                                                                                                                <  5:661086/96‑1 (MQ=255)
                                                                                                                                                                                                                        cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                                                                                                                                                             >  5:311900/1‑147 (MQ=255)
                                                                                                                                                                                                                        cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                                                                                                                                                             <  6:311900/147‑1 (MQ=255)
                                                                                                                                                                                                                                cGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                                                     <  2:107350/259‑1 (MQ=255)
                                                                                                                                                                                                                                cGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTaca                                                                                                                                            <  4:56461/171‑1 (MQ=255)
                                                                                                                                                                                                                                  aCGTTATGCTGGAATGGGTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                                                                                                                                                >  8:916688/1‑118 (MQ=255)
                                                                                                                                                                                                                                  aCGTTATGCTGGAATGGGTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCa                                                                                                                                                                                                <  7:916688/118‑1 (MQ=255)
                                                                                                                                                                                                                                  aCGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTg                                                                                                                                                                                             <  7:682860/121‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtg                                                                                                                                                                 <  8:970231/148‑1 (MQ=255)
                                                                                                                                                                                                                                          cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                                                                                                                                             >  8:35577/1‑97 (MQ=255)
                                                                                                                                                                                                                                          cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                                                                                                                                             <  7:35577/97‑1 (MQ=255)
                                                                                                                                                                                                                                            ggAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                                                                                                                                               >  7:327731/1‑93 (MQ=255)
                                                                                                                                                                                                                                             gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaa                                                                                                                                                                                                                                                                 >  8:961538/1‑42 (MQ=255)
                                                                                                                                                                                                                                             gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACt                                                                                                                                                                                              <  7:75388/109‑1 (MQ=255)
                                                                                                                                                                                                                                             gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                                                                                                                                                        <  7:997947/147‑1 (MQ=255)
                                                                                                                                                                                                                                               aTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgcg                      >  2:227919/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   cTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgat                                                                                                                                                                                                                                  >  8:642718/1‑67 (MQ=255)
                                                                                                                                                                                                                                                       ggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACg                                                                                                                                                >  3:81350/1‑145 (MQ=255)
                                                                                                                                                                                                                                                        gTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATGTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTaa                                                                                                                                                     <  6:334038/139‑1 (MQ=255)
                                                                                                                                                                                                                                                           cctttgcctttATCTCTGACATAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTaca                                                                                                                                            <  1:81350/145‑1 (MQ=255)
                                                                                                                                                                                                                                                           cctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                                                                                                                                                        >  1:135215/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                 cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTaca                                                                                                                                            <  1:45933/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                                                                                                                                      <  6:751707/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                                                                                                                                      >  5:751707/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                   tttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                                                                                                                                                                                            >  6:5150/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                   tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                                                                                                                                                                                            >  5:395529/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                   tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAgg  <  2:476366/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tttATCTCTGACAGAGATAAATGAATATTTAGCCACCCATCCGGTATTGATTGAGCGCGATGaa                                                                                                                                                                                                                     <  8:700487/64‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                                                                                                                                                      >  8:832639/1‑94 (MQ=255)
                                                                                                                                                                                                                                                                    ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                                                                                                                                                        >  8:73985/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                    ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGCATTTTCAGGAGCAGTGTGTGTAgttgt                                                                                                                                                        <  7:73985/124‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                                                                                                                                                                                                    <  8:178515/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                                                                                                                                                   <  7:574226/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                                                                                                                                                   >  8:574226/1‑96 (MQ=255)
                                                                                                                                                                                                                                                                         tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTac                                                                                                                                             >  6:185225/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                          ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGc                                                                                                                                                                                                             <  7:554828/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                                                                                                                                                           <  7:16927/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                                                                                                                                                           >  8:16927/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                          ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAgg                                                                                                                                                                          <  5:900083/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAAt                                                              >  3:5088/1‑207 (MQ=255)
                                                                                                                                                                                                                                                                            gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgt                                                                                                                                                                  >  6:913687/1‑106 (MQ=255)
                                                                                                                                                                                                                                                                            gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgt                                                                                                                                                                  <  5:913687/106‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                        
ATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGG  >  NC_000913/1425909‑1426442

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