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 50 [0] [1] 50 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)49 (0.750) 44/476 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

TACCATCAGTGGTTTGCCTCTTCGACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGAAAAA  >  NC_000913/4330020‑4330518
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
aaccttattccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgaa                                                                                                                                                                                                                                                          <  1:209317‑M2/2‑1 (MQ=255)
aaccttattccgcgcaactgtcggggcagcataacgacaccaacgtgctggcgtttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgaa                                                                                                                                                                                                                                                          <  2:195197‑M2/2‑1 (MQ=255)
       ttccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCAtg                                                                                                                                                                                                                                                   >  1:331165‑M2/243‑251 (MQ=255)
                  tttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTg                                                                                                                                                                                                                                         >  2:352155‑M2/232‑250 (MQ=255)
                    tcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAAGCCATGTGTGAAGTTGAt                                                                                                                                                                                                                                       <  2:202668‑M2/21‑1 (MQ=255)
                    tcgcggaagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGAt                                                                                                                                                                                                                                       <  1:222452‑M2/21‑1 (MQ=255)
                           agcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAAtt                                                                                                                                                                                                                               <  1:365865‑M2/29‑1 (MQ=255)
                            gcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAAttt                                                                                                                                                                                                                              <  1:352155‑M2/30‑1 (MQ=255)
                            gcataacgacaccaacgtgctggcttgtggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgatgaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAAttt                                                                                                                                                                                                                              <  2:316355‑M2/30‑1 (MQ=255)
                                      accaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgaggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGt                                                                                                                                                                                                                    >  1:57657‑M2/212‑251 (MQ=255)
                                               atggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggattcgtggctgggcgcacagaacgaatgcggtcgtcatcaacaacgcgtggaggcgattacgtcgatagagcagcggagaaattgagattcatccactacatgcctggatgagtcatgattaatgtggctagagtttctttttgaggttgAATTCCATGAGTGAAGTTGATCACAAATTTAATCACTGGGAGGGGaaa                                                                                                                                                                                                            >  2:559‑M2/203‑250 (MQ=255)
                                                         gttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatcccctacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAg                                                                                                                                                                                                                  >  2:15322‑M2/193‑234 (MQ=255)
                                                         gttcacgagtggttggactcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAg                                                                                                                                                                                                                  <  1:15322‑M2/42‑1 (MQ=255)
                                                                            cgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcagggatgagtaatgatcaatgtggatagattttctttttgagggtgAAATCCATGTGTGAAGTTGATCACAAATTTAAACCCTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGt                                                                                                                                                                                <  2:424638‑M2/76‑1 (MQ=255)
                                                                                   gcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGtt                                                                                                                                                                       >  1:80089‑M2/167‑251 (MQ=255)
                                                                                       aaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGAt                                                                                                                                                                   <  1:85318‑M2/89‑1 (MQ=255)
                                                                                                       gttgctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAAcc                                                                                                                                                   <  1:364884‑M2/105‑1 (MQ=255)
                                                                                                           ctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCgg                                                                                                                                                 >  1:309777‑M2/143‑249 (MQ=255)
                                                                                                                 ccacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgatgttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTa                                                                                                                                         <  2:103016‑M2/115‑1 (MQ=255)
                                                                                                                          gaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAAccc                                                                                                                                >  1:179655‑M2/128‑251 (MQ=255)
                                                                                                                             ggcggtcgtcatcaactacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcattgatgagtaatgattaatggggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTACTTAACCGGAGGGTGTAAGCAAACCCGCt                                                                                                                             <  2:57657‑M2/127‑1 (MQ=255)
                                                                                                                                           acaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgaggaatgattaatgcggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAATCACTGGTAGGGTAAAAAGGACATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGAAGGGTGTAAGAAAACCCGCTACACTTGTTAc                                                                                                                  >  2:377077‑M2/111‑248 (MQ=255)
                                                                                                                                            caacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACagag                                                                                                              <  2:179655‑M2/142‑1 (MQ=255)
                                                                                                                                             aacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCAt                                                                                                                                                                                                    <  1:335834‑M2/56‑1 (MQ=255)
                                                                                                                                             aacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCAt                                                                                                                                                                                                    >  2:335834‑M2/109‑164 (MQ=255)
                                                                                                                                                 cgtggaggcgattacggcaatagagcagcggagaaattgagattcatccacttcttgaatggatgagtaatgattaatgtgaatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTGTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGc                                                                                                         <  2:348027‑M2/147‑1 (MQ=255)
                                                                                                                                                  gtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGc                                                                                                         >  1:323653‑M2/104‑250 (MQ=255)
                                                                                                                                                    ggaggcgattacggaaatagagcagaggagaaattgagattcatccactacttgcatggatgagaaatgattcatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCAAAAATTTAAACACTGGGAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCAGCTAAGCTTGTTTAAGAGATTGCATc                                                                                                      >  2:461888‑M2/102‑251 (MQ=255)
                                                                                                                                                           attacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAAt                                                                                               >  2:466075‑M2/95‑251 (MQ=255)
                                                                                                                                                              acggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTcc                                                                                            <  1:417421‑M2/160‑1 (MQ=255)
                                                                                                                                                                      agagcagcggagaaattgagattcatccactacttgcatggattagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGGCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGGGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTcccc                                                                                    <  2:180104‑M2/168‑1 (MQ=255)
                                                                                                                                                                        agcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCtt                                                                                  >  2:410071‑M2/82‑251 (MQ=255)
                                                                                                                                                                               gagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGGTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGc                                                                           <  1:382040‑M2/177‑1 (MQ=255)
                                                                                                                                                                                agaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGcc                                                                          >  2:397440‑M2/74‑251 (MQ=255)
                                                                                                                                                                                 gaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCg                                                                         >  2:103786‑M2/73‑251 (MQ=255)
                                                                                                                                                                                    attgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTcg                                                                      >  1:377800‑M2/70‑251 (MQ=255)
                                                                                                                                                                                        agattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGATGATCACAAATTTAAACACTGGTAGGGTAAAAATGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTTCATTCCTTAACCGGAGGGTGTAAGCAAGCCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCt                                                                  <  1:61012‑M2/186‑1 (MQ=255)
                                                                                                                                                                                             catccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGAttt                                                             >  2:113090‑M2/61‑251 (MQ=255)
                                                                                                                                                                                                  actacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCTGCCGTCGCGCTGATTTTTCTg                                                        <  2:111059‑M2/196‑1 (MQ=255)
                                                                                                                                                                                                        tgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGtt                                                   <  2:36187‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                        tgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGt                                                    <  1:66542‑M2/200‑1 (MQ=255)
                                                                                                                                                                                                                   gtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAAAACTGGCAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGACAAATGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCACCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATggg                                       >  2:166918‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                                            aatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTg                              >  1:326449‑M2/30‑251 (MQ=255)
                                                                                                                                                                                                                             gtgtgggtagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACAGTGGTAGGGTAAACAGGTCATTAACTGCCCAATTCAGGAGTCAACTGGTTTGATTGTACATTCATTAACCGGAGGGTGTAAGAAAACCCTCTACGCTTGTTACAGCGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGGTTGCGGAAATGTGCCAACTCTGc                             >  2:428399‑M2/29‑251 (MQ=255)
                                                                                                                                                                                                                             atgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGc                             >  1:53988‑M2/29‑251 (MQ=255)
                                                                                                                                                                                                                                   atagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGaa                       <  1:128616‑M2/229‑1 (MQ=255)
                                                                                                                                                                                                                                           tctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATg                <  1:405958‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                                                     gttgAAATCCAGGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGGCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGAGTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCGACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGaa     <  2:194229‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGGCGCGCTGATTTTTCTGGCGTTTGCGGACATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGaa     >  2:493589‑M2/4‑250 (MQ=255)
                                                                                                                                                                                                                                                        gAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGaaaaa  <  1:178689/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
TACCATCAGTGGTTTGCCTCTTCGACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGAAAAA  >  NC_000913/4330020‑4330518

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