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 60 [0] [0] 60 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)59 (0.750) 54/492 0.2 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

TACCATCAGTGGTTTGCCTCTTCGACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGAAAAAA  >  NC_000913/4330020‑4330519
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
aaccttatgccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgaa                                                                                                                                                                                                                                                           <  1:395701‑M2/2‑1 (MQ=255)
     tattccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATcc                                                                                                                                                                                                                                                       >  1:368534‑M2/245‑250 (MQ=255)
      attccgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgaggaatgattaatgtggatagagtttctttttgaggttgAAATCCAt                                                                                                                                                                                                                                                     >  2:71166‑M2/244‑251 (MQ=255)
         ccgcgcaactttcgcggcagcataacgacaccaacgtgctggctgttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatcctcttcgtgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCAtgt                                                                                                                                                                                                                                                   <  2:490599‑M2/10‑1 (MQ=255)
          cgcgcaactttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCAtgtgt                                                                                                                                                                                                                                                 <  2:341372‑M2/12‑1 (MQ=255)
              caactttcgcggaagcataacgacatcatcgtgctgtcttttggttcacgagtggttggcctcgaactggtaaaaatgattgtggatgcgtggctgggcgcaaagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGa                                                                                                                                                                                                                                               <  2:554213‑M2/14‑1 (MQ=255)
                   ttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTg                                                                                                                                                                                                                                          <  1:544908‑M2/19‑1 (MQ=255)
                   ttcgcggcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTg                                                                                                                                                                                                                                          <  1:316431‑M2/19‑1 (MQ=255)
                    gcgcggcagcataacgacaccaacgtgctggcttttggttcaagagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggctgtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgatgaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGAt                                                                                                                                                                                                                                        <  2:105590‑M2/21‑1 (MQ=255)
                         gcagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACaa                                                                                                                                                                                                                                   >  1:510217‑M2/225‑250 (MQ=255)
                          cagcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAAt                                                                                                                                                                                                                                 <  1:539469‑M2/28‑1 (MQ=255)
                           agcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAAtt                                                                                                                                                                                                                                <  1:87483‑M2/29‑1 (MQ=255)
                           agcataacgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAAtt                                                                                                                                                                                                                                >  1:216042‑M2/223‑251 (MQ=255)
                                  cgacaccaacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAacac                                                                                                                                                                                                                         <  2:629942‑M2/36‑1 (MQ=255)
                                        caacgtgctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGt                                                                                                                                                                                                                     >  1:32375‑M2/210‑249 (MQ=255)
                                               ctggcttttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTaaa                                                                                                                                                                                                             >  1:491578‑M2/203‑250 (MQ=255)
                                                     tttggttcacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTc                                                                                                                                                                                                       <  1:309070‑M2/54‑1 (MQ=255)
                                                           tcacgagtggttgggctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAAc                                                                                                                                                                                                 <  2:606090‑M2/60‑1 (MQ=255)
                                                            cacgagtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacaaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACt                                                                                                                                                                                                <  1:34890‑M2/61‑1 (MQ=255)
                                                                 gtggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCaa                                                                                                                                                                                          <  1:94836‑M2/67‑1 (MQ=255)
                                                                   ggttggcctcgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAAtt                                                                                                                                                                                        <  2:374531‑M2/69‑1 (MQ=255)
                                                                            cgaactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcagtagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTc                                                                                                                                                                                >  1:332317‑M2/174‑250 (MQ=255)
                                                                              aactggcaaaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattccggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagcgtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAACCACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAAc                                                                                                                                                                             >  2:318103‑M2/172‑251 (MQ=255)
                                                                                      aaaatgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTg                                                                                                                                                                      >  1:553550‑M2/164‑250 (MQ=255)
                                                                                         atgattgtggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGGTTGATTg                                                                                                                                                                  <  1:47868‑M2/91‑1 (MQ=255)
                                                                                                tggatgcgtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTc                                                                                                                                                           >  1:153240‑M2/154‑251 (MQ=255)
                                                                                                   atgcgtggctgggcgcacagtacgcaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCtt                                                                                                                                                        >  2:528642‑M2/151‑251 (MQ=255)
                                                                                                       gtggctgggcgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAAcc                                                                                                                                                    >  1:449890‑M2/147‑251 (MQ=255)
                                                                                                                cgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGgtgt                                                                                                                                           <  2:332857‑M2/114‑1 (MQ=255)
                                                                                                                cgcacagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGgtgt                                                                                                                                           >  1:509955‑M2/138‑251 (MQ=255)
                                                                                                                   acagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTa                                                                                                                                          >  1:130016‑M2/135‑249 (MQ=255)
                                                                                                                   acagtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAg                                                                                                                                        <  1:259020‑M2/117‑1 (MQ=255)
                                                                                                                     agtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGc                                                                                                                                       >  1:578606‑M2/133‑250 (MQ=255)
                                                                                                                      gtacgaaggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCaa                                                                                                                                     <  1:117698‑M2/120‑1 (MQ=255)
                                                                                                                        acgaaggcggtcgtcatcagcaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAAc                                                                                                                                   <  2:190940‑M2/122‑1 (MQ=255)
                                                                                                                           aaggcggttgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTCGGGTAAAAAGGTAATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGAAAACCCg                                                                                                                                >  2:280198‑M2/127‑251 (MQ=255)
                                                                                                                            aggcggtcgtcatcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGc                                                                                                                               <  1:321339‑M2/126‑1 (MQ=255)
                                                                                                                                       atcaacaacgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCCGGCGTCAACTGGTTTGATTGTTCATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACCCttgat                                                                                                                     >  2:425442‑M2/115‑248 (MQ=255)
                                                                                                                                         caacagcgcgtggaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACa                                                                                                                  <  1:189513‑M2/139‑1 (MQ=255)
                                                                                                                                                     gaggcgattacggcaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCGACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGATATTGCATcc                                                                                                      <  2:244428‑M2/151‑1 (MQ=255)
                                                                                                                                                            ttacgacaatagagcagcggagaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttatttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAATACTGTTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACACGCTACGCTTGTTACAGAGATTGCATCCTGCAAtt                                                                                               >  2:458386‑M2/94‑251 (MQ=255)
                                                                                                                                                                                agaaattgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGcc                                                                           >  1:201926‑M2/74‑251 (MQ=255)
                                                                                                                                                                                    attgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTc                                                                                                                                                                                                       <  2:247632‑M2/54‑1 (MQ=255)
                                                                                                                                                                                    attgagattcatccactacttgcatggatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTc                                                                                                                                                                                                       >  1:247632‑M2/70‑123 (MQ=255)
                                                                                                                                                                                                         gcatggatgcgtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGGAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTg                                                  >  2:430000‑M2/49‑251 (MQ=255)
                                                                                                                                                                                                              gatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGa                                              >  1:619038‑M2/44‑250 (MQ=255)
                                                                                                                                                                                                              gatgagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGa                                              >  1:448145‑M2/44‑250 (MQ=255)
                                                                                                                                                                                                                 gagtaatgattaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGATTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATg                                          >  1:17604‑M2/41‑251 (MQ=255)
                                                                                                                                                                                                                        gcttaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAAc                                   >  2:260795‑M2/34‑251 (MQ=255)
                                                                                                                                                                                                                          ttaatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAActc                                 >  2:174522‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                                            aatgtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTg                               <  2:263855‑M2/222‑1 (MQ=255)
                                                                                                                                                                                                                               gtggatagagtttctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCg                             >  2:408615‑M2/27‑250 (MQ=255)
                                                                                                                                                                                                                                           tctttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATg                 <  1:493306‑M2/236‑1 (MQ=255)
                                                                                                                                                                                                                                             tttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTg              >  1:509318‑M2/13‑251 (MQ=255)
                                                                                                                                                                                                                                              ttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCATCTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGc                <  2:592876‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                              ttttgaggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGc                <  1:318103‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgaggttgAAATCCATGTGTGAAGTTGATCACAAAGTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGTCGGCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACTCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGACAACTCTGCGAGGAAAGCTATGCTGaaa           <  2:329193‑M2/242‑1 (MQ=255)
                                                                                                                                                                                                                                                   aggttgAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAgg        <  1:543917‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                                     gttgAAATCCATGGATGAAGTTGATCACAAATTTAAACACTGGTAGTGTAAAAAGGTCGTTAACTGCCCAATTCAGGGGTCAACGGGTTTGATTGTACATTCCGTAGCCGGAGGGTGTAAGCAGACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTGGGCGGCCGTCGCGGTGATGTTTCTGGCGATTGCGGTAATGGGCCAACTCTGCGAGGAAAGGTATGCTGAAAAGGaa      <  2:430322‑M2/247‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaaTCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAATAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAACCCCGCTACGCTTGTTACAGAGATTTCATCCTTCAATTCCCGCTCCCCTTTCTCGGCCCTCGAGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAACAGGaaaaaa  <  2:401362/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
TACCATCAGTGGTTTGCCTCTTCGACCGTTCCGACGGGCGGGGAGCATGCTGATTTTTTAAGTTGGCTGGGAAAGTGGACCACGGCAGATAAACAACCCGTTTGCTGGTCAGTGACCCAACGCTGGCAAACTGTCGCGCTGGGGATGCCACGACTCTGTTCAGCGCAGCGTCTGGCGGGGGCAATGCTCGAGGAAATCTTCTCTGTAAATTTGGCGTAAATAATCAGTTACATCAATGAGTCCTAAACGAAATCCATGTGTGAAGTTGATCACAAATTTAAACACTGGTAGGGTAAAAAGGTCATTAACTGCCCAATTCAGGCGTCAACTGGTTTGATTGTACATTCCTTAACCGGAGGGTGTAAGCAAACCCGCTACGCTTGTTACAGAGATTGCATCCTGCAATTCCCGCTCCCCTTTTGCGGCCGTCGCGCTGATTTTTCTGGCGTTTGCGGAAATGGGCCAACTCTGCGAGGAAAGCTATGCTGAAAAGGAAAAAA  >  NC_000913/4330020‑4330519

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