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 44 [0] [1] 45 [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)44 (0.590) 42/486 0.5 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

CACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCACCAGCAAAACGGCAACCAACGGCCACACAGACAAATTAAAAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATAT  >  NC_000913/2084796‑2085281
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
cACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGACGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgtt                                                                                                                                                                                                                                             >  1:598218‑M1/1‑246 (MQ=255)
 aCCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGGTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttc                                                                                                                                                                                                                                            <  1:592656‑M1/251‑7 (MQ=255)
           aGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccc                                                                                                                                                                                                                                  >  2:13011‑M1/1‑235 (MQ=255)
            gAAATTGCAGAATATTGGCCGCTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCAAGGCTAAAAGCGCAATCAAACCGCTGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTTGCGACTTTGAGTTTTGGCTTTTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgatcatcccg                                                                                                                                                                                                                                  >  1:477542‑M1/1‑234 (MQ=255)
                       atatTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttg                                                                                                                                                                                                                      >  1:99954‑M1/1‑223 (MQ=255)
                                                    cccAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttac                                                                                                                                                                                         >  1:602956‑M1/1‑194 (MQ=255)
                                                        aaCATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccgtttttgtatagataaaccatcagcggatagtttacctga                                                                                                                                                                                     <  1:416083‑M1/251‑62 (MQ=255)
                                                            agcagccggttgCATCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgtgcgtgcatcccgacgcctttttggtatagataaaccatcagctgatagtttacctgaagaa                                                                                                                                                                                 <  2:76801‑M1/246‑66 (MQ=255)
                                                                    gTTGCTTCACTCAGGGGCCAGGCTCAATGCGCAATCAAACCGATGAGTACAGCTGTGTCACAGGGATGCCAGCGTTTTTCAAACAGACTATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGTGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCCttgttcgttcatcccgactccttttttgtatagatcaaccatcagctgctagtttacctgaagaatatag                                                                                                                                                                            >  2:666665‑M1/1‑178 (MQ=255)
                                                                       gCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTGCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttgtttgtatagataaaccatcagctgatagttttcctgaagaatatagagaag                                                                                                                                                                       <  2:451130‑M1/250‑76 (MQ=255)
                                                                                            aaaaGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACGTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCTGAGCGTTATCAttgttcgttcatctcgactcctgttttgtatagagaaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttccca                                                                                                                                                  <  1:276254‑M1/250‑97 (MQ=255)
                                                                                             aaaGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatt                                                                                                                                                <  2:571002‑M1/251‑99 (MQ=255)
                                                                                                gcgcAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGCGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttgg                                                                                                                                             >  2:237590‑M1/1‑150 (MQ=255)
                                                                                                 cgcAATCAAACCGATGAGTACAGCTGGGACAAAGGGATGACAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGGTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccagttggt                                                                                                                                            <  2:151058‑M1/251‑103 (MQ=255)
                                                                                                       cAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATCTGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagctaaaccatccgctgatagtttacctgaagaatatagagaagtaattacttaacattttcccatttggtactatc                                                                                                                                      >  2:310355‑M1/1‑143 (MQ=255)
                                                                                                            cGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTGTCAAACAGAATAGGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCATAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagcggatagtttacctgaagaatagagagaaggacttacttaacattttcccatttggtactatctaacc                                                                                                                                 <  1:444244‑M1/251‑114 (MQ=255)
                                                                                                                    aCAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttca                                                                                                                         <  2:151833‑M1/251‑122 (MQ=255)
                                                                                                                            gACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaag                                                                                                                 <  2:652018‑M1/251‑130 (MQ=255)
                                                                                                                            gACAAAGGCATGCCAGCGTCTTTCAACCAGACTATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTAGCGACTTTGAGCTTTGGCTTCTTCAACTCTTTGATCCCCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtctagataacccatcagctggtagtttacctgaagaatatagagaagtccttacttaacattttcccatttggtactatctaaccccttttcactatta                                                                                                                   >  2:430151‑M1/1‑122 (MQ=255)
                                                                                                                               aaaGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaag                                                                                                              <  2:368415‑M1/251‑133 (MQ=255)
                                                                                                                                    gATGCCAGTGTTTTTCCAACTGAAAATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGCCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAAAGCGGTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctga                                                                                                                                                                                     >  1:516500‑M1/1‑114 (MQ=255)
                                                                                                                                    gATGCCAGCGTTTTTCAAACAGAATATGAGCTATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactcctttgttgtagagataaaccatcagctgatagtttacctga                                                                                                                                                                                     <  2:516500‑M1/175‑62 (MQ=255)
                                                                                                                                     aTGCCAGCGTTTGTCAAACAGAATATGAGCGATCCCGGTTCGTCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCTCCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagttcttacttaacattttcccatttggtactatctaaccccttttctctattaagaagtaatgc                                                                                                        <  1:318246‑M1/251‑139 (MQ=255)
                                                                                                                                         cAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctac                                                                                                    <  1:571426‑M1/251‑143 (MQ=255)
                                                                                                                                                     aaCAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGAATTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagt                                                                                        >  2:266516‑M1/1‑97 (MQ=255)
                                                                                                                                                       cAGAATCTGAGCGATCCCAGTTCGGCGCGGCGGTAAGGTAGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAACGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccctcagctgatagtttacctgtagactatcgagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattcagaagtaatgcctactctgactccagtcg                                                                                      >  2:449537‑M1/1‑95 (MQ=255)
                                                                                                                                                       cAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcg                                                                                      >  1:100483‑M1/1‑95 (MQ=255)
                                                                                                                                                            tatGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaataagtaatgcctactatgactcaagtcgcga                                                                                   >  2:488271‑M1/1‑90 (MQ=255)
                                                                                                                                                                   gATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTCTTGGCTTCTTAAACTCTTTAATCCCCACCCAGAGCGTTCTCAttgttcgttcctcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattc                                                                          >  1:36157‑M1/1‑83 (MQ=255)
                                                                                                                                                                                     cggTAAGGTCGCGACTGTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcggtcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgc                                                        <  2:701‑M1/251‑187 (MQ=255)
                                                                                                                                                                                             tCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgct                                                >  1:351327‑M1/1‑57 (MQ=255)
                                                                                                                                                                                                       gAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactagtaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctca                                       <  1:292766‑M1/250‑204 (MQ=255)
                                                                                                                                                                                                         gTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatc                                    >  2:23598‑M1/1‑45 (MQ=255)
                                                                                                                                                                                                         gTTTTGGCTTCTTCAACTCTTCCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactcctttattgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactatcaagaaggaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatc                                    <  1:332835‑M1/251‑207 (MQ=255)
                                                                                                                                                                                                                       aaCTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagctaaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatat                      >  1:422713‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                       aaCTCTTTCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatat                      <  1:13011‑M1/251‑221 (MQ=255)
                                                                                                                                                                                                                             ttCATCACCACCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatatgctgg                 <  1:23598‑M1/250‑226 (MQ=255)
                                                                                                                                                                                                                                  caccacCCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatatgctggagcaca           <  1:233178‑M1/251‑232 (MQ=255)
                                                                                                                                                                                                                                   cccacccAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatatgctggagcacat          <  2:100483‑M1/250‑233 (MQ=255)
                                                                                                                                                                                                                                    ccaccCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatatgctggagcacatc         <  1:158389‑M1/251‑234 (MQ=255)
                                                                                                                                                                                                                                      accCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagta                                                                                                            >  2:481117‑M1/1‑16 (MQ=255)
                                                                                                                                                                                                                                      accCAGAGCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagta                                                                                                            <  1:481117‑M1/150‑135 (MQ=255)
                                                                                                                                                                                                                                          agagCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatatgctggagcacatccaggct   <  1:237590‑M1/251‑240 (MQ=255)
                                                                                                                                                                                                                                           gagCGTTATCAttgttcgttcatcccgactccttttttgtatagataaaccatcagctgatagtttacctgaagaatatagagaagtacttacttaacattttcccatttggtactatctaaccccttttcactattaagaagtaatgcctactatgactcaagtcgcgaagaaaattctggtgacgtgcgcactgccgtacgctaacggctcaatccacctcggccatatgctggagcacatccaggctg  >  1:276445‑M1/1‑11 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
CACCCGCGACCAGAAATTGCAGAATATTGGCCGTTGGGGAAGTGATTCCTAACCCAAACATGCGTCCGGTTGCTTCACTCAGGGGCCAGGCTAAAAGCGCAATCAAACCGATGAGTACAGCTGTGACAAAGGGATGCCAGCGTTTTTCAAACAGAATATGAGCGATCCCGGTTCGGCGCGGCGGTAAGGTCGCGACTTTGAGTTTTGGCTTCTTCAACTCTTTCATCACCACCCAGAGCGTTATCACCAGCAAAACGGCAACCAACGGCCACACAGACAAATTAAAAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAACCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATAT  >  NC_000913/2084796‑2085281

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: