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 47 [0] [0] 47 [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)47 (0.660) 43/488 0.4 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

TCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTGG  >  NC_000913/2193918‑2194404
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
gaccagaaattgcagaatattggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgcgttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCAt                                                                                                                                                                                                                                              >  1:367868‑M2/240‑251 (MQ=255)
  ccagaaattgcagaatattggccgttggggcagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaagggcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttttcaTTGTTCGTTCATcc                                                                                                                                                                                                                                            >  2:423461‑M2/238‑251 (MQ=255)
  ccagaaattgcagaatattggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATcc                                                                                                                                                                                                                                            >  1:161057‑M2/238‑251 (MQ=255)
    agaaattgcagaatattggccgttggggaagtgattcctaacccaaacattcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcggttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCg                                                                                                                                                                                                                                          <  2:197222‑M2/16‑1 (MQ=255)
      aaattgcagaatattggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcttttttcaaacagaatatgagcgagcccggttcggcgcggcggtaaggtcgctactttgagttttggcttcttcaactctttcatcgccacccagagcagtatcaTCGTTCTTTCACCCCGAc                                                                                                                                                                                                                                        <  2:387952‑M2/18‑1 (MQ=255)
      aaattgcagaatattggccgttggggaagtgattactaacccaaacatgcgtccggttgcttcactcaggggccagtctaaaagcgcaatcaaaacgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggagggaagggcgcgaatttgagttttggagtcttcaactatttaagcaacaccaagagcgttaacaTTGTTCGTTAATTCCGAc                                                                                                                                                                                                                                        >  2:482618‑M2/234‑251 (MQ=255)
               aatattggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcgacgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccaccccgagcgttatcaTTGTTCGTTCATCCCGACTCCtttttt                                                                                                                                                                                                                               >  2:570831‑M2/225‑251 (MQ=255)
               aatagtggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCtttttt                                                                                                                                                                                                                               <  2:18631‑M2/27‑1 (MQ=255)
                  attggccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGta                                                                                                                                                                                                                            <  1:122889‑M2/30‑1 (MQ=255)
                       ccgttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTGCATCCCGACTCCTTTTTTGTATAGAt                                                                                                                                                                                                                       <  2:415579‑M2/35‑1 (MQ=255)
                       ccgttggggaagtgattcctaacccaaacatgcgtccgggtgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagctttatcaTTGGTCGTTCATCCCGACTCCTTTTTTGTATAGAt                                                                                                                                                                                                                       <  2:592621‑M2/35‑1 (MQ=255)
                         gttggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagggttatcaTTGTTCGTTCATCCCGACTCCTTTTTTTTATAGATa                                                                                                                                                                                                                      >  2:417183‑M2/215‑250 (MQ=255)
                            ggggaagtgattcctaacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAAc                                                                                                                                                                                                                   <  1:423921‑M2/39‑1 (MQ=255)
                                           aacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTgttcg                                                                                                                                                                                                                                                   <  2:357670‑M2/7‑1 (MQ=255)
                                           aacccaaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTgttcg                                                                                                                                                                                                                                                   >  1:357670‑M2/197‑203 (MQ=255)
                                                aaacatgcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaatcagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgaggtttggcttcttcaactctttcatcacctcccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTGTGTAGAGATAAACCATCAGCTGATAGTTTACCTg                                                                                                                                                                                              <  2:212770‑M2/60‑1 (MQ=255)
                                                  acatgcgtccggttgcttcactcagaggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcttttttcaaacagaatatgagcgaacccggttccgcgcggcggtaaggtcgcgactttgagttttggcttctgcaaccctttcatcaccacccagagcgttatcaTTGTTCGTTCATccc                                                                                                                                                                                                                                           <  1:392838‑M2/15‑1 (MQ=255)
                                                  acatgcgtacggctgcttcactcagggcacaggctaaaagcgaaatcaaaccgatgagtacagctgcgacaaaggggtgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggcaaggccgcgactttgagttttggcttcttcagatctttcatcaccacccagagcgttatcaTTGTTCGTTCATccc                                                                                                                                                                                                                                           >  2:392838‑M2/190‑204 (MQ=255)
                                                      gcgtccggttgcttcactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAat                                                                                                                                                                                        >  1:192674‑M2/186‑251 (MQ=255)
                                                      gcgtacggttgcttcactcaggggccaggctaaaagctcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAat                                                                                                                                                                                        <  1:343713‑M2/66‑1 (MQ=255)
                                                                    cactcaggggccaggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGtactt                                                                                                                                                                          <  2:161057‑M2/80‑1 (MQ=255)
                                                                         aggggccaggctaaaagcgcaatcaaaccgatgagtaaagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagagttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTacttactt                                                                                                                                                                      >  1:240230‑M2/167‑250 (MQ=255)
                                                                           gggccaggctcaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTCTAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGCGAAGTACTTACTTaa                                                                                                                                                                    >  2:181923‑M2/165‑250 (MQ=255)
                                                                                aggctaacagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTTGTTCATCCCGACTCCTTTTTTGTATAGATAAACCCTCCGCTGATAGTTTTCCTGGAGAATATAGAGAAGTACTTACTTAAAAtttt                                                                                                                                                              >  2:223924‑M2/160‑251 (MQ=255)
                                                                                aggctaaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACAtttt                                                                                                                                                              <  2:321469‑M2/92‑1 (MQ=255)
                                                                                    taaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTccc                                                                                                                                                           >  1:491315‑M2/156‑250 (MQ=255)
                                                                                    taaaagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTccc                                                                                                                                                           <  2:491315‑M2/95‑1 (MQ=255)
                                                                                       aagcgcaatcaaaccgatgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCAttt                                                                                                                                                       <  1:251661‑M2/99‑1 (MQ=255)
                                                                                         gcgcaatccacccgattagcactgctgtcacaaagggatgaccgcgtttttccaaccgaatatgagcgatcccggttcggcgcggcggtacggtcgcgactttgagttttggctttttcaactctttcatctcccccccgagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGCTACACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTgg                                                                                                                                                     >  2:379752‑M2/151‑251 (MQ=255)
                                                                                                       atgagtacagctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAAcccc                                                                                                                                       <  1:423461‑M2/115‑1 (MQ=255)
                                                                                                                gctgtgacaaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCActt                                                                                                                              >  1:387085‑M2/128‑250 (MQ=255)
                                                                                                                        aaagggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTaaga                                                                                                                        >  1:178632‑M2/120‑249 (MQ=255)
                                                                                                                          agggatgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggttaggtcgcgactttgagctttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCCTCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGt                                                                                                                     >  2:48250‑M2/118‑250 (MQ=255)
                                                                                                                              atgccagcgtttttcaaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGCCTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCCAACCCCTTTTCCCTCTTAAGAAGTAATGc                                                                                                                >  2:191345‑M2/114‑251 (MQ=255)
                                                                                                                                 ccagcgtttttcaaacagaatatgagcgatcccggttacgcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccacagcgttatcaTTGTTCGCTCATCCCGCCTCCTTTTTTGTATAGATAAACCCTCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTAATTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCct                                                                                                              >  1:399399‑M2/111‑250 (MQ=255)
                                                                                                                                            caaacagaatatgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTGCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATTCCTACTATGACTCa                                                                                                   <  2:72022‑M2/151‑1 (MQ=255)
                                                                                                                                                      atgagcgatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCgaag                                                                                         <  2:192674‑M2/161‑1 (MQ=255)
                                                                                                                                                            gatcccggttcggcgcggcggtaaggtcgcgactttgcgttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTc                                                                                  >  1:608413‑M2/84‑251 (MQ=255)
                                                                                                                                                            gatcccggttcggcgcggcggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAattt                                                                                  >  2:452666‑M2/84‑250 (MQ=255)
                                                                                                                                                                               ggtaaggtcgcgactttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATGCTGGTGACGTGCGCACTGcc                                                               <  1:190940‑M2/187‑1 (MQ=255)
                                                                                                                                                                                   aggtcgcgactttgagttttggcttcttcaactctttcatccccacccagagcgttatccTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGACGAAAATTCTGGTGACGTGCGCACTGCCGTAc                                                           >  2:30686‑M2/61‑251 (MQ=255)
                                                                                                                                                                                          gactttcagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACg                                                    >  1:38452‑M2/54‑251 (MQ=255)
                                                                                                                                                                                              ttgagttttggcttcttcaactctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCt                                                 <  1:570831‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                                  ctctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGc                            >  2:465865‑M2/30‑251 (MQ=255)
                                                                                                                                                                                                                  ctctttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGc                            >  1:168381‑M2/30‑251 (MQ=255)
                                                                                                                                                                                                                      ttcatcaccacccagagcgttatcaTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTgg                         >  2:515539‑M2/26‑250 (MQ=255)
                                                                                                                                                                                                                                             caTTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTAGTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAAGTCTGGGGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCGCCTCGGCCATATGCTGGTGCACATCCAGGCTGATGTCTgg  <  2:524842‑M2/248‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
TCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGTACGCTAACGGCTCAATCCACCTCGGCCATATGCTGGAGCACATCCAGGCTGATGTCTGG  >  NC_000913/2193918‑2194404

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: