Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,196,474 Δ29,182 bp [yehF]yohC 21 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2196474 2225655 29182 45 [0] [0] 44 [yehF]–yohC [yehF],yehI,yehK,yehL,yehM,yehP,yehQ,yehR,yehS,btsR,btsS,mlrA,yohO,yehW,yehX,yehY,osmF,bglX,dld,pbpG,yohC

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 21964730 (0.000)44 (0.600) 41/492 0.3 100% intergenic (+140/‑1) metG/yehF methionine‑‑tRNA ligase/PF05406 family protein YehF
?NC_000913 2225656 = 0 (0.000)intergenic (‑25/‑145) yohC/yohD putative inner membrane protein/DedA family protein YohD

AAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCGGAAGGTGAAACCGTGACTTTGCTGGGAGGCGTTGCG  >  NC_000913/2225407‑2225899
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
gcgcttcggtatctctgaaggcatggtgatggctgccggtcctggcgggaaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctcgctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaga                                                                                                                                                                                                                                                    <  1:451835‑M2/2‑1 (MQ=255)
      cggtatctctgaaggcatggtgatggctgcnggtcctggcgggaaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTg                                                                                                                                                                                                                                               <  2:293925‑M2/7‑1 (MQ=255)
          atctctgaaggcatggtgatggctgccggtcctggcgggaaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaagtcgctgcatcgacagcgccttttctttataaattcctaaagttgtgttcttgcgatgttgtctctctctaacccgcataaatactggtagcatctgcatgcaactggataaaattacagggatgcagaGAACCTGACAg                                                                                                                                                                                                                                           <  2:234879‑M2/11‑1 (MQ=255)
                  aggcatggtgatggctgccggtcctggcggtaaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGtat                                                                                                                                                                                                                                  <  2:30585‑M2/20‑1 (MQ=255)
                       gggttatggctgccggtcctggcgggaaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatccccctccaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactgggagcatctgcattcaactgtataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTc                                                                                                                                                                                                                             <  2:174860‑M2/25‑1 (MQ=255)
                              ggctgccggtcctggcgggaaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTtgt                                                                                                                                                                                                                      <  1:547131‑M2/32‑1 (MQ=255)
                                           ggtgggcaggatgttttcctgctaagcccggatgccggtgctaatccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaatttgttttcttgcgattgtgtctctctctaaccctcagaaatactggtagcatctgcattcaactggataaaataacagggattcagaGAAGCTGACAGTTCAAGTATAAGTCAGCTTGGGATTAGTTTgtgt                                                                                                                                                                                                         <  2:531695‑M2/45‑1 (MQ=255)
                                              gggaaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattactgggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCg                                                                                                                                                                                                      <  2:63861‑M2/48‑1 (MQ=255)
                                                 aaagatattttcctgctaagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGccc                                                                                                                                                                                                   <  1:259455‑M2/51‑1 (MQ=255)
                                                                  aagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTc                                                                                                                                                                                  >  1:56636‑M2/184‑251 (MQ=255)
                                                                  aagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTc                                                                                                                                                                                  <  1:498583‑M2/68‑1 (MQ=255)
                                                                  aagcccggatgccggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTc                                                                                                                                                                                  >  1:476372‑M2/184‑251 (MQ=255)
                                                                              cggtgctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTcc                                                                                                                                                                      <  1:457816‑M2/80‑1 (MQ=255)
                                                                                  gctaaaccgggtcatcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgactctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCttt                                                                                                                                                                   <  1:48617‑M2/83‑1 (MQ=255)
                                                                                                tcaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaaatacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGAt                                                                                                                                                    >  1:538341‑M2/154‑251 (MQ=255)
                                                                                                acaggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCGTGTGATTATTTTTTGTTCGCCGATGTAATTGTCACTGTCTGATGAACGTCCTGTTTCAGCGTGAATGAt                                                                                                                                                    <  2:444053‑M2/98‑1 (MQ=255)
                                                                                                 caggtgaaataatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcatcaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGAtt                                                                                                                                                   >  2:12347‑M2/153‑251 (MQ=255)
                                                                                                 caggccacataatcccccttcaaggcgatgcctcgacagcgccttttctttctaaattcctcaagttgttttcttgcgatcttggctctctagaacccgcataaatactggaagcatctgcattcaactggctaaatttacagggatgcagaGAACATGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTGCGCCCATGTCATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTTCATGAtt                                                                                                                                                   >  2:202962‑M2/153‑251 (MQ=255)
                                                                                                           aatcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATTTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGgtgt                                                                                                                                         <  1:63468‑M2/109‑1 (MQ=255)
                                                                                                            atcccccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatatgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGgtgt                                                                                                                                         >  1:244373‑M2/142‑250 (MQ=255)
                                                                                                                cccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCt                                                                                                                                    <  1:411577‑M2/114‑1 (MQ=255)
                                                                                                                cccttcaaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCt                                                                                                                                    >  1:224383‑M2/138‑251 (MQ=255)
                                                                                                                    tcaaggctctgcttcggtagcgccttgtctttataaattccgaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaacgtgataaaattacagggatgcagaTAACCTTACAGTTCAAGTATAAGTCAGCTTGTTATTATTCTTTGTTCGCCCATGTAATTTTCACTGTCTGTCGAACGTCCTTTTGCAGCGTGAATGATTCACAGTGTGTATGCTTa                                                                                                                                  <  2:320507‑M2/116‑1 (MQ=255)
                                                                                                                      aaggcgctgcatcgacagcgccttttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTgg                                                                                                                              >  2:6858‑M2/132‑251 (MQ=255)
                                                                                                                                   gacagctccttttctttataaattcctaaagttggtttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGgag                                                                                                                 <  2:404644‑M2/133‑1 (MQ=255)
                                                                                                                                             tttctttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCACAGGGAGACAGACTGGc                                                                                                       >  2:527378‑M2/109‑251 (MQ=255)
                                                                                                                                                  ttataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACATACTGGCCTAt                                                                                                   <  1:423204‑M2/147‑1 (MQ=255)
                                                                                                                                                    ataaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTGTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGa                                                                                                <  1:262533‑M2/150‑1 (MQ=255)
                                                                                                                                                      aaattcctaaagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAGCCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGGCCTTTTTCAGCGTGAATGATTCTCAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGAtc                                                                                              <  1:197902‑M2/152‑1 (MQ=255)
                                                                                                                                                           cctaaagttgtttgcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaTAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGATCGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCaa                                                                                          <  1:332457‑M2/156‑1 (MQ=255)
                                                                                                                                                               aagttgttttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATacac                                                                                     <  2:381354‑M2/161‑1 (MQ=255)
                                                                                                                                                                      tttcttgcgattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACATTTCAAGTATAAGTCAGCTTGTTATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTAtct                                                                               <  1:413911‑M2/167‑1 (MQ=255)
                                                                                                                                                                              gattttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATg                                                                      >  2:138115‑M2/76‑251 (MQ=255)
                                                                                                                                                                                 tttgtctctctctcacccgcataaatactggtagcatctgcattcacctggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCATATGGCTCTCAATACACTTATCTCCCAATAAGGtt                                                                   >  2:454521‑M2/73‑251 (MQ=255)
                                                                                                                                                                                 tttgtctctctctaacccgcataaatactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGtt                                                                   >  1:233453‑M2/73‑251 (MQ=255)
                                                                                                                                                                                                         atactggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGcc                                           >  1:573620‑M2/49‑251 (MQ=255)
                                                                                                                                                                                                           actggtagcatctgcattaaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTg                                         >  1:161234‑M2/47‑251 (MQ=255)
                                                                                                                                                                                                              ggtagcatctgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCg                                      <  1:380349‑M2/208‑1 (MQ=255)
                                                                                                                                                                                                                  gcatctgcattcaactggataaaattacaggnatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATAGGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCGGAAg                                  <  2:273797‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                                       tgcattcaactggataaaattacagggatgcagaGAACCTGACAGTTCAAGTCTAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTAACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCCAAGAGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCGGAAGGTGa                              >  2:508590‑M2/35‑250 (MQ=255)
                                                                                                                                                                                                                                      aaaattacagggatgcagaTAACCTGAAAGTTCAATTATACGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCGCAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCGGAAGGTGAAACCGTGACTTTGCTg              <  2:563807‑M2/232‑1 (MQ=255)
                                                                                                                                                                                                                                           tacagggatgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCGGAAGGTGAAACCGTGACTTTGCTTGGAgg         >  2:77241‑M2/15‑251 (MQ=255)
                                                                                                                                                                                                                                            acagggaggcagaCAACCTGACAGTTCAAGGATAAGTCAGATTGTGATTATTTTTTGTTCGTCCATGTAATGTTCACTGCCTGATGAACGTCATTGTTCAGCGTGATTGATTCACAGGGTGTATGCTTACTGGCAACCATAGGGAGACAGACTGGCCTATGGATCTCAATACTCTTATCGCACAATATGGTTATGGCGCGCTGGGGAGCGGTAGCCTGGCGGAAGGTGAAACCGTGACTTTGCTGGGAgg         <  2:538341‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                                  atgcagaGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCGGAAGGTGAAACCGTGACTTTGCTGGGAGGCGTTgcg  >  1:224115‑M2/8‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   
AAACAAGACAGCCAGCGCCAGTCCGGTAAACCAGGATAACTTCAGGATAGTGCCATCGCCAAAATTCCAGCCAATCTGTGTAGTGCCAATGAAGGCGCAAATCACCGGAATCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCCCATGTAATTTTCACTGTCTGATGAACGTCCTTTTTCAGCGTGAATGATTCACAGGGTGTATGCTTACTGGCAACCAAAGGGAGACAGACTGGCCTATGGATCTCAATACACTTATCTCACAATATGGTTATGCCGCGCTGGTGATCGGTAGCCTGGCGGAAGGTGAAACCGTGACTTTGCTGGGAGGCGTTGCG  >  NC_000913/2225407‑2225899

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