Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,620,828 43 bp→AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT intergenic (+2/‑13) deoB → / → deoD phosphopentomutase/purine nucleoside phosphorylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4620828 4620870 43 62 [0] [0] 63 deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46208270 (0.000)47 (0.710)
+AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT
32/458 0.7 100% intergenic (+1/‑56) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase
?NC_000913 4620871 = 0 (0.000)intergenic (+45/‑12) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

TTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTT  >  NC_000913/4620575‑4621153
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                          
gatgtctctgctgcgcgatgacgacatcctgatcctcactgctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataa                                                                                                                                                                                                                                                                                           >  1:350441‑M2/297‑298 (MQ=255)
        tgctgcgcgattacgacatcctgatcctcactgctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgggaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacaggtagagggtgaaaacatgttgatAAGGATAAAAc                                                                                                                                                                                                                                                                                  <  2:309875‑M2/11‑1 (MQ=255)
        tgctgcgcgatgacgacatcctgatcctcactgctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAAc                                                                                                                                                                                                                                                                                  <  2:630238‑M2/11‑1 (MQ=255)
        tgctgcgcgatgacgacatcctgatcctcactgctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgagAAGAATAAAAc                                                                                                                                                                                                                                                                                  <  2:690126‑M2/11‑1 (MQ=255)
                  tgacgacatcctgatcctcactgctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTAccc                                                                                                                                                                                                                                                                       <  1:366035‑M2/22‑1 (MQ=255)
                                   tcactgctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgctAAGGCTAAAACAATGGCTACCCCACACATTAATGCAGaa                                                                                                                                                                                                                                                      >  2:308233‑M2/262‑300 (MQ=255)
                                      ctgctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctaccacagggagagggttaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAAt                                                                                                                                                                                                                                                    >  1:31632‑M2/259‑299 (MQ=255)
                                        gctgaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAAt                                                                                                                                                                                                                                                    >  1:694730‑M2/257‑297 (MQ=255)
                                           gaccacggttgtgatccgacctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGa                                                                                                                                                                                                                                              >  1:425898‑M2/254‑300 (MQ=255)
                                                             acctggaccggtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcacccaggccctctccaacagggagagggtgaaaacatgttgctAAGGATAAAACAATGGCTAACCCACCCATTAATGCAGAAATGGGCGATTTCGCTGACGTAGtttt                                                                                                                                                                                                                            >  2:65824‑M2/236‑300 (MQ=255)
                                                             acctggaccagtactgaccacacgcgtgaacacattccggtactggtttacggcacgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtatttctgatatggaatatggcaaagccacgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAGCACTGGCTAACCCACACATTAATGCAGAGATGGGCGATTTCGCTGACGtagtctt                                                                                                                                                                                                                            >  2:65516‑M2/236‑297 (MQ=255)
                                                                       gtactgaccacacgcgtgaacacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGc                                                                                                                                                                                                                  >  1:298990‑M2/226‑300 (MQ=255)
                                                                                acacgcgggaacacatgccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacaggtagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtg                                                                                                                                                                                                         <  2:36016‑M2/84‑1 (MQ=255)
                                                                                          acacattccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCAAGGCGAACAGCTGCGTGCGAAGt                                                                                                                                                                                               >  2:276256‑M2/207‑300 (MQ=255)
                                                                                              attccggtactggtttacggcccgaaagtaaaaccgggctcactggttcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccngccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTatat                                                                                                                                                                                           <  2:497666‑M2/98‑1 (MQ=255)
                                                                                               ttccggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATAtt                                                                                                                                                                                          >  1:660771‑M2/202‑300 (MQ=255)
                                                                                               ttccggtactggtttacggcacgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTAACCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATAtt                                                                                                                                                                                          >  1:660756‑M2/202‑300 (MQ=255)
                                                                                                   ggtactggtttacggcccgaaagtaaaaccgggctcactgggtcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGc                                                                                                                                                                                        <  1:655740‑M2/101‑1 (MQ=255)
                                                                                                                              accgggctcactggttcaccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcagagccatgttctgataatgggtgcggcgtgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCATGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCgtgaag                                                                                                                                                           <  2:178212‑M2/130‑1 (MQ=255)
                                                                                                                                       actggttcaccgggaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcacccctgccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCTAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGaacaac                                                                                                                                                   <  2:110904‑M2/138‑1 (MQ=255)
                                                                                                                                               accgtgaaaccttcgcngatattggtcagactctgncaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAANTGAACAACGTTCGCgg                                                                                                                                           <  1:586192‑M2/146‑1 (MQ=255)
                                                                                                                                               accgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGt                                                                                                                                          <  1:1280‑M2/147‑1 (MQ=255)
                                                                                                                                                ccgtgaaaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTa                                                                                                                                         >  1:472643‑M2/153‑300 (MQ=255)
                                                                                                                                                          ttcgcggatattgtgcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacaggtagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGc                                                                                                                                 <  2:163842‑M2/156‑1 (MQ=255)
                                                                                                                                                          ttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcacccctgccctctcccacagggagagggttaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGc                                                                                                                                 <  2:31632‑M2/156‑1 (MQ=255)
                                                                                                                                                          ttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACNCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGc                                                                                                                                 <  2:507910‑M2/156‑1 (MQ=255)
                                                                                                                                                           ttgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctggtaatgggtgcggcttgatgccctcaccccgtccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACTTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCt                                                                                                                                <  2:180321‑M2/157‑1 (MQ=255)
                                                                                                                                                                 atattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctcccccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGt                                                                                                                        <  2:684870‑M2/165‑1 (MQ=255)
                                                                                                                                                                       gtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTAc                                                                                                                  >  1:46310‑M2/130‑300 (MQ=255)
                                                                                                                                                                              tctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGc                                                                                                            <  1:25547‑M2/177‑1 (MQ=255)
                                                                                                                                                                                                       gatatggaatatggcaaagccatgttctgataatggttgctgcttaatgcccgcaccccggccctctcccacagggagagggtgaaaacacgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCATAAGTGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTTAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACgg                                                                                  <  2:168575‑M2/203‑1 (MQ=255)
                                                                                                                                                                                                       gatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCAc                                                                                    >  1:419641‑M2/98‑298 (MQ=255)
                                                                                                                                                                                                       gatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCAc                                                                                    >  2:480910‑M2/98‑298 (MQ=255)
                                                                                                                                                                                                        atatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacctatcgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGAAGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGCGAACAACGTTCGCGGTATGCTGGGCTTCAACGGTACTTACAAAGGACGCAAAATTTCCGTAATGGGTCACGGt                                                                                 >  2:297719‑M2/97‑300 (MQ=255)
                                                                                                                                                                                                        atatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGt                                                                                 >  2:297710‑M2/97‑300 (MQ=255)
                                                                                                                                                                                                            ggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATGAATGCAGAAATTGGCGATTTCGCTGACGTAGTTTTGATGCCATGCGACCCGCTGCGTGCGAAGGATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATgg                                                                             <  2:64729‑M2/208‑1 (MQ=255)
                                                                                                                                                                                                            ggaatatgacaaagccatgttcgaataatgggtgcggcttgatgcccccaccccggccccctcccacagggagagggtgaaaacatgttgagAAGGATAAAACAACGGCTACCCCACACATTAAGGCAGAAATGGGCGATTTCGCGGACGTAGTTTTGATTACATGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATgg                                                                             <  2:64569‑M2/208‑1 (MQ=255)
                                                                                                                                                                                                             gaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTAt                                                                               >  1:688147‑M2/92‑297 (MQ=255)
                                                                                                                                                                                                             gaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATggg                                                                            <  2:639159‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                             gaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCAACGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCAAGGTAt                                                                               >  1:95477‑M2/92‑297 (MQ=255)
                                                                                                                                                                                                             gaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGATGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTAt                                                                               >  2:463427‑M2/92‑297 (MQ=255)
                                                                                                                                                                                                                   ggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATccc                                                                      <  1:365963‑M2/215‑1 (MQ=255)
                                                                                                                                                                                                                       aagccatgttctgataatgggtgcggcttgatgccctcacccctgccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTcc                                                                  <  1:648422‑M2/219‑1 (MQ=255)
                                                                                                                                                                                                                         gccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCACGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCccgtact                                                                 >  2:690685‑M2/80‑296 (MQ=255)
                                                                                                                                                                                                                         gccatgttctgataatggatgcggctcgatgccctcaccccggccctctcccagagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCACGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCAGAAAAATTTCCGTAATGGGTCACGGTATgggtatccagtact                                                                 >  2:12759‑M2/80‑293 (MQ=255)
                                                                                                                                                                                                                                   gataatgggtgcggcttgatgccctcacccctgccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGTCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCCGACGTAGTTGTGATGCCAGGCGACCCACTGCGTGCGAAGTATATTGCTTAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTac                                                      <  2:405995‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                                   gataatgggtgcggcttgatgccctcaccccggccctctaccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTac                                                      >  1:87395‑M2/70‑300 (MQ=255)
                                                                                                                                                                                                                                      aatgggtgcggcttgatgccctcaccccggccctttcccacaggtagagggtgaaaacatgttgatAAGGAGAAAACAATGGCTACCCCACACCTTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCTTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACAcc                                                   <  2:26344‑M2/234‑1 (MQ=255)
                                                                                                                                                                                                                                          ggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAg                                               <  1:472995‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                                gcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTACCCGGTACTTACAAAGGCCGCAAAATTTCTGTAATGGGTCACGGTATGGGTATCCAGTCCTGCTCCATCTACACCAAAGAACt                                           >  2:22439‑M2/57‑298 (MQ=255)
                                                                                                                                                                                                                                                  ttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATc                                       >  1:195034‑M2/55‑300 (MQ=255)
                                                                                                                                                                                                                                                        ccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGCAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAACTTTCCGTCATGGGTCCCGGTCTGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAACGAt                                 >  2:53028‑M2/49‑300 (MQ=255)
                                                                                                                                                                                                                                                        ccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTCATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAt                                 >  1:451735‑M2/49‑300 (MQ=255)
                                                                                                                                                                                                                                                        ccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAt                                 >  1:76839‑M2/49‑300 (MQ=255)
                                                                                                                                                                                                                                                        ccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCAACGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCAGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAt                                 >  1:76351‑M2/49‑300 (MQ=255)
                                                                                                                                                                                                                                                              ccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCAGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGc                           >  1:280857‑M2/43‑300 (MQ=255)
                                                                                                                                                                                                                                                                      ctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGNCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAACGATTTCGGCGTgaagaa                   >  2:540984‑M2/35‑300 (MQ=255)
                                                                                                                                                                                                                                                                      ctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTgaagaa                   >  1:233356‑M2/35‑300 (MQ=255)
                                                                                                                                                                                                                                                                        ctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGaaaa                 >  1:62472‑M2/33‑300 (MQ=255)
                                                                                                                                                                                                                                                                        ctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCACGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGaaaa                 >  1:194360‑M2/33‑300 (MQ=255)
                                                                                                                                                                                                                                                                          cccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAAtt               >  2:565540‑M2/31‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      gggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGt   >  1:30738‑M2/19‑300 (MQ=255)
                                                                                                                                                                                                                                                                                       ggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGtt  <  1:129909‑M2/283‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                          
TTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTT  >  NC_000913/4620575‑4621153

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