Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,638,037 C→A 18.6% A206S (GCG→TCG)  bamB ← BamABCDE complex OM biogenesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,638,0370CA18.6% 212.7 / 38.0 91A206S (GCG→TCG) bamBBamABCDE complex OM biogenesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base A (8/9);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

AGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGC  >  NC_000913/2637939‑2638136
                                                                                                  |                                                                                                   
aGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCa                                                                                                    <  1:1869094/100‑1 (MQ=255)
 gAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCAc                                                                                                   <  1:481055/100‑1 (MQ=255)
  aaCCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCAcc                                                                                                  <  1:2403110/100‑1 (MQ=255)
   aCCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCa                                                                                                 >  1:427980/1‑100 (MQ=255)
   aCCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCa                                                                                                 >  1:1395724/1‑100 (MQ=255)
    ccGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCaa                                                                                                >  1:1751031/1‑100 (MQ=255)
     cGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCaaa                                                                                               >  1:981203/1‑100 (MQ=255)
      ggTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCTCCCCCCACGACGGCCGCACCaaaa                                                                                              >  2:267963/1‑100 (MQ=255)
       gTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAg                                                                                             >  1:2299876/1‑100 (MQ=255)
         cGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGc                                                                                            >  1:2250293/1‑99 (MQ=255)
          gCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCg                                                                                          >  1:522500/1‑100 (MQ=255)
            cTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGtt                                                                                        >  2:83203/1‑100 (MQ=255)
                 aaaTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTTCGCTGACGCGACCATTATCGCCCCCCACGACGGAAACACCAAAAGCAGTTGTccg                                                                                   >  1:669929/1‑98 (MQ=255)
                  aaTCCGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGgc                                                                                  <  2:1485724/100‑1 (MQ=255)
                   aTACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGgcg                                                                                 <  1:707065/100‑1 (MQ=255)
                    tACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGgcgc                                                                                >  2:434723/1‑100 (MQ=255)
                    tACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGgcgc                                                                                <  1:2255276/100‑1 (MQ=255)
                     aCGCTGCTGCCAAATCATCTGGCCCTGTTCTATTTGTATTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCa                                                                               <  1:296631/100‑1 (MQ=255)
                     aCGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCa                                                                               <  1:1531948/100‑1 (MQ=255)
                     aCGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCa                                                                               >  1:988254/1‑100 (MQ=255)
                      cGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCCTTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAg                                                                              >  1:1392570/1‑100 (MQ=255)
                      cGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAg                                                                              >  1:688644/1‑100 (MQ=255)
                        ctgctgCCAAATCATCTGGCTCTGTTGCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGAc                                                                            <  1:1862288/100‑1 (MQ=255)
                          gctgcCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGACCGAACCAAAAGCCGTTGACGGCGCAGACt                                                                           >  2:773729/1‑99 (MQ=255)
                           ctgcCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCg                                                                         <  2:504270/100‑1 (MQ=255)
                            tgcCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGc                                                                        <  1:1424069/100‑1 (MQ=255)
                             gcCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGcc                                                                       >  1:1311795/1‑100 (MQ=255)
                                aaaTCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACg                                                                    <  2:1220891/100‑1 (MQ=255)
                                 aaTCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGc                                                                   <  2:2179561/100‑1 (MQ=255)
                                     atcTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAag                                                               <  2:111594/100‑1 (MQ=255)
                                       cTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAagag                                                             >  1:1325434/1‑100 (MQ=255)
                                         ggCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAagagag                                                           >  1:1073409/1‑100 (MQ=255)
                                          gCCCTGTTCCATCAGCACTGCGCGGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGc                                                          <  1:335213/100‑1 (MQ=255)
                                           cccTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCg                                                         <  2:870780/100‑1 (MQ=255)
                                            ccTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGa                                                        <  2:1732217/100‑1 (MQ=255)
                                             cTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGaa                                                       <  1:2120287/100‑1 (MQ=255)
                                              tGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAg                                                      <  1:561712/100‑1 (MQ=255)
                                               gctCCAACAGCACTGCGCTGACGCGACCATTATCGCCCCCTACGACGGCCGCACCAAAATCCGTTTTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAgg                                                     <  2:457838/98‑1 (MQ=255)
                                               gTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAgg                                                     >  1:323446/1‑100 (MQ=255)
                                               gTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAgg                                                     <  1:1958398/100‑1 (MQ=255)
                                                ttCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGTCTCGCCACGCAAAGAGAGCGAAGGc                                                    <  2:1485332/100‑1 (MQ=255)
                                                ttCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTTGCCACGCAAAGAGAGCGAAGGc                                                    <  2:1389852/100‑1 (MQ=255)
                                                ttCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGc                                                    <  2:2299849/100‑1 (MQ=255)
                                                ttCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGc                                                    <  1:2073110/100‑1 (MQ=255)
                                                 tCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCa                                                   <  1:483744/100‑1 (MQ=255)
                                                 tCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCa                                                   >  2:881675/1‑100 (MQ=255)
                                                  ccATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGAAGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCat                                                  >  2:2223080/1‑100 (MQ=255)
                                                  acATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGACGGCat                                                  >  2:1964829/2‑100 (MQ=255)
                                                   cATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAACGAGAGCGAAGGCat                                                  >  1:1972575/1‑99 (MQ=255)
                                                   cATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCata                                                 <  1:1611675/100‑1 (MQ=255)
                                                   cATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCata                                                 >  2:896393/1‑100 (MQ=255)
                                                    aTCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCatat                                                >  2:1883184/1‑100 (MQ=255)
                                                    aTCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCatat                                                >  1:1315178/1‑100 (MQ=255)
                                                    aTCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCatat                                                >  2:183646/1‑100 (MQ=255)
                                                    aTCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCatat                                                >  2:1076693/1‑100 (MQ=255)
                                                     tCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATc                                               >  2:1884131/1‑100 (MQ=255)
                                                       aGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGa                                             >  1:2200187/1‑100 (MQ=255)
                                                        gCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAg                                            >  2:1080097/1‑100 (MQ=255)
                                                           cTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAACAGCCGTTGTCGGCGCAGATTCGCCACGCAAAGAGAGCGAAGGCATATCGAgg                                           <  2:2293443/98‑1 (MQ=255)
                                                           cTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGtt                                         >  2:19187/1‑100 (MQ=255)
                                                            tGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTa                                        <  1:2185782/100‑1 (MQ=255)
                                                            tGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTa                                        >  1:1365265/1‑100 (MQ=255)
                                                             gcgcTGGCGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTaa                                       <  1:880970/100‑1 (MQ=255)
                                                              cgcTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAAc                                      <  1:2040408/100‑1 (MQ=255)
                                                               gcTGACGTGACCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAAc                                      <  1:83205/99‑1 (MQ=255)
                                                                  gACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTc                                  >  1:551943/1‑100 (MQ=255)
                                                                    cgcgACCATTCTCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCa                                >  2:2336832/1‑100 (MQ=255)
                                                                      cgACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCAtt                              >  2:662327/1‑100 (MQ=255)
                                                                        aCCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTg                            <  2:170706/100‑1 (MQ=255)
                                                                        aCCATTATCGCCCCCCACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTg                            <  2:28789/100‑1 (MQ=255)
                                                                         ccATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTg                            <  1:552483/99‑1 (MQ=255)
                                                                          cATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGAc                          <  2:613549/100‑1 (MQ=255)
                                                                           aTTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACa                         <  2:301560/100‑1 (MQ=255)
                                                                           aTTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACa                         <  2:695430/100‑1 (MQ=255)
                                                                               tCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAgcgc                     <  2:427977/100‑1 (MQ=255)
                                                                                 gCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCg                   >  2:2016252/1‑100 (MQ=255)
                                                                                 gCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCg                   >  1:484697/1‑100 (MQ=255)
                                                                                   cccccACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAACGAGCGCGAAGGCATATCGAGGATAACTGTCCATTTGACAGCGCCGTc                 >  1:497185/1‑100 (MQ=255)
                                                                                   cccccACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTc                 >  2:337832/1‑100 (MQ=255)
                                                                                   cccccACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTc                 <  1:1812013/100‑1 (MQ=255)
                                                                                    ccccACGACGTCCGCACCAAAAGCCGTGGTCGGCGCATACTTGCCACGCAAAGAGAGCGAAAGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCa                <  2:1481153/100‑1 (MQ=255)
                                                                                     cccACGACGGCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAg               <  1:19188/100‑1 (MQ=255)
                                                                                      ccACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGc              >  2:43482/1‑100 (MQ=255)
                                                                                           acgGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGt         <  1:183649/100‑1 (MQ=255)
                                                                                           acgGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGt         <  1:849936/100‑1 (MQ=255)
                                                                                             ggCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTc       <  1:285050/100‑1 (MQ=255)
                                                                                              gCCGAACCAAAAGCCGTTGTCGGCGCAGACTCGCTACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCa      <  1:705476/100‑1 (MQ=255)
                                                                                               ccGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAg     >  2:2048218/1‑100 (MQ=255)
                                                                                               ccGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAg     >  1:2062727/1‑100 (MQ=255)
                                                                                               ccGAACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAg     >  1:1147009/1‑100 (MQ=255)
                                                                                                cGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCa      <  1:2193324/98‑1 (MQ=255)
                                                                                                  cACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAgcgc  >  1:1514392/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGAACCGGTCGCCTGGGAAATACGCTGCTGCCAAATCATCTGGCCCTGTTCCATCAGCACTGCGCTGACGCGACCATTATCGCCCCCCACGACGGCCGCACCAAAAGCCGTTGTCGGCGCAGACTCGCCACGCAAAGAGAGCGAAGGCATATCGAGGTTAACTGTCCATTTGACAGCGCCGTCAGCTTCGTTCAGCGC  >  NC_000913/2637939‑2638136

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: