Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 3,992,517 Δ5 bp 100% coding (1365‑1369/2547 nt) cyaA → adenylate cyclase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3992517 3992521 5 104 [0] [0] 103 cyaA adenylate cyclase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3992522 =0 (0.000)97 (1.060) 66/494 NT 100% coding (1370/2547 nt) cyaA adenylate cyclase
?NC_000913 = 3992516 0 (0.000)coding (1364/2547 nt) cyaA adenylate cyclase

Only 100 of 103 total aligned reads displayed.

AGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTG  >  NC_000913/3992238‑3992808
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                              
gacgaagctcgtctggcaatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCgg                                                                                                                                                                                                                                                                                            <  1:125651‑M1/5‑1 (MQ=255)
    aagctcgtctggcaatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTa                                                                                                                                                                                                                                                                                        <  1:610931‑M1/9‑1 (MQ=255)
    aagctcgtctggcaatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTa                                                                                                                                                                                                                                                                                        <  2:640664‑M1/9‑1 (MQ=255)
    aagctcgtctggcaatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTa                                                                                                                                                                                                                                                                                        <  2:739789‑M1/9‑1 (MQ=255)
          gtctggcaatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGtt                                                                                                                                                                                                                                                                                  >  1:94467‑M1/275‑289 (MQ=255)
            ctgggaatgctcagtaaccgggctaactggaagattggtcggtggcggggggcggacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgcttggcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgcgccGCCGGGCTAACCGTTc                                                                                                                                                                                                                                                                                 >  2:1037‑M1/273‑288 (MQ=255)
            ctggcaatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTc                                                                                                                                                                                                                                                                                 >  1:170718‑M1/273‑288 (MQ=255)
                caatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGt                                                                                                                                                                                                                                                                             >  2:264153‑M1/269‑288 (MQ=255)
                 aatgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGtt                                                                                                                                                                                                                                                                            >  1:739789‑M1/268‑288 (MQ=255)
                   tgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         >  2:367726‑M1/266‑289 (MQ=255)
                   tgctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         >  1:740268‑M1/266‑289 (MQ=255)
                     ctcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         >  2:629983‑M1/264‑287 (MQ=255)
                      tcgataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTAt                                                                                                                                                                                                                                                                       >  2:673853‑M1/263‑288 (MQ=255)
                      tcgataaccgcgctaactggaagatggatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTAt                                                                                                                                                                                                                                                                       >  1:611618‑M1/263‑288 (MQ=255)
                      tccataaccgcgctaatgggagaattgatcaggggggcgagtcgcacaacgtgttgcgctacggtagtctgcagggctaccgtaatctcatcctctttgcgcgtcgcaaaaaccttagcgtcgccgccagtccgcatgatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTAt                                                                                                                                                                                                                                                                       >  2:309088‑M1/263‑288 (MQ=255)
                         ataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         >  1:842303‑M1/260‑283 (MQ=255)
                         ataaccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         <  2:842303‑M1/24‑1 (MQ=255)
                             ccgcgctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGtata                                                                                                                                                                                                                                                                <  1:412879‑M1/33‑1 (MQ=255)
                                  ctaactggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTCTCTGTATAACcgc                                                                                                                                                                                                                                                           <  2:659099‑M1/38‑1 (MQ=255)
                                      ctggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGcc                                                                                                                                                                                                                                                      >  2:137588‑M1/247‑289 (MQ=255)
                                       tggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCa                                                                                                                                                                                                                                                     <  1:171116‑M1/44‑1 (MQ=255)
                                       gggaagagtgatcaggtgcgtggggcgcacaacgagttgctcgacgcgatgatgcagagctagcgtaatctgctccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGcc                                                                                                                                                                                                                                                      >  2:509965‑M1/246‑288 (MQ=255)
                                        ggaagattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCaa                                                                                                                                                                                                                                                    >  2:144278‑M1/245‑289 (MQ=255)
                                           agattgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAata                                                                                                                                                                                                                                                 >  1:866780‑M1/242‑289 (MQ=255)
                                               tgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         <  1:146781‑M1/24‑1 (MQ=255)
                                               tgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         >  2:146781‑M1/238‑261 (MQ=255)
                                               tgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctcccggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTggtatctgtataaccgccgccccaatcatga                                                                                                                                                                                                                                             <  2:216800‑M1/52‑15 (MQ=255)
                                               tgatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgaactctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCACATATtgg                                                                                                                                                                                                                                             <  1:337052‑M1/52‑2 (MQ=255)
                                                gatcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAg                                                                                                                                                                                                                                            <  2:842122‑M1/53‑1 (MQ=255)
                                                 atcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTTCGGTTGGTATCTGTACAACCGCGCGCCAAATAttggg                                                                                                                                                                                                                                            <  2:675185‑M1/53‑3 (MQ=255)
                                                 atcaggtgcgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAg                                                                                                                                                                                                                                            <  1:166697‑M1/53‑1 (MQ=255)
                                                  gaaggggcgggcggcgcacaacggggtgctgacagcagtgctgcagagctaccgtaacctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTc                                                                                                                                                                                                                                          >  2:500701‑M1/235‑289 (MQ=255)
                                                         cgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggctatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgactctgacctttctttaggtgccGCCGGGCTAACCGTTCAGGCTGGTATCTGTATACCCGCGCGCCAAATATTGAGtcgatgct                                                                                                                                                                                                                                    <  2:721542‑M1/61‑4 (MQ=255)
                                                         cgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         <  2:803800‑M1/24‑1 (MQ=255)
                                                         cgtgaggcgcacaacgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGt                                                                                                                                                                                                                                                                         >  1:803800‑M1/228‑251 (MQ=255)
                                                                      acgagttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGc                                                                                                                                                                                                                      >  1:842122‑M1/215‑289 (MQ=255)
                                                                         agttgctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTg                                                                                                                                                                                                                    >  2:778642‑M1/212‑288 (MQ=255)
                                                                             gctcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAata                                                                                                                                                                                                               >  1:675185‑M1/208‑289 (MQ=255)
                                                                               tcgacgcgatgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggtaccgaatttgacctttatttatgtgccGCCGGTCTAACCGTTCAGGTTGGTACTTGTATCACCGCCCTCCCAATATTTCATCGATACTCAGCCATCAGCCTCTTGAAtata                                                                                                                                                                                                             <  2:384109‑M1/84‑1 (MQ=255)
                                                                                        tgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTa                                                                                                                                                                                                      <  1:767875‑M1/91‑1 (MQ=255)
                                                                                        tgatgcagagctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTAc                                                                                                                                                                                                     >  1:773130‑M1/197‑288 (MQ=255)
                                                                                             cagaggtagcgtaagctgatccggtttgcgcggcgcaatacccttagcgtctcggccagtccgcaggatgtcggcgtgctgacgggtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCAGCAGCCATCAGCCGCTGGAATATAACCGTTACCTGaa                                                                                                                                                                                                >  2:495811‑M1/192‑288 (MQ=255)
                                                                                                 gctcccgtagcctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATa                                                                                                                                                                                              >  2:166697‑M1/188‑286 (MQ=255)
                                                                                                  ctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTc                                                                                                                                                                                                                                                                                 <  2:253105‑M1/16‑1 (MQ=255)
                                                                                                  ctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTc                                                                                                                                                                                                                                                                                 >  1:253105‑M1/187‑202 (MQ=255)
                                                                                                  ctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTAt                                                                                                                                                                                                                                                                       >  1:34078‑M1/187‑212 (MQ=255)
                                                                                                  ctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTAt                                                                                                                                                                                                                                                                       <  2:34078‑M1/26‑1 (MQ=255)
                                                                                                  ctaccgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATaa                                                                                                                                                                                             <  1:232651‑M1/100‑1 (MQ=255)
                                                                                                   taccgtaatctgatccgcattgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtg                                                                                                                                                                                         >  1:411683‑M1/186‑289 (MQ=255)
                                                                                                    accgtaatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtgg                                                                                                                                                                                        <  1:281919‑M1/105‑1 (MQ=255)
                                                                                                         aatctgatccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGc                                                                                                                                                                                    >  1:811041‑M1/180‑288 (MQ=255)
                                                                                                               atccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGc                                                                                                                                                                                                                      <  1:866178‑M1/75‑1 (MQ=255)
                                                                                                                tccggtgtgcgcgtagcattaacctgagcgtctcctccagtcggaagtagataggcgtgctggagcgtaagcggtaggctgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCa                                                                                                                                                                             >  2:478124‑M1/173‑288 (MQ=255)
                                                                                                                tccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCa                                                                                                                                                                             >  1:216800‑M1/173‑288 (MQ=255)
                                                                                                                tccgctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCa                                                                                                                                                                             >  1:640664‑M1/173‑288 (MQ=255)
                                                                                                                gccgcgttgcgctgcgcgagaaccttcgcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatcttgcctttatttatgtgccGCCGTGCTAACCGTTGAGGTTGGTATTTGTATAACCGCGCGCCAAATATTGAGTCAATCATCAGCCATAAGCCGCTGGAATATAACCGTTACCTGAATAAACTGATGGCGTGGGCAt                                                                                                                                                                            >  2:337052‑M1/173‑289 (MQ=255)
                                                                                                                    ctttgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGtt                                                                                                                                                                        <  1:67653‑M1/121‑1 (MQ=255)
                                                                                                                       tgcgcgtcgcaataaccttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTa                                                                                                                                                                      <  1:154138‑M1/123‑1 (MQ=255)
                                                                                                                                        ttagcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgc                                                                                                                                                    <  1:290808‑M1/141‑1 (MQ=255)
                                                                                                                                          agcgtctccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTcgcg                                                                                                                                                   <  2:272122‑M1/142‑1 (MQ=255)
                                                                                                                                                tccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                                                                            >  2:154138‑M1/141‑289 (MQ=255)
                                                                                                                                                tccgccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                                                                            >  2:226452‑M1/141‑289 (MQ=255)
                                                                                                                                                    ccagtccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGt                                                                                                                                        >  1:465443‑M1/137‑289 (MQ=255)
                                                                                                                                                        tccgcaggatatcggcgtgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTatat                                                                                                                                    <  1:21304‑M1/157‑1 (MQ=255)
                                                                                                                                                                         tgctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCTCGCGCCAAATATTGAGTCGATCATCAGCCATCCGCCGCTGGAATATAACCGTTATCTGAATAAACTGGTGGCGTTGGCCTGGTTTAACGGCCTGCTGACCTTGCGCACCCGCGTGTATATTAAAGGTAACGGCATTg                                                                                                                   <  2:545861‑M1/174‑1 (MQ=255)
                                                                                                                                                                          gctgacgcgtaagctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTg                                                                                                                   >  2:767875‑M1/115‑288 (MQ=255)
                                                                                                                                                                                     agctgtatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTa                                                                                                       >  2:67653‑M1/104‑289 (MQ=255)
                                                                                                                                                                                          tatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAAcc                                                                                                                                                                                                          >  1:423698‑M1/99‑185 (MQ=255)
                                                                                                                                                                                          tatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAAcc                                                                                                                                                                                                          <  2:423698‑M1/87‑1 (MQ=255)
                                                                                                                                                                                          tatgccgcgtttgaagcattaccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGtt                                                                                                   <  1:388096‑M1/190‑1 (MQ=255)
                                                                                                                                                                                                 cgttcgaagcattaccaggtaaagtgacgctggtacacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTCCAGGTTGGTATCTGTATACCCGCGCGCCAACTATTGAGTCGACCATCAGCCATCAGACGCTGGAATATAACCGTTACCTGACTAAACTGTTCGCGTAGGCATGGCTTAATGCCCTGCAGACCTGGTGCCTCCGCTTGTATACTAAAGGCAAGGTCAGTGTCGATTTGCctaaattgatgcgca                                                                                           <  2:565119‑M1/198‑8 (MQ=255)
                                                                                                                                                                                                             taccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCTCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGTCTAAGTTGCAGGAGATGGTCGCCGACg                                                                               <  1:539870‑M1/210‑1 (MQ=255)
                                                                                                                                                                                                             taccaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACg                                                                               <  1:208896‑M1/210‑1 (MQ=255)
                                                                                                                                                                                                              accaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAata                                                                                                                                                                                                               >  2:588476‑M1/79‑160 (MQ=255)
                                                                                                                                                                                                              accaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAata                                                                                                                                                                                                               <  1:588476‑M1/82‑1 (MQ=255)
                                                                                                                                                                                                              accaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTaa                                                                                                                           <  2:214179‑M1/166‑1 (MQ=255)
                                                                                                                                                                                                              accaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTaa                                                                                                                           >  1:214179‑M1/79‑244 (MQ=255)
                                                                                                                                                                                                              accaggtaaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACgt                                                                              <  1:628534‑M1/211‑1 (MQ=255)
                                                                                                                                                                                                                    taaagtgacgctggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCg                                                                          >  1:527270‑M1/73‑287 (MQ=255)
                                                                                                                                                                                                                               tggtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCg                                                              >  2:290808‑M1/62‑288 (MQ=255)
                                                                                                                                                                                                                                 gtaaacccgcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTTCAGGAAATGGTCGGCGACGTGTCGCACCTtttccccgc                                                            <  2:782585‑M1/229‑4 (MQ=255)
                                                                                                                                                                                                                                   aaaccccctgatttcacccgatgtctcggaaccgaatctgacctttatttatgggccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCAGCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTgcg                                                         >  2:412879‑M1/58‑289 (MQ=255)
                                                                                                                                                                                                                                         gcagatttcacccgatctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTAcc                                                   >  2:388096‑M1/52‑289 (MQ=255)
                                                                                                                                                                                                                                                       atctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                     <  1:202354‑M1/252‑1 (MQ=255)
                                                                                                                                                                                                                                                        tctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                     >  2:208896‑M1/37‑288 (MQ=255)
                                                                                                                                                                                                                                                        tcgctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                     >  1:734782‑M1/37‑288 (MQ=255)
                                                                                                                                                                                                                                                         ctctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATGGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGGTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAg                                   >  2:171116‑M1/36‑289 (MQ=255)
                                                                                                                                                                                                                                                           ttcggaacctaatctgaccttgatttattggcgGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgc                                 >  2:122192‑M1/34‑289 (MQ=255)
                                                                                                                                                                                                                                                           ctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgc                                 >  1:384109‑M1/34‑289 (MQ=255)
                                                                                                                                                                                                                                                           ctcggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgc                                 >  1:174680‑M1/34‑289 (MQ=255)
                                                                                                                                                                                                                                                              ggaaccgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgcgc                               <  1:714022‑M1/257‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cgaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGCCCCCGAAGGCGCCCTACa                         <  2:684915‑M1/264‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gaatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCAGCTGCGCTTACCTGCGCCGACACCGAAGGCGCCCTACa                         <  1:512642‑M1/264‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAt                                                                                                              <  2:511403‑M1/179‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAt                                                                                                              >  1:511403‑M1/24‑202 (MQ=255)
                                                                                                                                                                                                                                                                     aatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                     <  2:40412‑M1/252‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aatctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                     >  1:40412‑M1/24‑275 (MQ=255)
                                                                                                                                                                                                                                                                        ctgacctttatttatgtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCgt                   >  2:714022‑M1/21‑289 (MQ=255)
                                                                                                                                                                                                                                                                                       gtgccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGAACCCGAAGGCGCTCTACAGCCCGGGTGAACTCcgcaa      <  2:189318‑M1/283‑3 (MQ=255)
                                                                                                                                                                                                                                                                                          ccGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTg  <  1:811700‑M1/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                              
AGCTCGTCTGGCAATGCTCGATAACCGCGCTAACTGGAAGATTGATCAGGTGCGTGAGGCGCACAACGAGTTGCTCGACGCGATGATGCAGAGCTACCGTAATCTGATCCGCTTTGCGCGTCGCAATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTG  >  NC_000913/3992238‑3992808

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