Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,085,042 Δ109,115 bp [tsuA][mrp] 106 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2085042 2194156 109115 49 [0] [0] 49 [tsuA]–[mrp] [tsuA],plaP,yoeI,yeeY,yeeZ,yoeB,yefM,hisL,hisG,hisD,hisC,hisB,hisH,hisA,hisF,hisI,wzzB,ugd,gnd,wbbL,insH7,wbbL,wbbK,wbbJ,wbbI,wbbH,glf,rfbX,rfbC,rfbA,rfbD,rfbB,galF,wcaM,wcaL,wcaK,wzxC,wcaJ,cpsG,cpsB,wcaI,gmm,fcl,gmd,wcaF,wcaE,wcaD,wcaC,wcaB,wcaA,wzc,wzb,wza,yegH,asmA,dcd,udk,dgcE,alkA,yegD,yegI,yegJ,pphC,yegL,sibA,ibsA,sibB,ibsB,mdtA,mdtB,mdtC,mdtD,baeS,baeR,yegP,trhP,cyaR,ogrK,yegR,yegS,gatR,insE5,insF5,gatR,gatD,gatB,gatA,gatZ,gatY,fbaB,yegT,yegU,yegV,yegW,yegX,thiD,thiM,rcnR,rcnA,rcnB,yehA,yehB,yehC,yehD,yehE,[mrp]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 20850410 (0.000)49 (0.880) 44/488 0.1 100% coding (485/1059 nt) tsuA thiosulfate transporter
?NC_000913 2194157 = 0 (0.000)coding (12/1110 nt) mrp P‑loop NTPase family protein Mrp

TAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTG  >  NC_000913/2193922‑2194403
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
agaaattgcagaagattggccggtggggaagtgagtcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgaggagtacagctgtgacaaagggatgccagcttttttcaaacagaatatgagcgagcccggttcggcgcggcggtaagggcgcgaggttgagttttgtcttcttcaactctttcagcaccacccagtgcgtgatgtTTGTTCTTTCGTCCCg                                                                                                                                                                                                                                         <  2:412219‑M2/16‑1 (MQ=255)
  aaattgcagaatattggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaagggcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGAc                                                                                                                                                                                                                                       <  1:265749‑M2/18‑1 (MQ=255)
                         gggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCa                                                                                                                                                                                                                <  1:267789‑M2/41‑1 (MQ=255)
                                       aacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGttt                                                                                                                                                                                                  <  1:12699‑M2/55‑1 (MQ=255)
                                             aacattcgtccggttgtttcactcagggtccaggctaaaagcgcaattataccgatgagtacagctgtgacaaagggatgcgaacgtttttcaaacagaataggagcgttcccggttctgcgcggcggtaaggtcgcgactttgagttttggcttcttcaactcttgcatcaccccccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTAACTga                                                                                                                                                                                            <  2:240865‑M2/61‑1 (MQ=255)
                                                  gcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTgaagaa                                                                                                                                                                                        <  2:136291‑M2/65‑1 (MQ=255)
                                                  gcgtccggtggcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGAGAAACCATCAGCTGATAGTTTACCTGAAGAat                                                                                                                                                                                       <  2:56441‑M2/66‑1 (MQ=255)
                                                    gtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAatat                                                                                                                                                                                     <  1:273823‑M2/68‑1 (MQ=255)
                                                        ggttgcgtcactcaggggccaggctgaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttgtcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCTACTCCTTATTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATagag                                                                                                                                                                                 <  2:148487‑M2/72‑1 (MQ=255)
                                                            gcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGt                                                                                                                                                                             <  1:368072‑M2/75‑1 (MQ=255)
                                                                     aggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTActtactta                                                                                                                                                                    <  1:329942‑M2/85‑1 (MQ=255)
                                                                                 aaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccccccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCAt                                                                                                                                                        >  2:175553‑M2/155‑251 (MQ=255)
                                                                                    agcgcaatcaaaccgatgagtacagctgtgacaacgggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatacTTGTTCGTTCATACCGACTCCTTTTTTGTATCGATAAACCATCAGCTGATAGTTTAACTGACGAATATAGAGAAGTACTTACTTAACATTTTCCCatttt                                                                                                                                                     >  2:304290‑M2/152‑250 (MQ=255)
                                                                                         aatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGCCTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGt                                                                                                                                                                             >  2:176038‑M2/147‑222 (MQ=255)
                                                                                         aatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGt                                                                                                                                                                             <  1:176038‑M2/76‑1 (MQ=255)
                                                                                                     gagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccaccccgagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTTTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATATAACCCCtt                                                                                                                                    >  2:243759‑M2/135‑251 (MQ=255)
                                                                                                         acagctgtgacaaagggatgccagcgtttttcaaacagaatctgcgcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGTTAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTg                                                                                                                                                     >  1:32651‑M2/131‑230 (MQ=255)
                                                                                                         acagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTg                                                                                                                                                     <  2:32651‑M2/100‑1 (MQ=255)
                                                                                                         acagctgtgacaaagggatgccagcggttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTc                                                                                                                                 <  2:311847‑M2/120‑1 (MQ=255)
                                                                                                               gtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTAtt                                                                                                                           <  1:137069‑M2/126‑1 (MQ=255)
                                                                                                               gtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTAtt                                                                                                                           >  1:277688‑M2/125‑250 (MQ=255)
                                                                                                                  acaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCCTCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTa                                                                                                                          >  2:423701‑M2/122‑248 (MQ=255)
                                                                                                                   caaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCAATATTaaga                                                                                                                       >  2:300644‑M2/121‑250 (MQ=255)
                                                                                                                    caagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTaaga                                                                                                                       >  2:153122‑M2/120‑249 (MQ=255)
                                                                                                                     aagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGt                                                                                                                    >  1:91282‑M2/119‑251 (MQ=255)
                                                                                                                        ggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAAt                                                                                                                 <  1:327597‑M2/136‑1 (MQ=255)
                                                                                                                          atgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatctccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGc                                                                                                               <  1:312248‑M2/138‑1 (MQ=255)
                                                                                                                              cagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCctac                                                                                                           <  1:175553‑M2/142‑1 (MQ=255)
                                                                                                                                 cgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTAt                                                                                                        <  1:7327‑M2/145‑1 (MQ=255)
                                                                                                                                     tttcaaacagaatatgagagatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGGCTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACt                                                                                                    >  1:361662‑M2/103‑251 (MQ=255)
                                                                                                                                             agaatctgagcgatcacggttcggcgcggcggtaaggtcacgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCCTCCCGACTCCTTTTTTGCATAGCTTAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTAATTACTTAACATTTTCCCATTTCGTACTATCTAACCCCTTTTCACTATTAAAAAGTATTGCCTACTATGACTCAAGAcgc                                                                                            >  2:376886‑M2/95‑251 (MQ=255)
                                                                                                                                             agaatatgagcgatcccggttcggcgcggcggtaaggtcgcgaatttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATCGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCCAGTcgc                                                                                            >  2:52998‑M2/95‑251 (MQ=255)
                                                                                                                                                  atgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCga                                                                                          >  2:221646‑M2/90‑248 (MQ=255)
                                                                                                                                                            ccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactatttcatcaccaccccgagcgctatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATACACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGt                                                                             >  2:428891‑M2/80‑251 (MQ=255)
                                                                                                                                                                        ggcggtaaggtcgcgactttgagttttggcttcttcaacgctttcatcaccacccagagcgttatcaTTGTTCGTTCATACCGACTCCTGTTTTGTATAGATAAGCCAGCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTGCTTACGTAGCATTTTCCCATTTGGTACTATCGAACCCCTTTTCACTATTAATAAGTAATGCCTACTAGGACGCAAGtcgctaagaaaattctattgacggtcgcact                                                                 <  2:170128‑M2/184‑16 (MQ=255)
                                                                                                                                                                          cggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGc                                                               <  2:100066‑M2/186‑1 (MQ=255)
                                                                                                                                                                                      gactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACg                                                   >  1:348051‑M2/54‑251 (MQ=255)
                                                                                                                                                                                           tgagttttggcttcttcaactctttcatcacaacccagagcgttattaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACACCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCCCTGCCGTACGCTAACGGCTc                                               >  2:324082‑M2/49‑250 (MQ=255)
                                                                                                                                                                                             agttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCTAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAAt                                            <  1:22499‑M2/205‑1 (MQ=255)
                                                                                                                                                                                                   ggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCt                                      <  1:77356‑M2/211‑1 (MQ=255)
                                                                                                                                                                                                     cttcttcaactctttcatcacaacccagcgcgttctcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGACGTAATGCCTACTATGCCTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAGCGGCTAAATCCACATCg                                    >  2:237565‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                      ttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCgg                                   >  1:132414‑M2/38‑251 (MQ=255)
                                                                                                                                                                                                         ttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGcc                                 >  1:411173‑M2/35‑250 (MQ=255)
                                                                                                                                                                                                             actctttcctcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATg                            >  2:147474‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                             actctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATg                            >  2:114547‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                               tctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCt                          >  1:195611‑M2/29‑251 (MQ=255)
                                                                                                                                                                                                                ctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTg                         >  2:57502‑M2/28‑251 (MQ=255)
                                                                                                                                                                                                                                       atcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTGATGCCTTCTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTTCGCACTGCCGTACGCTAACGGCGCAATCCACCTCGGCCATATGCGGGAGCACATCCAGGCTGATGTCt   <  1:259992‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                       atcaTTGTTCGTTAATCCCGACTCCTTTTTTGTATAGATAAAAAATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACCTTTTCCCATTTGGCACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTg  >  1:284999‑M2/5‑251 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
TAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTG  >  NC_000913/2193922‑2194403

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: