Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,465,279 T→C intergenic (‑272/‑22) mlaA ← / → yfdC ABC transporter maintaining OM lipid asymmetry, OM lipoprotein component/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,465,2790TC100.0% 257.8 / NA 68intergenic (‑272/‑22)mlaA/yfdCABC transporter maintaining OM lipid asymmetry, OM lipoprotein component/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/35);  total (33/35)

CATAGCCTGGTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCTCC  >  NC_000913/2465144‑2465405
                                                                                                                                       |                                                                                                                              
cATAGCCTGGTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCt                                                                                                                             <  2:1285128/139‑1 (MQ=255)
cATAGCCTGGTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCt                                                                                                                             >  2:1359716/1‑139 (MQ=255)
     ccTGGTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTc                                                                                                                        >  2:1219834/1‑139 (MQ=255)
     ccTGGTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTc                                                                                                                        >  2:820308/1‑139 (MQ=255)
     ccTGGTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGC‑‑TACGCGCTAATTCaa                                                                                                                      >  1:145781/1‑139 (MQ=255)
        ggTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTAGATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaag                                                                                                                     >  1:1044029/1‑139 (MQ=255)
        ggTAATCGCCGCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaaga                                                                                                                    >  2:1565493/1‑139 (MQ=255)
         gTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaaga                                                                                                                    >  2:1642704/1‑139 (MQ=255)
             tCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaagagaag                                                                                                                <  1:1207435/139‑1 (MQ=255)
             tCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaagagaag                                                                                                                <  2:1082157/139‑1 (MQ=255)
             tCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaagagaag                                                                                                                <  1:1359716/139‑1 (MQ=255)
             tCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaagagaag                                                                                                                <  1:820308/139‑1 (MQ=255)
               gccgccCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaagagaagag                                                                                                              >  2:1053439/1‑139 (MQ=255)
                  gccCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGAcc                                                                                                           <  2:746870/139‑1 (MQ=255)
                        tAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGAc                                                                                                     >  1:713019/1‑139 (MQ=255)
                          ggCTACTTGATTCCATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGATCATAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACaa                                                                                                   <  1:1180379/139‑1 (MQ=255)
                              aCTTGATTCTTTAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGAc                                                                                               >  1:685254/1‑139 (MQ=255)
                               cTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACa                                                                                              >  1:358179/1‑139 (MQ=255)
                                ttGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaa                                                                                             >  1:131331/1‑139 (MQ=255)
                                  gATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaaaa                                                                                           >  1:1598702/1‑139 (MQ=255)
                                  gATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaaaa                                                                                           <  2:1259319/139‑1 (MQ=255)
                                  gATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaaaa                                                                                           >  1:625490/1‑139 (MQ=255)
                                  gATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaaaa                                                                                           <  1:1053439/139‑1 (MQ=255)
                                       tataGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaa                                                                                             >  1:1646169/1‑132 (MQ=255)
                                       tataGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaa                                                                                             <  2:1646169/132‑1 (MQ=255)
                                       tataGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGAt                                                                                      <  1:1206511/139‑1 (MQ=255)
                                         taGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCa                                                                                    >  1:618887/1‑139 (MQ=255)
                                            aaaCAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAaca                                                                                 <  2:586217/139‑1 (MQ=255)
                                              aCAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAAcaca                                                                               >  1:1170988/1‑139 (MQ=255)
                                                  aaaaaGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGa                                                                                                <  2:1237685/118‑1 (MQ=255)
                                                  aaaaaGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGa                                                                                                >  1:1237685/1‑118 (MQ=255)
                                                     aaGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGcgacga                                                                        <  2:618887/139‑1 (MQ=255)
                                                       ggCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGaaa                                                                      <  1:1642704/139‑1 (MQ=255)
                                                            aCGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGaa                                                                 >  2:188905/1‑139 (MQ=255)
                                                                    ttttCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGagag                                                         <  2:131331/139‑1 (MQ=255)
                                                                      ttCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCg                                                       <  2:124972/139‑1 (MQ=255)
                                                                      ttCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCg                                                       <  2:1598702/139‑1 (MQ=255)
                                                                      ttCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCg                                                       <  2:1044029/139‑1 (MQ=255)
                                                                      ttCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCg                                                       <  1:155964/139‑1 (MQ=255)
                                                                       tCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACATAATTGATCAACACAGcga                                                                           <  1:992866/118‑1 (MQ=255)
                                                                       tCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGcga                                                                           >  2:992866/1‑118 (MQ=255)
                                                                       tCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCga                                                      >  1:1730309/1‑139 (MQ=255)
                                                                            aTTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGcgacga                                                                        >  2:867605/1‑116 (MQ=255)
                                                                            aTTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGcgacga                                                                        <  1:867605/116‑1 (MQ=255)
                                                                                        aCTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCa                                     <  1:1131317/139‑1 (MQ=255)
                                                                                         cTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaaaa                                                                                           <  2:1711102/84‑1 (MQ=255)
                                                                                         cTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACaaaa                                                                                           >  1:1711102/1‑84 (MQ=255)
                                                                                         cTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaa                                    >  1:868892/1‑139 (MQ=255)
                                                                                              gTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCa                                                                                                          <  2:1382355/64‑1 (MQ=255)
                                                                                              gTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCa                                                                                                          >  1:1382355/1‑64 (MQ=255)
                                                                                              gTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCaaaaaaa                               >  2:808701/1‑139 (MQ=255)
                                                                                                  gTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaataga                           <  2:1730309/139‑1 (MQ=255)
                                                                                                   ttGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaagagaaga                                                                                                               <  2:1263470/54‑1 (MQ=255)
                                                                                                   ttGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaagagaaga                                                                                                               >  1:1263470/1‑54 (MQ=255)
                                                                                                    tGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAaaatagaaa                         <  2:358179/139‑1 (MQ=255)
                                                                                                      cGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaa                                                                                                                      <  1:545359/44‑1 (MQ=255)
                                                                                                      cGAGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCaa                                                                                                                      >  2:545359/1‑44 (MQ=255)
                                                                                                        aGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGAt                                                                                      <  1:607843/74‑1 (MQ=255)
                                                                                                        aGCAAAAAAGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGAt                                                                                      >  2:607843/1‑74 (MQ=255)
                                                                                                           aaaaaaGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGa                  <  2:685254/139‑1 (MQ=255)
                                                                                                             aaaaGCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAg                >  1:1020236/1‑139 (MQ=255)
                                                                                                                 gCCGCTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCg            <  2:1170988/139‑1 (MQ=255)
                                                                                                                     cTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCga                                                      <  1:761717/93‑1 (MQ=255)
                                                                                                                     cTACGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCga                                                      >  2:761717/1‑93 (MQ=255)
                                                                                                                       aCGCTTTAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGActccc      >  2:1416273/1‑139 (MQ=255)
                                                                                                                           tttAGCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGTGCGGCAAAAAAATAGAAATAGATGAAGACCGActcccctcc  <  1:1565493/139‑1 (MQ=255)
                                                                                                                               gCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCg                                                       >  1:1481047/1‑82 (MQ=255)
                                                                                                                               gCTATACGCGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCg                                                       <  2:1481047/82‑1 (MQ=255)
                                                                                                                                       |                                                                                                                              
CATAGCCTGGTAATCGCCGCCCTGTAGGCTACTTGATTCTATAGAAACAGAAAAAGGCAAACGCAGCCTTTTCTACATTTTCAGAGTAACTCCCGTCAGTTGCGAGCAAAAAAGCCGCTACGCTTTAGCTATACGTGCTAATTCAAGAGAAGAGACCATGGACAACGACAAAATTGATCAACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCTCC  >  NC_000913/2465144‑2465405

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: