Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,130,005 G→A 5.8% L57L (TTG→TTA metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,0050GA 5.8% 564.2 / 16.6 173L57L (TTG→TTAmetLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base G (76/86);  new base A (6/4);  total (83/90)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

Only 100 of 173 total aligned reads displayed.

GGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAGTGGTGGGCCACGGGGAAGTATGGTCGGCACGTCTGATGTCTGCGGTACTTAATCAACAAGGGCTGCCAG  >  NC_000913/4129764‑4130276
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                               
       aaGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCtgtg                                                                                                                                                                                                                                         >  1:2254/1‑275 (MQ=255)
              cGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGATGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGAATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTGGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTATCTGCGCAt                                                                                                                                                                                                                                    >  2:466584/1‑273 (MQ=255)
                 agacttgaatGCCGGACTTGAAAATTGCTTCCGGGCTGCAAAAGAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGt                                                                                                                                                                                                                               <  1:84310/266‑1 (MQ=255)
                                                                            tgtgATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACt                                                                                                                                                                                                                                                                          >  2:753767/1‑173 (MQ=255)
                                                                                    cgcAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGAAGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTACGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCTGTCTGCTACCCGCTGAAGAAGCCg                                                                                                                                                              >  2:377482/1‑273 (MQ=255)
                                                                                              ggCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCCCCACTAACCAGTTGa                                                                                                                                                                                                                                                                                >  4:159001/1‑149 (MQ=255)
                                                                                               gCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGAt                                                                                                                                                                                                                                                                               <  8:647280/149‑1 (MQ=255)
                                                                                                  aaaGGTCGTCAGCTACATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTaa                                                                                                                                                                                                                                                                            >  8:581022/1‑149 (MQ=255)
                                                                                                   aaGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAAc                                                                                                                                                                                                                                                                           >  8:252090/1‑149 (MQ=255)
                                                                                                   aaGGTCGTCAGCTGCATAAATTTGGTCGCAGTAGTCTGGCTGACGTGCAGCGTTATTTGCGTGTCGCGGGCATTATGCCGGAGTACTCCCAGCCCGACGCCACGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGACTAACTGGTTGAAACt                                                                                                                                                                                                                                                                >  2:177879/1‑160 (MQ=255)
                                                                                                    aGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAAc                                                                                                                                                                                                                                                                           <  5:122718/148‑1 (MQ=255)
                                                                                                     ggTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACt                                                                                                                                                                                                                                                                          >  5:292058/1‑148 (MQ=255)
                                                                                                     ggTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTg                                                                                                                                                                                                                                                                         <  7:754146/149‑1 (MQ=255)
                                                                                                     ggTCGTCAGCTACATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTg                                                                                                                                                                                                                                                                         <  7:581022/149‑1 (MQ=255)
                                                                                                      gtcgtcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTg                                                                                                                                                                                                                                                                         <  5:481417/148‑1 (MQ=255)
                                                                                                          tcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGt                                                                                                                                                                                                                                                                       <  7:825390/146‑1 (MQ=255)
                                                                                                          tcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGt                                                                                                                                                                                                                                                                       <  7:395790/146‑1 (MQ=255)
                                                                                                          tcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGt                                                                                                                                                                                                                                                                       >  8:825390/1‑146 (MQ=255)
                                                                                                          tcAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGt                                                                                                                                                                                                                                                                       >  8:395790/1‑146 (MQ=255)
                                                                                                                gCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTaa                                                                                                                                                                                                                                                              <  8:809142/149‑1 (MQ=255)
                                                                                                                 cATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTaa                                                                                                                                                                                                                                                              >  6:498382/1‑148 (MQ=255)
                                                                                                                  aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGTCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGc                                                                                                                                                                                                                                                            <  8:347163/149‑1 (MQ=255)
                                                                                                                  aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGa                                                                                                                                                                                                                                                                                >  6:607213/1‑129 (MQ=255)
                                                                                                                  aTAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAg                                                                                                                                                                                                                                                             <  6:295362/148‑1 (MQ=255)
                                                                                                                   tAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGc                                                                                                                                                                                                                                                            >  5:331282/1‑148 (MQ=255)
                                                                                                                           gtggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGtt                                                                                                                                                                                                                                                                      >  4:746709/1‑130 (MQ=255)
                                                                                                                           gtggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGa                                                                                                                                                                                                                                                    <  6:616285/148‑1 (MQ=255)
                                                                                                                           gtggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGCAGAAGCCGATAGCCTCATCAGCGCTTTTGTCAGCGACCTTGAGCGCCt                                                                                                                      >  2:735371/1‑274 (MQ=255)
                                                                                                                            tggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                                                                                                                                                  <  8:22984/149‑1 (MQ=255)
                                                                                                                            tggCAGTAGTCTGGCTGATGTGAAGTGTTAATTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                                                                                                                                                  <  8:50883/149‑1 (MQ=255)
                                                                                                                             ggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGtta                                                                                                                                                                                                                                                                     <  8:145367/129‑2 (MQ=255)
                                                                                                                             ggCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTAATCAACTGGtta                                                                                                                                                                                                                                                                     >  7:145367/1‑128 (MQ=255)
                                                                                                                               cAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGTGNATTATGNCGTAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCCCATCAGGTTCAACAACCGCTGCGt                                                                                                                                                                                                               <  3:39686/181‑1 (MQ=255)
                                                                                                                                   agtCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGtctctc                                                                                                                                                                                                                                           <  7:845402/149‑1 (MQ=255)
                                                                                                                                    gtCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGtctctc                                                                                                                                                                                                                                           <  6:544647/148‑1 (MQ=255)
                                                                                                                                        ggCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTgcg                                                                                                                                                                                                                                       >  5:550725/1‑148 (MQ=255)
                                                                                                                                         gCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGa                                                                                                                                                                                                                                                                                >  4:355492/1‑106 (MQ=255)
                                                                                                                                          cTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCa                                                                                                                                                                                                                                     >  5:17959/1‑148 (MQ=255)
                                                                                                                                           tGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCAt                                                                                                                                                                                                                                    <  6:283147/148‑1 (MQ=255)
                                                                                                                                             aTGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTACCCAGTTGATTAACTGGTTGAAACTCAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCTTCGCTATCAGTGCGATCTGATTAGCGGTCTACTACCCGCTGAAGAAGCCCATAGCCTCATTAGCGCATTTGTCCGCGCCCTTGAGCGCCCGGCGGCGCTGCTcagcag                                                                                                    >  4:344260/1‑269 (MQ=255)
                                                                                                                                             aTGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTaa                                                                                                                                                                                                                                                                            <  2:355494/106‑1 (MQ=255)
                                                                                                                                              tgtgAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCCCCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCCGACCGATCGTCTCTCTGCGCATCAGGTTAAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGACGAAGCCGTTAGCCTCATCTGCGCTTTTGTCTGGGACCCTGAGCGCCTGGCGGCGCTtctccccagag                                                                                                  >  4:570271/1‑268 (MQ=255)
                                                                                                                                                tgAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGt                                                                                                                                                                                                                               <  5:182061/148‑1 (MQ=255)
                                                                                                                                                   aGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaa                                                                                                                                                                                                                           <  7:226872/149‑1 (MQ=255)
                                                                                                                                                   aGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCCTTAGCGCTTTTGTCAGCGACCTTGAGCCCCTGGCCGCCCTGCCCGACAGCGGTAtt                                                                                             >  4:660240/1‑275 (MQ=255)
                                                                                                                                                   aGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCCGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTAtt                                                                                             >  2:593417/1‑275 (MQ=255)
                                                                                                                                                      gtTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaaca                                                                                                                                                                                                                         >  6:494658/1‑148 (MQ=255)
                                                                                                                                                       ttATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaacaa                                                                                                                                                                                                                        <  6:17959/148‑1 (MQ=255)
                                                                                                                                                        tATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAAc                                                                                                                                                                                                                      <  7:252182/149‑1 (MQ=255)
                                                                                                                                                          tttGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGc                                                                                                                                                                                                                    >  8:790015/1‑149 (MQ=255)
                                                                                                                                                            tGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGc                                                                                     <  2:660242/274‑1 (MQ=255)
                                                                                                                                                            tGCGTGTCGCGGGCATTATGGCGGAGTAATCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCa                                                                                    <  2:344262/275‑1 (MQ=255)
                                                                                                                                                               gtgtCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATAGTTTCCCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGAGCGATCTGATTAGCGGTCTGCTACCTGCTGACGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCGTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAggt                                                                                 >  2:264340/1‑273 (MQ=255)
                                                                                                                                                                tgtCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGaa                                                                                                                                                                                                                                                                   <  5:865224/96‑1 (MQ=255)
                                                                                                                                                                  tCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGc                                                                                                                                                                                                            >  7:831068/1‑149 (MQ=255)
                                                                                                                                                                       gtcATTTTGGGGGATAATTCGCAGCCTTANGTTATGTTGGTGGTTTCGGCCGCCGNTAGACACACTAAGAAGTTGATTACGTGGTTGAAACTAGGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCNGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGATGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGa                                                                         <  3:142128/273‑1 (MQ=255)
                                                                                                                                                                       ggCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCa                                                                                                                                                                                                       <  7:287787/149‑1 (MQ=255)
                                                                                                                                                                            tATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCAt                                                                                                                                                                                                                                    >  4:184848/1‑115 (MQ=255)
                                                                                                                                                                               ggCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                                                                                                                                                  >  8:182065/1‑98 (MQ=255)
                                                                                                                                                                               ggCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATc                                                                                                                                                                                                                                                  <  7:182065/98‑1 (MQ=255)
                                                                                                                                                                               ggCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGAt                                                                                                                                                                                                >  6:841849/1‑148 (MQ=255)
                                                                                                                                                                               ggCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGAt                                                                                                                                                                                                <  5:421273/148‑1 (MQ=255)
                                                                                                                                                                               ggCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATc                                                                                                                                                                                               >  8:357023/1‑149 (MQ=255)
                                                                                                                                                                               ggCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATc                                                                                                                                                                                               <  8:752805/149‑1 (MQ=255)
                                                                                                                                                                                gCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAgg                                                                                                                                                                                                                                <  2:184849/115‑1 (MQ=255)
                                                                                                                                                                                  ggAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaagaa                                                                                                                                                                  >  1:175323/1‑175 (MQ=255)
                                                                                                                                                                                    aGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACt                                                                                                                                                                                                                                                                          >  2:200424/1‑69 (MQ=255)
                                                                                                                                                                                      tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTa                                                                                                                                                                                         <  5:384302/148‑1 (MQ=255)
                                                                                                                                                                                      tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCg                                                                                                                                                              <  3:175322/175‑1 (MQ=255)
                                                                                                                                                                                      tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAGGGGTGGGCCACgg                                                          >  2:114366/1‑275 (MQ=255)
                                                                                                                                                                                        ctctCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGtt                                                                                                                                                                                                                                                                      <  4:200423/69‑1 (MQ=255)
                                                                                                                                                                                        ctctCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTAAAACTAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTTGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCCGTGTATGCGGAAGTGGTGGGCCAAGGGGa                                                       >  2:158284/1‑275 (MQ=255)
                                                                                                                                                                                          ctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCgg                                                                                                                                                                                     >  6:735886/1‑148 (MQ=255)
                                                                                                                                                                                            cAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATc                                                                                                                                                                                                        <  6:227535/127‑1 (MQ=255)
                                                                                                                                                                                               ccTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCa                                                                                    <  4:31862/240‑1 (MQ=255)
                                                                                                                                                                                                   aCGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGAAGTGTATGCGGAAGTGGTGGGCCACGGGGAAGTATGGt                                               >  2:674560/1‑273 (MQ=255)
                                                                                                                                                                                                      atatGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCACTATCAGTGCGATCTGATTAGc                                                                                                                                                                                       >  8:175795/1‑134 (MQ=255)
                                                                                                                                                                                                         tgatgGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTga                                                                                                                                                                      <  5:621918/148‑1 (MQ=255)
                                                                                                                                                                                                            tggtggTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaaga                                                                                                                                                                   <  5:262749/148‑1 (MQ=255)
                                                                                                                                                                                                              gtggtTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGaagaag                                                                                                                                                                 <  5:692276/148‑1 (MQ=255)
                                                                                                                                                                                                              gtggtTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGaagaag                                                                                                                                                                 <  5:627190/148‑1 (MQ=255)
                                                                                                                                                                                                              gtggtTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGc                                                                                                                                                                >  7:354365/1‑149 (MQ=255)
                                                                                                                                                                                                                  tttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGAt                                                                                                                                                                                                <  7:807201/113‑1 (MQ=255)
                                                                                                                                                                                                                   ttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTg                                                                                                                                                                                             >  6:458835/1‑115 (MQ=255)
                                                                                                                                                                                                                    tCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAg                                                                                                                                                          <  8:296918/149‑1 (MQ=255)
                                                                                                                                                                                                                      cgccgccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGc                                                                                                                                                         <  6:624653/148‑1 (MQ=255)
                                                                                                                                                                                                                          gccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAGGCCAGACCGATCGTCTCTCTGCGCATCAGGTGCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCgg                                                                          <  1:171374/223‑1 (MQ=255)
                                                                                                                                                                                                                            cgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTa                                                                                                                                                  <  8:275306/149‑1 (MQ=255)
                                                                                                                                                                                                                                 gCACAACTAACCAGTGGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGGGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAGTGGTGGGCCACGGGGAAGTATGGTCGGCACGTCTGATGTCTGCGGTACTTAATcaa               <  2:570273/275‑1 (MQ=255)
                                                                                                                                                                                                                                  caccacTAACCAGTTGATTAACTGGTTGCAACTAAGCCAGACCGAACGTCTCTCTGCGCATCAGGTTCAACAAACGCt                                                                                                                                                                                                                   >  4:815244/1‑78 (MQ=255)
                                                                                                                                                                                                                                    ccacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtttt                                                                                                                                          <  8:487236/149‑1 (MQ=255)
                                                                                                                                                                                                                                     cacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtt                                                                                                                                            <  7:5089/146‑1 (MQ=255)
                                                                                                                                                                                                                                     cacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTg                                                                                                                                         >  8:600112/1‑149 (MQ=255)
                                                                                                                                                                                                                                          aCCAGTTAATCAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGc                                                                                                                                    >  8:687234/1‑149 (MQ=255)
                                                                                                                                                                                                                                               ttGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCt                                                                                                                                                                                                                   >  1:51709/1‑65 (MQ=255)
                                                                                                                                                                                                                                               ttGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCt                                                                                                                               >  7:773544/1‑149 (MQ=255)
                                                                                                                                                                                                                                               ttGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGCCAGCGGTATTAACGACGCAGTGTATGCGGAAGTGGTGGGCCACGGGGAAGTATGGTCGGCACGTCTGATGTCTGCGGTACTTAATCAACAAGGGCTGCCAg  >  2:436323/1‑274 (MQ=255)
                                                                                                                                                                                                                                               ttGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCCCTCTGCGCCTCAGGTTCAACAAACGCTGCCTCGCTATCAGTGCGAACTGATTAGCGCTCTGCTACCCGCTGAAGAAGCCGATAGCCGCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGACTGCGCTGCTCGACAGCGGTATTAACGACGCAGTACATGCGGCAGTGGTGGGCCACGGGGGAGTATGGCACGCACGTCCCCAGTCAGTGTTACTTACTCAACAAggtcctgcag  >  4:401084/1‑266 (MQ=255)
                                                                                                                                                                                                                                                tGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCg                                                                                                                                                              <  5:6832/117‑1 (MQ=255)
                                                                                                                                                                                                                                                tGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCg                                                                                                                                                              >  6:6832/1‑117 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                               
GGTATTGAAGATGGCGAAGATTTAATTGCCGACCTGGAAAATGGCTTCCGGGCTGCAAACAAGGGGTAAAAATGAGTGTGATTGCGCAGGCAGGGGCGAAAGGTCGTCAGCTGCATAAATTTGGTGGCAGTAGTCTGGCTGATGTGAAGTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAGTGGTGGGCCACGGGGAAGTATGGTCGGCACGTCTGATGTCTGCGGTACTTAATCAACAAGGGCTGCCAG  >  NC_000913/4129764‑4130276

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