Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,313,864 Δ16,405 bp yjdPrdcB 18 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4313864 4330268 16405 74 [0] [0] 73 yjdP–rdcB yjdP,phnP,phnO,phnN,phnM,phnL,phnK,phnJ,phnI,phnH,phnG,phnF,phnD,phnC,yjdN,yjdM,rdcA,rdcB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 43138630 (0.000)73 (0.950) 67/492 0.1 100% intergenic (+64/‑5) rpiB/yjdP allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B/protein YjdP
?NC_000913 4330269 = 0 (0.000)intergenic (+31/‑233) rdcB/proP regulator of diguanylate cyclase RdcB/osmolyte:H(+) symporter ProP

ACCATCAGTGGTTTGCCTCTTCGACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAG  >  NC_000913/4330021‑4330512
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
accttattccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgaaa                                                                                                                                                                                                                                                   <  1:433035‑M2/3‑1 (MQ=255)
 ccttattccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggaaagagtttctttttgaggttgAAAt                                                                                                                                                                                                                                                  >  1:85202‑M2/248‑251 (MQ=255)
        ccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCAtgtg                                                                                                                                                                                                                                           >  1:55263‑M2/241‑251 (MQ=255)
            gctactttcgcggcagcataacgactccaacttgctggcttttggttcacgagtggttggactcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgatttcggcaatagagcagcggagaaattgagattcatcctctacttgcatgtatgagtaatgattaatgtggatagagtttctttttgaggtttAAATCCATGTGTGaa                                                                                                                                                                                                                                       <  2:87882‑M2/15‑1 (MQ=255)
            gcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGa                                                                                                                                                                                                                                        >  1:503540‑M2/237‑250 (MQ=255)
               actttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGtt                                                                                                                                                                                                                                    <  1:204672‑M2/18‑1 (MQ=255)
                  ttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTg                                                                                                                                                                                                                                   <  1:243266‑M2/19‑1 (MQ=255)
                  ttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGAt                                                                                                                                                                                                                                 >  2:189900‑M2/231‑251 (MQ=255)
                    cgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATc                                                                                                                                                                                                                                <  1:387943‑M2/22‑1 (MQ=255)
                     gcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATcac                                                                                                                                                                                                                              <  1:39735‑M2/24‑1 (MQ=255)
                           gcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAAttt                                                                                                                                                                                                                        <  1:152249‑M2/30‑1 (MQ=255)
                              taacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtaagaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTaa                                                                                                                                                                                                                      >  2:620520‑M2/219‑250 (MQ=255)
                                acgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAac                                                                                                                                                                                                                    >  1:376167‑M2/217‑250 (MQ=255)
                                  gacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACt                                                                                                                                                                                                                 >  2:439424‑M2/215‑251 (MQ=255)
                                  gacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACt                                                                                                                                                                                                                 >  1:292638‑M2/215‑251 (MQ=255)
                                   acaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtgtatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTg                                                                                                                                                                                                                <  1:239945‑M2/38‑1 (MQ=255)
                                   acaccaacgtgctggctgttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaaggttgatagagtttctttttgaggttgAAATCCAGGTGTGAAGTTGATCACAAATTTAAACACTg                                                                                                                                                                                                                <  2:206233‑M2/38‑1 (MQ=255)
                                    caccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTgg                                                                                                                                                                                                               >  1:590471‑M2/213‑251 (MQ=255)
                                                ggcttttgggtcccgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcgtccactacttgaatggatgagtaatgattaatgtggatagaggttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGAaaaa                                                                                                                                                                                                     >  2:606802‑M2/201‑249 (MQ=255)
                                                      tggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCAt                                                                                                                                                                                              <  2:109773‑M2/56‑1 (MQ=255)
                                                       ggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCAtt                                                                                                                                                                                             <  1:457879‑M2/57‑1 (MQ=255)
                                                       ggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTa                                                                                                                                                                                            <  1:163709‑M2/58‑1 (MQ=255)
                                                        gttcacgagtggatggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTa                                                                                                                                                                                            >  2:61521‑M2/193‑250 (MQ=255)
                                                             cgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGcc                                                                                                                                                                                      <  1:131454‑M2/64‑1 (MQ=255)
                                                                gtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCaa                                                                                                                                                                                   <  1:499883‑M2/67‑1 (MQ=255)
                                                                       gcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGc                                                                                                                                                                            <  2:292638‑M2/74‑1 (MQ=255)
                                                                           cgaactggcaaaagtgattgcggatgcgtagctgggcgcacagaacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggctagagtttctttttgaggttgAAATCCATGTGTGTAGTTGATCACAAATTTAACCACTGGTAGGGAAAAAAGGTCATTAACTGCCCAATTCAAGCGTc                                                                                                                                                                         >  2:450265‑M2/174‑250 (MQ=255)
                                                                           cgaacgtgcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCa                                                                                                                                                                        <  2:247579‑M2/78‑1 (MQ=255)
                                                                            gaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAACGGTCATTAACGGCCCAATTCAGGCGTCa                                                                                                                                                                        >  1:191771‑M2/173‑250 (MQ=255)
                                                                                        atgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTg                                                                                                                                                           <  1:30851‑M2/91‑1 (MQ=255)
                                                                                         tgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgcttacggcaatagagcagcggagaaattgagattcctccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGATGATCACCAATTTAAACACTGGTAGGGAAAAAAGGTCATTAACTGCCAAATTCAGGCGTCAACTGGGTTGATTGt                                                                                                                                                          >  2:31422‑M2/160‑251 (MQ=255)
                                                                                                        ggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCg                                                                                                                                            <  1:449108‑M2/106‑1 (MQ=255)
                                                                                                                   cagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGc                                                                                                                                >  2:148501‑M2/134‑251 (MQ=255)
                                                                                                                        cgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCAGAGGGTGTAAGCAAAcc                                                                                                                           <  1:475133‑M2/123‑1 (MQ=255)
                                                                                                                         gaaggcggtcgccctaaaccccgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcctggatgagtaatgaataatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTACACACTGGTAGGGAAAAAAGGTCATTCACTGCCCAATTCAGGCGTCAACTGGATTGATTGAACATTCCTTAACCGAAGGGTGTAAGCAAAccc                                                                                                                          >  2:550879‑M2/128‑251 (MQ=255)
                                                                                                                          aaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCg                                                                                                                         >  1:68887‑M2/127‑251 (MQ=255)
                                                                                                                           aggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGc                                                                                                                        >  1:578848‑M2/126‑251 (MQ=255)
                                                                                                                             gcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTa                                                                                                                      <  1:550879‑M2/128‑1 (MQ=255)
                                                                                                                              cggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTAc                                                                                                                     >  2:621847‑M2/123‑251 (MQ=255)
                                                                                                                              cggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatggggatagagtttctttttgaggtttAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCGTAACCGGAGGGTGTAAGCAAACCCGCTAc                                                                                                                     <  2:107039‑M2/129‑1 (MQ=255)
                                                                                                                                  cgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCtt                                                                                                                 <  1:148501‑M2/133‑1 (MQ=255)
                                                                                                                                  cgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCtt                                                                                                                 >  1:378183‑M2/119‑251 (MQ=255)
                                                                                                                                  cgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCt                                                                                                                  <  1:614712‑M2/132‑1 (MQ=255)
                                                                                                                                  cgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCt                                                                                                                  <  2:401949‑M2/132‑1 (MQ=255)
                                                                                                                                      atcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGtt                                                                                                              >  1:627512‑M2/115‑250 (MQ=255)
                                                                                                                                        caacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTAc                                                                                                            <  1:620520‑M2/138‑1 (MQ=255)
                                                                                                                                               gcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTTCAGAGATTg                                                                                                    <  2:109253‑M2/146‑1 (MQ=255)
                                                                                                                                                      ggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTg                                                                                             <  1:198794‑M2/153‑1 (MQ=255)
                                                                                                                                                       gcgattacggcaagagggcagcggagaaattgagatgcatccactaattgcatggatgagtaatgattaatgtggatagagtttctttgtgaggttgAAATCCATGTGTGAAGTGGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACTCTGGTTACAGAGATTGCATCCTGc                                                                                            <  2:133397‑M2/154‑1 (MQ=255)
                                                                                                                                                              cggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTccc                                                                                     >  1:96890‑M2/91‑251 (MQ=255)
                                                                                                                                                                gcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGc                                                                                   >  1:489173‑M2/89‑251 (MQ=255)
                                                                                                                                                                  aatagagctgcggagacattgagattcatccactacttgcatggatgagtaatgattaatgtggatagattttctttttgaggttgAAGTCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGACATTAACTGCCCAATTCAGGCGCCAACTGGTTTGATTGTACATTCCTTAAACGGAGGGGGTAAGCAAACACGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTc                                                                                 >  1:350251‑M2/87‑251 (MQ=255)
                                                                                                                                                                        gcagcggagaaattgagattcatccactacatgcatggatgagtaatgattaatttggatagagtttctttgtgaggttgAAAGCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTATCCGGAGTGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCttt                                                                           <  2:565497‑M2/171‑1 (MQ=255)
                                                                                                                                                                             ggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCgg                                                                      >  1:590174‑M2/76‑251 (MQ=255)
                                                                                                                                                                                 aaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACg                                                                                                                    <  2:224922‑M2/130‑1 (MQ=255)
                                                                                                                                                                                 aaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACg                                                                                                                    >  1:224922‑M2/72‑201 (MQ=255)
                                                                                                                                                                                 aaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGt                                                                  >  2:88326‑M2/72‑251 (MQ=255)
                                                                                                                                                                                  aattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTc                                                                 <  1:215315‑M2/181‑1 (MQ=255)
                                                                                                                                                                                            catccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCTCGCTGAttt                                                       <  1:297257‑M2/191‑1 (MQ=255)
                                                                                                                                                                                            catccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGAttt                                                       <  1:621847‑M2/191‑1 (MQ=255)
                                                                                                                                                                                             atccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGCtttt                                                      >  2:48656‑M2/60‑251 (MQ=255)
                                                                                                                                                                                               ccgctacttgcattgtgtagtaatgattaatgtggatagaatttcttttttaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACTTTCCTTAACCGGAGGGTGTAAGCATACCCGCTACGCTTGTTACAGAGATTTCATCCTGCAATTCCCGCTCCCCTTTTGCGTCCGTCGCGCTGAttttt                                                     <  2:619737‑M2/193‑1 (MQ=255)
                                                                                                                                                                                                 cctacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTg                                                  <  2:470933‑M2/196‑1 (MQ=255)
                                                                                                                                                                                                       tgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGttt                                            >  1:325675‑M2/50‑251 (MQ=255)
                                                                                                                                                                                                                gagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATg                                   >  1:357812‑M2/41‑251 (MQ=255)
                                                                                                                                                                                                                gagtaatgattaatgaggatagagcttctttttgcggctgAAATCAATGTGTGAAGGTGATCACAAATTTAACCACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTAAGGCGCCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGAGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATg                                   >  2:521626‑M2/41‑251 (MQ=255)
                                                                                                                                                                                                                    aatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTGGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGcc                               <  2:503909‑M2/215‑1 (MQ=255)
                                                                                                                                                                                                                                   tagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGaa                 >  1:248376‑M2/22‑250 (MQ=255)
                                                                                                                                                                                                                                           ctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGGCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCACTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCt        >  2:157583‑M2/14‑251 (MQ=255)
                                                                                                                                                                                                                                            tttttgaggttgAAAGCCATGTGTGAAGTTGATCACAAAGTTAAACACTGGTAGGGTAAAAAGTTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTGAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCt        <  2:192826‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                             ttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTg       <  1:88326‑M2/239‑1 (MQ=255)
                                                                                                                                                                                                                                              tttgaggttgAAATCCATGTGTGAAGTTGATAACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGACAACTGTTTTGATTGTCCATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGTTACCCTTTTGCGGCCGTCGCTCTGTTTTTTCTTGCGTTTGCGGAAATGGGCTAACTCTGCGAGGATAGCTATGCTGa      >  2:283998‑M2/11‑250 (MQ=255)
                                                                                                                                                                                                                                                 gaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCtgaaagg  >  2:37663‑M2/8‑249 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
ACCATCAGTGGTTTGCCTCTTCGACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAG  >  NC_000913/4330021‑4330512

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