Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,870,441 G→A *355* (TAG→TAA) 
G430G (GGC→GGT) 		 cbrA →
dgoT ←		 colicin M resistance protein, FAD‑binding protein, putative oxidoreductase
D‑galactonate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,870,4410GA100.0% 245.3 / NA 75*355* (TAG→TAA) 
G430G (GGC→GGT) 		cbrA
dgoT		colicin M resistance protein, FAD‑binding protein, putative oxidoreductase
D‑galactonate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/40);  total (35/40)

CTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGG  >  NC_000913/3870305‑3870571
                                                                                                                                        |                                                                                                                                  
cTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAAcc                                                                                                                                  >  1:727327/1‑139 (MQ=255)
     tGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAcgc                                                                                                                             >  1:319954/1‑139 (MQ=255)
     tGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAcgc                                                                                                                             >  2:540882/1‑139 (MQ=255)
          cgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTc                                                                                                                        >  2:409867/1‑139 (MQ=255)
          cgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTc                                                                                                                        >  2:1579997/1‑139 (MQ=255)
          cgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGACTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCaa                                                                                                                                >  2:427301/1‑131 (MQ=255)
          cgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGACTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCaa                                                                                                                                <  1:427301/131‑1 (MQ=255)
           gcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTca                                                                                                                       <  2:566067/139‑1 (MQ=255)
             cACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTcaca                                                                                                                     <  2:853893/139‑1 (MQ=255)
                        tGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  2:727327/139‑1 (MQ=255)
                        tGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  1:540882/139‑1 (MQ=255)
                        tGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  2:1558103/139‑1 (MQ=255)
                          cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  2:165769/139‑1 (MQ=255)
                                  aCTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  1:902721/129‑1 (MQ=255)
                                  aCTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          >  2:902721/1‑129 (MQ=255)
                                      ttCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAagaga                                                                                            <  1:638814/139‑1 (MQ=255)
                                      ttCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAagaga                                                                                            >  1:168780/1‑139 (MQ=255)
                                          gCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcgc                                                                                        <  1:754306/139‑1 (MQ=255)
                                                tAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAAt                                                                                  >  2:1362819/1‑139 (MQ=255)
                                                         gCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        >  1:1569540/1‑108 (MQ=255)
                                                         gCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  2:1569540/108‑1 (MQ=255)
                                                          ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAc                                                                                                                               >  1:837896/1‑84 (MQ=255)
                                                          ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAAc                                                                                                                               <  2:837896/84‑1 (MQ=255)
                                                          ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcg                                                                        >  2:1024224/1‑139 (MQ=255)
                                                          ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcg                                                                        >  2:580689/1‑139 (MQ=255)
                                                          ccGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcg                                                                        >  2:1485928/1‑139 (MQ=255)
                                                             aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  1:312142/104‑1 (MQ=255)
                                                             aTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        >  2:312142/1‑104 (MQ=255)
                                                             aTGCATGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                     <  1:1469182/139‑1 (MQ=255)
                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  1:584369/139‑1 (MQ=255)
                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  1:1362819/139‑1 (MQ=255)
                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  2:963977/139‑1 (MQ=255)
                                                              tGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  1:1024224/139‑1 (MQ=255)
                                                                   gACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                               <  1:322774/139‑1 (MQ=255)
                                                                   gACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                               <  2:1605692/139‑1 (MQ=255)
                                                                         aCCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  1:623197/92‑1 (MQ=255)
                                                                         aCCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        >  2:623197/1‑92 (MQ=255)
                                                                             gCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTg                                                     >  1:705/1‑139 (MQ=255)
                                                                              cTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGc                                                    <  2:319954/139‑1 (MQ=255)
                                                                                      aaGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                            <  2:519016/139‑1 (MQ=255)
                                                                                       aGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                               <  1:1467837/119‑1 (MQ=255)
                                                                                       aGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                               >  2:1467837/1‑119 (MQ=255)
                                                                                            aTTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGt                                      >  1:772187/1‑139 (MQ=255)
                                                                                                   gcACTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGCGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTg                               >  1:1036376/1‑139 (MQ=255)
                                                                                                     aGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          >  2:776153/1‑62 (MQ=255)
                                                                                                     aGTGGTGTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcag                                                                                                          <  1:776153/62‑1 (MQ=255)
                                                                                                       tggtgTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGcc                           <  2:772187/139‑1 (MQ=255)
                                                                                                       tggtgTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGcc                           <  1:1579997/139‑1 (MQ=255)
                                                                                                        ggtgTGGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCa                          <  1:1485928/139‑1 (MQ=255)
                                                                                                          tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        <  1:637052/59‑1 (MQ=255)
                                                                                                          tgtgGCACATATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACcagca                                                                                                        >  2:637052/1‑59 (MQ=255)
                                                                                                                 caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    >  2:1591612/1‑88 (MQ=255)
                                                                                                                 caTATTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  1:1591612/88‑1 (MQ=255)
                                                                                                                   tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGTGCGCCAATCAACGCGAcggcgg                                                                    <  2:929603/86‑1 (MQ=255)
                                                                                                                   tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGTGCGCCAATCAACGCGAcggcgg                                                                    >  1:929603/1‑86 (MQ=255)
                                                                                                                   tatTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCccacc               >  2:1541994/1‑139 (MQ=255)
                                                                                                                    atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                            >  1:16303/1‑109 (MQ=255)
                                                                                                                    atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                            <  2:16303/109‑1 (MQ=255)
                                                                                                                    atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGt                       >  2:1069703/1‑130 (MQ=255)
                                                                                                                    atTCCACAGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGt                       <  1:1069703/130‑1 (MQ=255)
                                                                                                                        cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcgcg                                                                                       <  1:431303/62‑1 (MQ=255)
                                                                                                                        cacaGTTGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAgcgcg                                                                                       >  2:431303/1‑62 (MQ=255)
                                                                                                                             ttGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTaa                                    >  1:1108786/1‑108 (MQ=255)
                                                                                                                             ttGAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTaa                                    <  2:1108786/108‑1 (MQ=255)
                                                                                                                               gAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                               >  1:136676/1‑79 (MQ=255)
                                                                                                                               gAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                               <  2:136676/79‑1 (MQ=255)
                                                                                                                               gAAGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATAAAAGAGAGCGAGCCAATCAACGCGAcggcg                                                                     <  2:1075949/73‑1 (MQ=255)
                                                                                                                               gAAGGATTAACCAACCCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcg                                                                     >  1:1075949/1‑73 (MQ=255)
                                                                                                                                aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  1:240641/73‑1 (MQ=255)
                                                                                                                                aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    >  2:240641/1‑73 (MQ=255)
                                                                                                                                aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCccaccac             >  1:803510/1‑128 (MQ=255)
                                                                                                                                aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCccaccac             <  2:803510/128‑1 (MQ=255)
                                                                                                                                aaGGATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACgg  >  2:376171/1‑139 (MQ=255)
                                                                                                                                  ggATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    <  1:7521/71‑1 (MQ=255)
                                                                                                                                  ggATTAACCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGAcggcgg                                                                    >  2:7521/1‑71 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGG  >  NC_000913/3870305‑3870571

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