Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,048,805 Δ48,085 bp [gfcB]ycdU 45 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1048805 1096889 48085 39 [0] [0] 39 [gfcB]–ycdU [gfcB],gfcA,insA4,insB4,cspH,cspG,ymcF,ymcE,gnsA,yccM,torS,torT,torR,torC,torA,torD,cbpM,cbpA,yccE,agp,yccJ,wrbA,ymdF,rutG,rutF,rutE,rutD,rutC,rutB,rutA,rutR,putA,ymdG,putP,efeO,efeB,phoH,pgaD,pgaC,pgaB,pgaA,dgcT,insF4,insE4,ycdU

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 10488040 (0.000)39 (0.770) 37/488 0.2 100% coding (529/645 nt) gfcB lipoprotein GfcB
?NC_000913 1096890 = 0 (0.000)intergenic (+61/+675) ycdU/serX uncharacterized protein YcdU/tRNA‑Ser

AACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCCAGCAGATTATTATTAAGTTTGACGG  >  NC_000913/1048557‑1049052
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
aacaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCCGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGatt                                                                                                                                                                                                                                                       >  2:210978‑M1/1‑248 (MQ=255)
 acaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattg                                                                                                                                                                                                                                                      >  1:140998‑M1/1‑247 (MQ=255)
   aTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgtt                                                                                                                                                                                                                                                    <  1:209746‑M1/251‑7 (MQ=255)
     tttCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATAAGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGCATGTGTTTTGCCACTGattgttg                                                                                                                                                                                                                                                   >  2:50571‑M1/1‑243 (MQ=255)
             gggCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaa                                                                                                                                                                                                                                          <  1:339932‑M1/251‑17 (MQ=255)
               gCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGATTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTGTTGCCACTGattgttgactgagtaacc                                                                                                                                                                                                                                        <  1:329165‑M1/251‑19 (MQ=255)
                   aCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtacccagac                                                                                                                                                                                                                                    >  2:425439‑M1/1‑229 (MQ=255)
                   aCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagac                                                                                                                                                                                                                                    >  1:281767‑M1/1‑229 (MQ=255)
                       aaaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagtt                                                                                                                                                                                                                                <  1:138730‑M1/251‑27 (MQ=255)
                        aaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttg                                                                                                                                                                                                                               >  1:409125‑M1/1‑224 (MQ=255)
                            gtCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatg                                                                                                                                                                                                                            <  1:300672‑M1/250‑31 (MQ=255)
                              cTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgc                                                                                                                                                                                                                         >  1:409806‑M1/1‑218 (MQ=255)
                                gttgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcac                                                                                                                                                                                                                       <  2:409806‑M1/251‑36 (MQ=255)
                                gttgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcac                                                                                                                                                                                                                       >  1:448133‑M1/1‑216 (MQ=255)
                                   gctcaccaggcaGATAAAGTGTCTCCGTTTCTTGCGCAAAGGCAGGGGGTGTCATTACGTAAAGTTCGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGTCGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattggtgactgagtaaccagacagttgatgttcacgat                                                                                                                                                                                                                    <  2:284773‑M1/246‑39 (MQ=255)
                                               gATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgattgcccctcgccc                                                                                                                                                                                                        <  2:391205‑M1/251‑51 (MQ=255)
                                                        tCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGGCCGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGGCATTTCTACGGGAATAACCCCTGCGCCTAACTTTTGCCGCCTTTGACGAACCTGGCCAGTAGCCGAATAGACCCAGAATGTGTTTTGCCACTGattggtgactgagtaaccaaacagttgatgtgcacgaattcccctagccctatcaaacg                                                                                                                                                                                               >  2:316592‑M1/1‑192 (MQ=255)
                                                           ccGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGGACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACAGTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgg                                                                                                                                                                                            <  2:391182‑M1/251‑63 (MQ=255)
                                                               tCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgg                                                                                                                                                                                        >  1:85929‑M1/1‑185 (MQ=255)
                                                                ccTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgg                                                                                                                                                                                        <  1:83902‑M1/250‑67 (MQ=255)
                                                                  ttGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtggtcggggg                                                                                                                                                                                     >  2:7834‑M1/1‑182 (MQ=255)
                                                                           ggCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGAGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtg                                                                                                                                                                            >  1:260662‑M1/1‑173 (MQ=255)
                                                                            gCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgt                                                                                                                                                                           <  1:397578‑M1/251‑80 (MQ=255)
                                                                                ggggTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgc                                                                                                                                                                       <  2:352743‑M1/251‑84 (MQ=255)
                                                                                       cATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagcca                                                                                                                                                                                                                                   >  2:3174‑M1/1‑161 (MQ=255)
                                                                                       cATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagaca                                                                                                                                                                                                                                   <  1:3174‑M1/184‑24 (MQ=255)
                                                                                                 aGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaa                                                                                                                                                        >  2:436164‑M1/1‑151 (MQ=255)
                                                                                                                aaaTACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtac                                                                                                                                       <  2:126945‑M1/251‑116 (MQ=255)
                                                                                                                          gTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttct                                                                                                                             <  1:360252‑M1/251‑126 (MQ=255)
                                                                                                                          gTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttct                                                                                                                             <  1:84370‑M1/251‑126 (MQ=255)
                                                                                                                             aTTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttcttt                                                                                                                           >  1:132630‑M1/1‑123 (MQ=255)
                                                                                                                                 aTTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttcgttaaata                                                                                                                      <  1:214966‑M1/251‑133 (MQ=255)
                                                                                                                                       ggTGCGGGTTTAAGAACCGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcgg                                                                                                                >  2:366449‑M1/1‑113 (MQ=255)
                                                                                                                                                           cATTTCTACGGGAATAACCCCTGCGCCTAGCATTTGCCGACTTTGACGAACCTGGTCAGTAGCCGAATCGACCCATAATGTGTTTTGCCACTGatttttgacttagtatcccgacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcaccctgtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctgtaaatattgcggtgagggggggattgactcg                                                                                             <  2:177464‑M1/250‑158 (MQ=255)
                                                                                                                                                           cATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGTCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgggggggggagtcacgcgg                                                                                            >  2:425221‑M1/1‑93 (MQ=255)
                                                                                                                                                                       attaaCCCCTGCGCCTAACATTGTCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggattgactcgcttcgctcgccct                                                                                <  2:140998‑M1/249‑171 (MQ=255)
                                                                                                                                                                            ccccTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggattgactcgcttcgctcgccctgcggg                                                                           <  1:425439‑M1/251‑176 (MQ=255)
                                                                                                                                                                                                                           gAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgagggggggggtgaaacgctttcctcgcccttcgggaagcccgccccccgcggttacggtctgtccaacttggtgtcccccgt                            >  1:201394‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                                     cTGattgttgactgagtaaccagacagttgatgtgcacgatttcccctcgccctaacagacgtgggcgggggcacccggtgtgtgcaatatacgaaaaaaaagcccgtactttcgtacgagctcttctttaaatatggcggtgaggggggggttgactcgcttcgccccccccgcgggcagcccgctcaccgcgttcacggtctgtccaactggctgtcgcccgttggcgagcccccgtggggggttctcc  >  1:353907‑M1/1‑3 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
AACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGACATTACACAGCACATTGCGACCCGCAATATTCAGCGTCTCTAATCCAGCAAAGGAAAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCCAGCAGATTATTATTAAGTTTGACGG  >  NC_000913/1048557‑1049052

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