Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,619,069 Δ18,342 bp rhsB[yhiN] rhsB, yhhH, yhhI, yhhJ, rbbA, yhiI, yhiY, yhiJ, yhiM, [yhiN]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3619069 3637410 18342 65 [0] [0] 65 rhsB–[yhiN] rhsB,yhhH,yhhI,yhhJ,rbbA,yhiI,yhiY,yhiJ,yhiM,[yhiN]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 36190680 (0.000)65 (0.890) 55/492 0.1 100% intergenic (+79/‑124) nikR/rhsB DNA‑binding transcriptional repressor NikR/rhs element protein RhsB
?NC_000913 3637411 = 0 (0.000)intergenic (‑1/‑231) yhiN/pitA putative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA

ACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGG  >  NC_000913/3618823‑3619316
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
ccgACCTCTCCGTCGCCACGCTGCATGTGCCCATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATCCCCTATCCGAgttgc                                                                                                                                                                                                                                                     <  2:393337‑M1/250‑6 (MQ=255)
acGACCTCCCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATTCCCTATCCGAgttgc                                                                                                                                                                                                                                                     <  2:302304‑M1/251‑6 (MQ=255)
   aCCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgct                                                                                                                                                                                                                                                  <  2:191181‑M1/251‑9 (MQ=255)
          cGTCGCCACGCTGCATGTGCACATCAACCACGATGACTGTCTGGAAATCTCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGGCGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctctta                                                                                                                                                                                                                                           <  2:403261‑M1/251‑16 (MQ=255)
                 ccgCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCGAAATACCCTATCCGAgttgcgctcctcttaaaaaaac                                                                                                                                                                                                                                    <  2:458705‑M1/250‑23 (MQ=255)
                 aCGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaac                                                                                                                                                                                                                                    <  1:361352‑M1/251‑23 (MQ=255)
                   gCTGCAAGTGCTCATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTGTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCGATTATTGGGTCCTTTGGTGAACGATCAACAAATGGCCACTAGCACTCCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcataaaaaaacgg                                                                                                                                                                                                                                  <  2:564181‑M1/251‑25 (MQ=255)
                      gCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGGACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacg                                                                                                                                                                                                                               >  1:377617‑M1/1‑224 (MQ=255)
                      gCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacg                                                                                                                                                                                                                               >  1:61689‑M1/1‑224 (MQ=255)
                       cATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgt                                                                                                                                                                                                                              <  1:550880‑M1/251‑29 (MQ=255)
                          gtgCAAATCAACCACGAAGGATTGCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACCCCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtga                                                                                                                                                                                                                            <  1:581116‑M1/250‑31 (MQ=255)
                            gCACATCAACCACGACGACTGTCTGGAAAGCGCCGTGTTGAAAGGGGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGAATTTGCAGTGCTTGCCGAAGGAAGATTGGGTCTATTATTGGGTCCTTTGGTGAACGATCACCAAAGGGCCACTAGCCCACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgcat                                                                                                                                                                                                                          <  2:377064‑M1/250‑33 (MQ=255)
                               caTCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgta                                                                                                                                                                                                                      <  1:392292‑M1/251‑37 (MQ=255)
                                    aCCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagat                                                                                                                                                                                                                 <  1:155773‑M1/251‑42 (MQ=255)
                                        cgacgaCTGTCTGGAAATCGCCGTGTTGAAANGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcct                                                                                                                                                                                                             <  2:272305‑M1/251‑46 (MQ=255)
                                        cgacgaCTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcct                                                                                                                                                                                                             <  1:409819‑M1/251‑46 (MQ=255)
                                         gacgatTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGTCGTTATCGCCCAGCGCGGCGTGCGGCACTGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCTACTAGCACACCAGATGCCTCGCTAAATACCCTATCCGAggtgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctc                                                                                                                                                                                                            <  2:341087‑M1/251‑47 (MQ=255)
                                         gacgacTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGACGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttcaaaaaacggccgtgaattgtaacaattcctc                                                                                                                                                                                                            >  2:214033‑M1/1‑205 (MQ=255)
                                         gacgacTGTCTGGAAAGCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctc                                                                                                                                                                                                            <  2:245582‑M1/251‑47 (MQ=255)
                                            cgacTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcag                                                                                                                                                                                                          <  2:326409‑M1/251‑49 (MQ=255)
                                                    tGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTGGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgc                                                                                                                                                                                                   <  2:407542‑M1/250‑56 (MQ=255)
                                                    tGGAAATCGCCGTGTTGAAAGGTGACATGGGGGACGTGCAGCATTTTGCCGATGACGTAATCGCCCAGCGCGGCGGGCGGCACGGGCATTTGCAGTGCTTGCCGACGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCACTAAATACCCTCTCCGAgttgcgctcctcttaaaaaaacgatcgagaattgtaaagattcatcaggtgacga                                                                                                                                                                                                   >  2:551276‑M1/1‑195 (MQ=255)
                                                     ggAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgc                                                                                                                                                                                                   <  1:184840‑M1/249‑56 (MQ=255)
                                                            gCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATGTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatc                                                                                                                                                                                          <  2:239992‑M1/251‑65 (MQ=255)
                                                            gCCGTGTTGAAAGGTGAAATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatc                                                                                                                                                                                          >  2:246583‑M1/1‑187 (MQ=255)
                                                             ccGTGTTGAAAGGGGACATGGGAGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagatacctcagtagtagcaccagcatc                                                                                                                                                                                          >  2:118382‑M1/1‑186 (MQ=255)
                                                                 ggtgaAAGGTGACAGGGGTGACGTGCAGCATTTTGCCGAAGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATGGAGTCTATTAGGGGGTCCTTTGGGGAAAGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATCCCCTATCCGAgtttcgctcctcttaaaaaaacggacgtgaattgtaaagagtcctcagtggacgcaccagcatcaatat                                                                                                                                                                                     <  2:30967‑M1/249‑70 (MQ=255)
                                                                        ggggacttgggtgacGTGCAGCATTTTGCCGATGACGTTAGCGCCCAGCTCGGCGGGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaa                                                                                                                                                                              <  2:199879‑M1/244‑77 (MQ=255)
                                                                        ggTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaa                                                                                                                                                                              <  1:511487‑M1/251‑77 (MQ=255)
                                                                         gTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaa                                                                                                                                                                             <  1:160145‑M1/251‑78 (MQ=255)
                                                                           gACATGGGTGACGTGCAGCATTTTGCCGATGACGCTATCGCCCAGCGCGGCGTGCGGCACGGGCCTTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctccgtggtcgcaccagcatcaattttactaaaagg                                                                                                                                                                           >  2:93045‑M1/1‑172 (MQ=255)
                                                                           gACATGGGGGACGGGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagaatcaatattactaaaagg                                                                                                                                                                           <  2:477713‑M1/251‑80 (MQ=255)
                                                                            aCATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaagga                                                                                                                                                                          <  1:470088‑M1/251‑81 (MQ=255)
                                                                              aTGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaag                                                                                                                                                                        <  1:381721‑M1/251‑83 (MQ=255)
                                                                                gggTGACGTGCAGCATTTTGCCGATGACGTTATCGACCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagt                                                                                                                                                                       >  1:159004‑M1/1‑167 (MQ=255)
                                                                                  gTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatt                                                                                                                                                                    >  1:507917‑M1/1‑165 (MQ=255)
                                                                                    gACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttg                                                                                                                                                                  <  1:288377‑M1/251‑89 (MQ=255)
                                                                                              aTTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattat                                                                                                                                                        >  1:363604‑M1/1‑153 (MQ=255)
                                                                                                     cGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagat                                                                                                                                                 <  1:295324‑M1/251‑106 (MQ=255)
                                                                                                      gATGACGTTATCGCCCACCGCGGCGTGAGGCACGGGCATTTGCAGTGCTTGCCGAAGGCAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagat                                                                                                                                                 >  1:450890‑M1/1‑145 (MQ=255)
                                                                                                       aTGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGCAGATTGAGTCTATAATTGGGTCCTTTGGTGAACGATCAAAAAAGGGCCACTCGCACACCTGATTCCTCTCTCAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatatcagata                                                                                                                                                >  1:322876‑M1/1‑144 (MQ=255)
                                                                                                         gACGTTATCGCCCAGCGCGGCGTGCGGCACGNGCATTTGCAGTGCTTGCCTAAGGAAGATTTATGCTATTATGGAGTCCTTGGGTGAACGATCTACAAAGGGGCACTATCTCTCCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcattggtctctccagcatcaatattactaaaaggaagtatttgcctaaattatataagataatt                                                                                                                                             <  2:273592‑M1/251‑110 (MQ=255)
                                                                                                             ttATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattt                                                                                                                                         >  1:128742‑M1/1‑138 (MQ=255)
                                                                                                                 cGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattatttttt                                                                                                                                      <  1:122774‑M1/250‑117 (MQ=255)
                                                                                                                 cGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattatttttt                                                                                                                                      <  1:509096‑M1/250‑117 (MQ=255)
                                                                                                                  gCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttga                                                                                                                                    <  1:39333‑M1/251‑119 (MQ=255)
                                                                                                                    ccAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagt                                                                                                                                  <  1:539046‑M1/251‑121 (MQ=255)
                                                                                                                             cgTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattgtttgagtgaaatccat                                                                                                                         <  1:411344‑M1/251‑130 (MQ=255)
                                                                                                                              gTGCGGCACTGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACCCCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaagtgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccata                                                                                                                        <  2:127092‑M1/251‑131 (MQ=255)
                                                                                                                                                      ttGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtct                                                                                                >  1:79849‑M1/1‑97 (MQ=255)
                                                                                                                                                           gaaggaagATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatatttctaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatt                                                                                           <  1:487288‑M1/251‑160 (MQ=255)
                                                                                                                                                              ggaagATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttca                                                                                        <  1:198185‑M1/251‑163 (MQ=255)
                                                                                                                                                                            attattGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccg                                                                          >  1:567182‑M1/1‑75 (MQ=255)
                                                                                                                                                                                ttGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCAAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcgtcaatattactaaaagga                                                                                                                                                                          <  2:493825‑M1/151‑81 (MQ=255)
                                                                                                                                                                                ttGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCAAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaagga                                                                                                                                                                          >  1:493825‑M1/1‑71 (MQ=255)
                                                                                                                                                                                              gAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTAACCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatg                                                        <  1:214033‑M1/251‑195 (MQ=255)
                                                                                                                                                                                                   aTCAACAAAGGGCCACTAGCACACCTGATTCNCCGCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaatgaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatgcgccg                                                   <  2:290357‑M1/251‑200 (MQ=255)
                                                                                                                                                                                                   aTCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatgcgccg                                                   >  1:254161‑M1/1‑52 (MQ=255)
                                                                                                                                                                                                      aacaacGGGCCACTAGCTCATCTGATTCCTCTCTAAATCCCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatttcagataattattttttgagtgaaatccatacagggggcaaatcaaacaaagtctatatttcactttgcctgcgccgcgaaagtcactgataatgcgccgcgt                                                >  2:42967‑M1/1‑49 (MQ=255)
                                                                                                                                                                                                         aaaGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatgcgccgcgttc                                              >  1:432426‑M1/1‑46 (MQ=255)
                                                                                                                                                                                                                    aGCACACCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatgcgccgcgttcatgtcctcaaa                                   <  1:573850‑M1/250‑216 (MQ=255)
                                                                                                                                                                                                                      cacaCCTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatgcgccgcgttcatgtcctcaaaatg                                <  2:377617‑M1/251‑219 (MQ=255)
                                                                                                                                                                                                                       ccaccTGATTCCTCTCTAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatgcgccgcgttcatgtcctcaaaatgg                               <  1:13212‑M1/250‑220 (MQ=255)
                                                                                                                                                                                                                                      ctAAATACCCTATCCGAgttgcgctcctcttaaaaaaacggacgtgaattgtaaagattcctcagtggtcgcaccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatccatacagggggcaaatcaaaaaaagtctatatttcactttgcccgcgccgcgaaagtcactgataatgcgccgcgttcatgtcctcaaaatggcgtaacgtcctatgc                >  2:541005‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                                    cGAgttgcgctactcttaaaaaaacggacgtgaattgtacagattcctcagtggtcgccccagcatcaatattactaaaaggaagtatttgcctgaattatataagataattattttttgagtgaaatacataccgggtgcaaatcaaaaaaagtctatttttcactttgcccgctccgcgaaagtcactgataatgcgccccgttcatgtcctcacaantgcggaacgccctatgctacctttgtttgct  >  2:279237‑M1/1‑3 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        
ACGACCTCTCCGTCGCCACGCTGCATGTGCACATCAACCACGACGACTGTCTGGAAATCGCCGTGTTGAAAGGTGACATGGGTGACGTGCAGCATTTTGCCGATGACGTTATCGCCCAGCGCGGCGTGCGGCACGGGCATTTGCAGTGCTTGCCGAAGGAAGATTGAGTCTATTATTGGGTCCTTTGGTGAACGATCAACAAAGGGCCACTAGCACACCTGATTCCTCTCTAAATACCCTATCCGAACCTTCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGG  >  NC_000913/3618823‑3619316

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