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 42 [0] [1] 43 [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)42 (0.750) 39/486 0.3 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

TTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTGGGT  >  NC_000913/2193911‑2194406
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
cacccgcgaccagaaattgcagaatgttggccggtggggaagtggttactaacccaaacatgcgtccggttgcttcacgcaggggccaggctaaaagggcaatcgaaccgatgagtacagctgtgacaaagggatgccagcgttttgcaaacagaatatgagcgatcccggttcggcgcggcggtaagggcgcgactttgagttttggcttctgcaactctttcatcaccgcccagagcgttatcaTTgtt                                                                                                                                                                                                                                                       <  2:450104‑M2/5‑1 (MQ=255)
      cgaccagaaattgcagaatattggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtacggtcgcgactttgagttttggcttcttcaactctttcatcccccccccgagcgttatcaTTgttcgttc                                                                                                                                                                                                                                                  >  2:119138‑M2/241‑250 (MQ=255)
                     gaatattggccgttgcggaagcaattcctcacccacactagcgtcaggttccttcactcaggggccaggctaaaagcgcaatccaaccgatgagtccagctgtgacacagggatgccagcgtttttaaaacagaatatgagcgatcccggtacggcgcggcggtaaggtcgcgtctttgagttttggcttcttcaactctttcctcaccacccagagcgttatcaTTGTTCGTTCCTCCCGACTCCttttt                                                                                                                                                                                                                                  >  2:81527‑M2/226‑251 (MQ=255)
                            ggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAg                                                                                                                                                                                                                           <  1:380967‑M2/33‑1 (MQ=255)
                               cgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgctgagtacagctgtgacaaagggatgccagcgtttttcacacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccaccccgagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGAt                                                                                                                                                                                                                         >  1:432449‑M2/216‑250 (MQ=255)
                                gttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatccTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATa                                                                                                                                                                                                                        >  2:72061‑M2/215‑250 (MQ=255)
                                   ggggaagtgattcctaacccaaacatgcgtccggttgcgtcactcagggggcaggctaaaagcgcaatcaatccgatgagtacagctgtgacaaagtgatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaagtcttgcatcaccacccatatcgttatcaTTGTTCGTTCATCCCGACTCCTTGTTTGTATAGATAAAc                                                                                                                                                                                                                     <  1:431967‑M2/39‑1 (MQ=255)
                                      gaagtgattcctcacccaaccatgcgtacggatgcttcactaaggggcccggacaaaagcgcaatcaaaccgctgagtacaactgtgccaaagggatgccagcgtttttcaaacagcatatgagcgatcccggttcggcgcggcggtaaggccgcgactttgagttttggcttcttccactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATc                                                                                                                                                                                                                 >  2:467472‑M2/209‑251 (MQ=255)
                                           gattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCt                                                                                                                                                                                                             <  2:62722‑M2/47‑1 (MQ=255)
                                               cctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatctccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAg                                                                                                                                                                                                        <  1:217981‑M2/52‑1 (MQ=255)
                                                                   ggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATag                                                                                                                                                                                      >  1:391645‑M2/180‑249 (MQ=255)
                                                                        cttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGt                                                                                                                                                                                >  1:479030‑M2/175‑250 (MQ=255)
                                                                                  gggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAAc                                                                                                                                                                     <  2:133817‑M2/87‑1 (MQ=255)
                                                                                     ccatgctaaaagcgcaatcaaaccgatgaggacagctgtgtcaaagggatgccagcttttttcaaacagaagatgagcgatccagtttcggctcggctgtaatgtcgcgattttgagttttggcttcttcacctctttcatcaccgcccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACGTGACAAATATAGATAAGTACTTACTTAACAtt                                                                                                                                                                  <  1:305136‑M2/90‑1 (MQ=255)
                                                                                     ccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATAGAGAGAAGTACTTACTTAACAtt                                                                                                                                                                  <  2:329724‑M2/90‑1 (MQ=255)
                                                                                               agcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTg                                                                                                                                                        >  1:450385‑M2/152‑251 (MQ=255)
                                                                                               agcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTGACTTAACATTTTCCCATTTg                                                                                                                                                        <  1:119138‑M2/100‑1 (MQ=255)
                                                                                                   caatccaaccgatgcgtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggacgcgactttgggttttggcttcttcaactctttcatcacaaccctgagcgttatccTTGTTCCTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGCATATAGACAAGTACTTTCTTAACCTTTTCCCATTTTGTAc                                                                                                                                                    >  2:223105‑M2/148‑251 (MQ=255)
                                                                                                             gatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAAccc                                                                                                                                          >  1:132168‑M2/138‑251 (MQ=255)
                                                                                                                gagtacagctgtgacaaagggatgccagcgtttttcaagcagactatgagcgatcccggttcggcgcagcggtaaggtcgcgcctttgagatttggcttcttcaactctttcatccccacccagagcgttctcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACACCtt                                                                                                                                       >  1:79148‑M2/135‑251 (MQ=255)
                                                                                                                  gtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCtttt                                                                                                                                     >  1:305616‑M2/133‑251 (MQ=255)
                                                                                                                          gtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTAtt                                                                                                                              <  1:166389‑M2/126‑1 (MQ=255)
                                                                                                                                        ccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccaccccgagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATCGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCct                                                                                                                >  2:425403‑M2/111‑250 (MQ=255)
                                                                                                                                         cagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTGTGGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTGGGTGCTATCTAACCCCTTTTCACTATTAAGAAGTAATGCctac                                                                                                              <  1:290487‑M2/142‑1 (MQ=255)
                                                                                                                                                          aatatgcgcgatccaggttcggcgcggcggtaaggtcgcgactttgcgttttggcttcttcaactatttcatccccacccagagcgttatccTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGCTAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTACCATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAATcgcg                                                                                              >  2:431573‑M2/93‑250 (MQ=255)
                                                                                                                                                            tatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCga                                                                                             >  2:19135‑M2/91‑249 (MQ=255)
                                                                                                                                                               gagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAATAAGTAATGCCTACTATGACCCAAGTCGCgaagaa                                                                                         >  2:468‑M2/88‑250 (MQ=255)
                                                                                                                                                                    atcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCt                                                                                   >  1:275444‑M2/83‑251 (MQ=255)
                                                                                                                                                                                cgcggcggtaaggtcgcgactttgaggtttggcttcttcaactctttcctcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCCAAGAAAATTCTGGTGACGTgcgc                                                                       >  2:221378‑M2/71‑251 (MQ=255)
                                                                                                                                                                                   ggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGAc                                                                                                        <  2:3700‑M2/148‑1 (MQ=255)
                                                                                                                                                                                   ggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGAc                                                                                                        >  1:3700‑M2/68‑215 (MQ=255)
                                                                                                                                                                                                   ctttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGc                                                    >  1:383459‑M2/52‑251 (MQ=255)
                                                                                                                                                                                                                cttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCg                                       >  1:280480‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                                  tcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCg                                       <  1:420212‑M2/213‑1 (MQ=255)
                                                                                                                                                                                                                       aactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAAGTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAAGCCACCTCGGCCatat                                <  2:56035‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                                       aactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCAGCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCatat                                <  2:150057‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                                                  caccacccagagcgttatcaTTGTTCTTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGcac                      >  1:43526‑M2/21‑250 (MQ=255)
                                                                                                                                                                                                                                          agagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCgaag                                                                                           >  2:14051‑M2/13‑173 (MQ=255)
                                                                                                                                                                                                                                          agagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCgaag                                                                                           <  1:14051‑M2/161‑1 (MQ=255)
                                                                                                                                                                                                                                          agagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCCCATCCAGGCt             >  1:274194‑M2/13‑251 (MQ=255)
                                                                                                                                                                                                                                                 tatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTCCCTGCAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCt      >  1:233899‑M2/6‑251 (MQ=255)
                                                                                                                                                                                                                                                    caTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTgg    <  1:72061‑M2/248‑1 (MQ=255)
                                                                                                                                                                                                                                                     aTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTACGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACGTCCAGGCTGATGTCTGGGt  >  1:131490/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
TTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTGGGT  >  NC_000913/2193911‑2194406

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: