Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,773,602 T→C V441A (GTT→GCT)  mtlA → mannitol‑specific PTS enzyme: IIA, IIB and IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,773,6020TC100.0% 265.1 / NA 74V441A (GTT→GCT) mtlAmannitol‑specific PTS enzyme: IIA, IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/36);  total (38/36)

GTAAGAAAATTCAGGATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCT  >  NC_000913/3773476‑3773736
                                                                                                                              |                                                                                                                                      
gTAAGAAAATTCAGGATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGcaca                                                                                                                            <  1:958913/139‑1 (MQ=255)
gTAAGAAAATTCAGGATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGcaca                                                                                                                            <  1:1201237/139‑1 (MQ=255)
    gAAAATTCAGGATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCat                                                                                                                        <  2:689483/139‑1 (MQ=255)
          tCAGGATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCg                                                                                                                  <  2:987732/139‑1 (MQ=255)
              gATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGa                                                                                                              >  1:1802042/1‑139 (MQ=255)
                   aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAg                                                                                                                                 <  2:193731/115‑1 (MQ=255)
                   aGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAg                                                                                                                                 >  1:193731/1‑115 (MQ=255)
                          tCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCg                                                                                                  <  2:1802042/139‑1 (MQ=255)
                             cAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACa                                                                                               <  1:1133361/139‑1 (MQ=255)
                                aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                            <  2:1596025/139‑1 (MQ=255)
                                aTTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                            <  1:1298729/139‑1 (MQ=255)
                                  ttCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGc                                                                                          >  2:138596/1‑139 (MQ=255)
                                     tGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGc                                                                                                                               <  2:282542/99‑1 (MQ=255)
                                     tGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGc                                                                                                                               >  1:282542/1‑99 (MQ=255)
                                             aCAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCGTCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAg                                                                               >  2:327943/1‑139 (MQ=255)
                                               aGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGcc                                                                             >  2:1844156/1‑139 (MQ=255)
                                                 cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCt                                                                                                    <  1:22984/114‑1 (MQ=255)
                                                 cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCt                                                                                                    >  2:22984/1‑114 (MQ=255)
                                                 cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                            <  1:1841632/122‑1 (MQ=255)
                                                 cgcgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                            >  2:1841632/1‑122 (MQ=255)
                                                   cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGca                                                                                                                              <  1:1712329/86‑1 (MQ=255)
                                                   cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGca                                                                                                                              >  2:1712329/1‑86 (MQ=255)
                                                   cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCg                                                                                                                  <  1:387437/98‑1 (MQ=255)
                                                   cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCg                                                                                                                  >  2:387437/1‑98 (MQ=255)
                                                   cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGAc                                                                         >  2:696248/1‑139 (MQ=255)
                                                   cgATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGAc                                                                         >  1:848968/1‑139 (MQ=255)
                                                            acCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGca                                                                                                                              >  2:1462630/1‑77 (MQ=255)
                                                            acCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGca                                                                                                                              <  1:1462630/77‑1 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCg                                                                                                                                    >  2:839619/1‑67 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCg                                                                                                                                    <  1:839619/67‑1 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCata                                                                                                                       >  1:1418618/1‑80 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCata                                                                                                                       <  2:1418618/80‑1 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGc                                                                                          <  1:1737926/109‑1 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGc                                                                                          >  2:1737926/1‑109 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                            >  1:1052724/1‑139 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                            >  2:1673117/1‑139 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                            >  2:1327816/1‑139 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                            >  2:552879/1‑139 (MQ=255)
                                                                gccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGtt                                                            >  2:1099003/1‑139 (MQ=255)
                                                                 ccgccAGATGTGGACCTCGTCGTCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtg                                                           <  1:327943/139‑1 (MQ=255)
                                                                 ccgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtg                                                           >  2:1021669/1‑139 (MQ=255)
                                                                  cgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACa                                                                                               <  1:554235/102‑1 (MQ=255)
                                                                  cgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACa                                                                                               >  2:554235/1‑102 (MQ=255)
                                                                  cgccAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTtgc                                                          >  2:1686857/1‑139 (MQ=255)
                                                                        aTGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCa                                                    >  1:609477/1‑139 (MQ=255)
                                                                                  cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCCCAGCATAGTTCGCTGACCCACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                          >  1:1649208/1‑139 (MQ=255)
                                                                                  cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGcc                                                                             >  1:1711978/1‑104 (MQ=255)
                                                                                  cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGcc                                                                             <  2:1711978/104‑1 (MQ=255)
                                                                                  cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                          >  1:1130905/1‑139 (MQ=255)
                                                                                  cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                          >  1:566436/1‑139 (MQ=255)
                                                                                  cGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACg                                          >  1:814739/1‑139 (MQ=255)
                                                                                   gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                            <  2:509213/88‑1 (MQ=255)
                                                                                   gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCg                                                                                            >  1:509213/1‑88 (MQ=255)
                                                                                   gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACgg                                         <  2:1052724/139‑1 (MQ=255)
                                                                                   gTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACgg                                         <  1:138596/139‑1 (MQ=255)
                                                                                    tcatcaCTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGATCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGc                                        >  2:1286523/1‑139 (MQ=255)
                                                                                         actcacCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACg                                   <  2:814739/139‑1 (MQ=255)
                                                                                             acCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGaa                                 <  1:780956/137‑1 (MQ=255)
                                                                                             acCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGaa                                 >  2:780956/1‑137 (MQ=255)
                                                                                             acCGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                               <  2:848968/139‑1 (MQ=255)
                                                                                               cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                             <  1:696248/139‑1 (MQ=255)
                                                                                               cGTGACCTGACCGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAgaga                             <  1:552879/139‑1 (MQ=255)
                                                                                                  gacctgaccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAg                          >  1:10617/1‑139 (MQ=255)
                                                                                                   acctgaccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCcac                                             >  2:1121270/1‑119 (MQ=255)
                                                                                                   acctgaccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCcac                                             <  1:1121270/119‑1 (MQ=255)
                                                                                                       gaccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                               >  1:868095/1‑129 (MQ=255)
                                                                                                       gaccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAaga                               <  2:868095/129‑1 (MQ=255)
                                                                                                         ccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCca                                              >  1:909452/1‑112 (MQ=255)
                                                                                                         ccGAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCca                                              <  2:909452/112‑1 (MQ=255)
                                                                                                           gAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                 <  1:1021669/139‑1 (MQ=255)
                                                                                                           gAACGCGCTATGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACa                 <  2:609477/139‑1 (MQ=255)
                                                                                                                     tGCGCCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGa       <  1:1413439/139‑1 (MQ=255)
                                                                                                                       cgcCAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACa     >  1:731933/1‑139 (MQ=255)
                                                                                                                          cAGGCTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCt  <  1:1673117/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                      
GTAAGAAAATTCAGGATGCAGGTCTGTCGCAGATTTCTGTTACTAACAGCGCGATCAACAACCTGCCGCCAGATGTGGACCTCGTCATCACTCACCGTGACCTGACCGAACGCGCTATGCGCCAGGTTCCGCAGGCACAGCATATTTCGCTGACCAACTTCCTCGACAGCGGCCTGTACACCAGCCTGACCGAACGTCTGGTTGCTGCCCAACGCCACACGGCAAACGAAGAGAAAGTAAAAGACAGCCTGAAAGACAGCT  >  NC_000913/3773476‑3773736

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