Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,235,905 16376 bp→AAGTAAAGGAAAACATCAACGTTCGCAAAAGAAAATGCCGATTTGCGGCCGC yjbEmalM yjbE, yjbF, yjbG, yjbH, yjbT, psiE, xylE, malH, malG, malF, malE, malK, lamB, malM

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4235905 4252280 16376 83 [0] [0] 83 yjbE–malM yjbE,yjbF,yjbG,yjbH,yjbT,psiE,xylE,malH,malG,malF,malE,malK,lamB,malM

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 42359040 (0.000)71 (0.950)
+AAGTAAAGGAAAACATCAACGTTCGCAAAAGAAAATGCCGATTTGCGGCCGC
62/384 0.1 100% intergenic (+497/‑2) pgi/yjbE glucose‑6‑phosphate isomerase/uncharacterized protein YjbE
?NC_000913 4252281 = 0 (0.000)intergenic (+1807/‑225) malM/ubiC maltose regulon periplasmic protein/chorismate lyase

TTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGT  >  NC_000913/4252032‑4252509
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
ttgtttttattttaagtttatgatttttattgttatttaaatataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgct                                                                                                                                                                                                                                      >  1:202963‑M2/250‑250 (MQ=255)
                tttatgatttttattgttatttaaatataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATc                                                                                                                                                                                                                     >  1:116669‑M2/234‑251 (MQ=255)
                tttatgatttttattgttatttaaatataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATc                                                                                                                                                                                                                     >  1:404551‑M2/234‑251 (MQ=255)
                    tgatttttattgttatttaaatataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCg                                                                                                                                                                                                                 >  1:640516‑M2/230‑251 (MQ=255)
                                     taaatataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGAc                                                                                                                                                                                                >  2:372038‑M2/213‑251 (MQ=255)
                                         tataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTg                                                                                                                                                                                            >  2:773396‑M2/209‑251 (MQ=255)
                                          ataagttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGa                                                                                                                                                                                           <  2:472875‑M2/44‑1 (MQ=255)
                                              gttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGcc                                                                                                                                                                                       >  1:812965‑M2/204‑251 (MQ=255)
                                               ttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggcaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGccc                                                                                                                                                                                      >  1:61067‑M2/203‑251 (MQ=255)
                                               ttgaaacttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGccc                                                                                                                                                                                      >  2:131997‑M2/203‑251 (MQ=255)
                                                     cttatatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACgg                                                                                                                                                                                >  1:609483‑M2/197‑251 (MQ=255)
                                                     cttatatttcatattcattccacttatcctaaaacgccatcgctaattccccgcgccgtacttcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaacttgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTACAGGAAAACATCACCGTCCCTGGCATCCTGGACGGTGATGCCCCTACgg                                                                                                                                                                                >  2:687749‑M2/197‑251 (MQ=255)
                                                        atatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTg                                                                                                                                                                             >  1:306279‑M2/194‑251 (MQ=255)
                                                        atatttgatattcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTg                                                                                                                                                                             <  2:194737‑M2/58‑1 (MQ=255)
                                                                 attcattccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGcc                                                                                                                                                                     <  1:522575‑M2/66‑1 (MQ=255)
                                                                  ttcattccaattatcctaaaacgccatcgctaattccccgcgccgtacttcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGcc                                                                                                                                                                     >  2:301835‑M2/184‑249 (MQ=255)
                                                                    cattccaattatcctaaaacgccatcgctaattacccgcgccgtaattcgcatgctttagttttgtatactcgatcccgcccgaaatctttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCGGGACGGAGATGCCCATACGGTTGCCCTCGCCAgcc                                                                                                                                                                 >  1:49156‑M2/182‑250 (MQ=255)
                                                                      ttccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACg                                                                                                                                                               >  2:440212‑M2/180‑251 (MQ=255)
                                                                        ccaattatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACggg                                                                                                                                                             >  1:195127‑M2/178‑251 (MQ=255)
                                                                           attatcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCAt                                                                                                                                                          >  1:698519‑M2/175‑251 (MQ=255)
                                                                            ttatcctaaaaagccatcgctacttcaccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaacttgtaagtacaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggcctcTAAGAAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGAGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATc                                                                                                                                                         >  2:557666‑M2/174‑251 (MQ=255)
                                                                             tatcctaaaacgacatagctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCg                                                                                                                                                        >  2:734128‑M2/173‑251 (MQ=255)
                                                                              atcctaaaacgccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCgg                                                                                                                                                       >  1:282028‑M2/172‑251 (MQ=255)
                                                                                        gccatcgctaattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTa                                                                                                                                              >  2:172914‑M2/162‑250 (MQ=255)
                                                                                                 aattccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTGGCCCTCGCCAGCACGGGCATCGGTAAAGGGTAAGGTTCaa                                                                                                                                      <  2:265220‑M2/97‑1 (MQ=255)
                                                                                                    tccccgcgccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATc                                                                                                                                  <  1:609798‑M2/101‑1 (MQ=255)
                                                                                                          cgccgtaattcgcatgcgttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGTTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTAc                                                                                                                           <  2:202963‑M2/108‑1 (MQ=255)
                                                                                                           gccgtaattcgcatgctttagttgtttatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACAtcgtttttcc                                                                                                                          <  2:770638‑M2/109‑4 (MQ=255)
                                                                                                           gccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTAcc                                                                                                                          <  2:849745‑M2/109‑1 (MQ=255)
                                                                                                           gccgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTAcc                                                                                                                          >  1:220487‑M2/143‑251 (MQ=255)
                                                                                                             cgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCAc                                                                                                                        >  2:567236‑M2/141‑251 (MQ=255)
                                                                                                             cgtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCAc                                                                                                                        >  1:474895‑M2/141‑251 (MQ=255)
                                                                                                              gtaattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACt                                                                                                                       >  2:283092‑M2/140‑251 (MQ=255)
                                                                                                               taattcgcatgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACtt                                                                                                                      >  2:528413‑M2/139‑251 (MQ=255)
                                                                                                                aattcgcatgcttgagttgtgtatactcgatcccgcccgaaaggtttttgggtaaaactccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgagttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGAGATGCCCCTACGTGTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTGCAACATCGTTTTTCCACtt                                                                                                                      <  2:298647‑M2/113‑1 (MQ=255)
                                                                                                                    cgcatgctttagtggtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGc                                                                                                                 <  1:523778‑M2/118‑1 (MQ=255)
                                                                                                                       atgctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAt                                                                                                              >  1:553578‑M2/131‑251 (MQ=255)
                                                                                                                         gctttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAtt                                                                                                             <  1:172914‑M2/122‑1 (MQ=255)
                                                                                                                            ttagttgtgtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAttgtt                                                                                                          <  2:553578‑M2/125‑1 (MQ=255)
                                                                                                                                    gtatactcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGttttt                                                                                                  <  2:346805‑M2/133‑1 (MQ=255)
                                                                                                                                        actcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGc                                                                                             >  1:605132‑M2/114‑251 (MQ=255)
                                                                                                                                         ctcgatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCg                                                                                            >  2:156516‑M2/113‑251 (MQ=255)
                                                                                                                                            gatcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTAt                                                                                         >  2:443762‑M2/110‑251 (MQ=255)
                                                                                                                                             atcccgcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTAtt                                                                                        <  2:358746‑M2/143‑1 (MQ=255)
                                                                                                                                                  gcccgaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTAtt                                                                                        >  1:562914‑M2/104‑246 (MQ=255)
                                                                                                                                                    ccgaaatgtttttgggtaaatctccattccttaaatgaagggaaattgttagaaaaggcaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCAc                                                                                 >  1:674969‑M2/102‑251 (MQ=255)
                                                                                                                                                      gaaatgtttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACtt                                                                               >  2:834082‑M2/100‑251 (MQ=255)
                                                                                                                                                       aaagggttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgagttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGGCGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATGGTTGCGTTTTTGTTGCGTATTAGATCACTTa                                                                              <  2:61067‑M2/153‑1 (MQ=255)
                                                                                                                                                            tttttgggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTg                                                                         >  2:670198‑M2/94‑251 (MQ=255)
                                                                                                                                                                 gggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTa                                                                    <  2:1863‑M2/163‑1 (MQ=255)
                                                                                                                                                                  ggtaaatctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTAc                                                                   <  1:654322‑M2/164‑1 (MQ=255)
                                                                                                                                                                     aaatctccattcattcaatgaagggaaattgtaagtaaatgaaaacatcaacgttcgcaaaagaaaatgccgatttgcggcctcTAAGTAAAGGAAAACATCACCGTGCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTGTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACAtc                                                                <  2:135883‑M2/167‑1 (MQ=255)
                                                                                                                                                                       atctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACAtctc                                                              >  2:80622‑M2/83‑251 (MQ=255)
                                                                                                                                                                       atctccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACAtctc                                                              >  1:816374‑M2/83‑251 (MQ=255)
                                                                                                                                                                           ccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACAtctc                                                              <  1:80622‑M2/169‑1 (MQ=255)
                                                                                                                                                                           ccattcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGt                                                          <  1:734128‑M2/173‑1 (MQ=255)
                                                                                                                                                                             attcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTa                                                         >  1:533463‑M2/77‑250 (MQ=255)
                                                                                                                                                                               tcattcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAac                                                      >  1:244648‑M2/75‑251 (MQ=255)
                                                                                                                                                                                 attcaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAacac                                                    >  1:161971‑M2/73‑251 (MQ=255)
                                                                                                                                                                                    caatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACttt                                                 <  1:184214‑M2/182‑1 (MQ=255)
                                                                                                                                                                                    caatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACttt                                                 <  1:283092‑M2/182‑1 (MQ=255)
                                                                                                                                                                                    caatgaagggaaattgtaagtaaaggaaaaaatcaacgttcgcaaaagaaaatgacgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTAAGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACttt                                                 >  1:634391‑M2/70‑251 (MQ=255)
                                                                                                                                                                                    aaatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAACGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTACGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACttt                                                 >  2:526220‑M2/70‑251 (MQ=255)
                                                                                                                                                                                     aatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACtttt                                                >  1:372335‑M2/69‑251 (MQ=255)
                                                                                                                                                                                     aatgaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACtttt                                                >  2:537405‑M2/69‑251 (MQ=255)
                                                                                                                                                                                        gaagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCt                                              <  1:90179‑M2/185‑1 (MQ=255)
                                                                                                                                                                                         aagggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAATGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGc                                            >  2:839272‑M2/65‑251 (MQ=255)
                                                                                                                                                                                          agggaaattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCg                                           >  2:668267‑M2/64‑251 (MQ=255)
                                                                                                                                                                                               aattgtaagtaaaggaaaacatcaccgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACa                                      <  2:372335‑M2/193‑1 (MQ=255)
                                                                                                                                                                                               aattgtaagtaaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATAc                                       >  2:152887‑M2/59‑250 (MQ=255)
                                                                                                                                                                                                   gtaagtaaaggcaaacatcaacgttcgaaaaagacaatgccgatttgcggccgaTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGATGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATAGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGc                                  >  2:175849‑M2/55‑251 (MQ=255)
                                                                                                                                                                                                         aaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAt                                                                                                              >  2:790549‑M2/49‑169 (MQ=255)
                                                                                                                                                                                                         aaaggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGAt                                                                                                              <  1:790549‑M2/121‑1 (MQ=255)
                                                                                                                                                                                                           aggaaaacatcaaggttcgcaaaagaaaatgccgctttgcggccgcTAAGTAAAGGAAAACATCACAGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTAACACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAAtgcctttcagt                          >  2:558667‑M2/47‑247 (MQ=255)
                                                                                                                                                                                                            ggaaaacatcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGtt                         >  1:107669‑M2/46‑251 (MQ=255)
                                                                                                                                                                                                                   atcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCCCCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAAc                  <  2:559779‑M2/213‑1 (MQ=255)
                                                                                                                                                                                                                   atcaacgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAAc                  >  1:74383‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                                        cgttcgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCAGCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGgaga             <  2:756821‑M2/218‑1 (MQ=255)
                                                                                                                                                                                                                            cgcaaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTc         <  1:670198‑M2/222‑1 (MQ=255)
                                                                                                                                                                                                                              caaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTc         <  2:459846‑M2/222‑1 (MQ=255)
                                                                                                                                                                                                                              caaaagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTc         >  1:459846‑M2/28‑249 (MQ=255)
                                                                                                                                                                                                                                 aagaaaatgccgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCAt    >  1:568853‑M2/25‑251 (MQ=255)
                                                                                                                                                                                                                                   gaaaatgacgatttgcggccgcTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCCCTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGt  >  2:847274‑M2/23‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                    
TTTATCCTGTTCTCCTCAGTCGATGATGGGCCAGTTATGGTGGTAAGCTCCCAGCGCTTACATGATATGTTGAATCCTACAAAAGATACCAATTGGAATTCCACGTATATCTACAAATCCAGACATGAGATGTTGCCTGTTAATCTTACTCAGGAAACACTTTTCAGCTCCAAATCTCATGGTAAATATGCGCTTTTCCCCATTTTTACTGCGAGTTGGCGAGCTCATCGTATAATGAATAAGGGTGTTTAAGTAAAGGAAAACATCACCGTTCCTGGCATCCTGGACGGTGATGCCCCTACGGTTGCCCTCGCCAGCACGGGCATCGGTAAAGCGTAAGGTTCAACATCGTTTTACCACTTCATGCGATTGTTGCGTTTTTGTTGCGTATTAGATCACTTAATTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGT  >  NC_000913/4252032‑4252509

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