Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,139,720 C→A 6.5% intergenic (+235/+39) yegH → / ← asmA inner membrane protein/suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,139,7200CA 6.5% 176.8 / 3.0 69intergenic (+235/+39)yegH/asmAinner membrane protein/suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (24/40);  new base A (3/2);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

AGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTG  >  NC_000913/2139574‑2139864
                                                                                                                                                  |                                                                                                                                                
aGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcgc                                                                                                                                                >  2:335412/1‑149 (MQ=255)
  ggtggtGGATATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCTTTCCATATTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGt                                                                                                                                              <  1:477531/149‑1 (MQ=255)
     ggtggAGATCGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCtt                                                                                                                                           <  1:680861/149‑1 (MQ=255)
      gtggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                          >  2:370475/1‑149 (MQ=255)
       tggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTAt                                                                                                                                         <  2:870865/149‑1 (MQ=255)
        ggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                                        <  1:32214/149‑1 (MQ=255)
        ggAGAATGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                                        >  1:658012/1‑149 (MQ=255)
          agaTTGAGCCACCCTGGCACGTTCACCGTCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAg                                                                                                                                      <  1:166809/149‑1 (MQ=255)
            aTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCAGGc                                                                                                                                    <  1:219905/149‑1 (MQ=255)
            aTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCAGGc                                                                                                                                    <  1:867526/149‑1 (MQ=255)
            aTTGAGCCACCCTGGCACGATCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGc                                                                                                                                    <  1:179745/149‑1 (MQ=255)
                  ccACCCTGGCACGTTCACCGGCTTGCGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAAGCGTCTTATCAGGCCTACAg                                                                                                                              >  1:231068/1‑149 (MQ=255)
                     cccTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGAc                                                                                                                           <  2:944686/149‑1 (MQ=255)
                     cccTGGCACGGTCACCGGCTTGAGCGATCCATAGTCGCGAGGAACATCAAGTGCACGGAACATCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCGTTACGCGTCTTATCAGGCCTACAGGAc                                                                                                                           <  1:84145/149‑1 (MQ=39)
                      gcTGGCACGTTAACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGTATCAGCCTTCGATCTTACTTATTCAGAATATAATGATGCCGCGACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACt                                                                                                                          <  2:800046/148‑1 (MQ=39)
                      ccTGGCACGTTCACCGGCCTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGACCGGAACAGCCTTCGACCTTACATATTCAGAATATAATGATGCCGCCCCGACCGCCTGCTGCGACGCTGTACGCGTCTTATCAGGCCTACAGGACt                                                                                                                          <  2:74329/149‑1 (MQ=255)
                            aCGTTCACCGGCTTGAGCGTTCCATACTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATACTGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCAc                                                                                                                    <  1:432851/149‑1 (MQ=255)
                             cGTTCACCGGCGTGAGCGTTACATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCAcc                                                                                                                   <  2:58114/149‑1 (MQ=255)
                                   ccGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                             >  1:823627/1‑149 (MQ=255)
                                   ccGGCTTGAGCGCTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGACTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTaa                                                                                                             <  1:855689/149‑1 (MQ=255)
                                      gCTTGCGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGTATCATCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACTCATGATGCGACGCTTTACGAGTCTTATCAGGCCTACAGGAATGGTCAACTGGTAACCt                                                                                                          >  2:1002251/1‑149 (MQ=255)
                                        ttGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTAc                                                                                                        <  2:99460/149‑1 (MQ=255)
                                          gAGCGTTCCATAGTAGCGAGGAACATCAAGTGCACGGAACAGCCTTCGCTCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCCCCTGGTAACCTCCAt                                                                                                      >  1:768294/1‑149 (MQ=255)
                                          gAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAt                                                                                                      >  1:966830/1‑149 (MQ=255)
                                          gAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAt                                                                                                      >  1:976926/1‑149 (MQ=255)
                                           aGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtc                                                                                                     >  1:532943/1‑149 (MQ=255)
                                           aGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtc                                                                                                     >  1:137768/1‑149 (MQ=255)
                                               ttCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttct                                                                                                 <  1:335412/149‑1 (MQ=255)
                                                 ccATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGAGGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTcc                                                                                               <  1:367440/149‑1 (MQ=255)
                                                  cATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCa                                                                                              <  2:208662/149‑1 (MQ=255)
                                                    tagtagCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGc                                                                                            <  2:617841/149‑1 (MQ=255)
                                                          cGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGACTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtc                                                                                      <  1:425078/149‑1 (MQ=255)
                                                            aGGAACATCAAGTGCACGGAACAGCCTTCGCTCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCGTATCAGGCCTACAGGACTGGGCACCTGGTCACCTACATCTTCGCCAGCAActtctt                                                                                    >  2:199269/1‑149 (MQ=255)
                                                             ggAACATCAAGTGAACGGAACAGCCTTCGATATTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACttcttc                                                                                   <  2:102451/149‑1 (MQ=255)
                                                               aaCATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTcac                                                                                 <  2:887936/149‑1 (MQ=255)
                                                                aCATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTcaca                                                                                >  1:73772/1‑149 (MQ=255)
                                                                aCATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTcaca                                                                                <  1:61211/149‑1 (MQ=255)
                                                                 cATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACAt                                                                               >  1:532349/1‑149 (MQ=255)
                                                                      aGTTAACGGAACATCCTTCGATCTTAAATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTg                                                                          <  2:245470/149‑1 (MQ=255)
                                                                      aGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTg                                                                          >  1:1009719/1‑149 (MQ=255)
                                                                       ttGAACGGAACAGCCTTCTATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTTGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGc                                                                         <  1:515923/148‑1 (MQ=255)
                                                                        tGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGcc                                                                        <  2:411548/149‑1 (MQ=255)
                                                                           aCGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCAtt                                                                     >  1:635343/1‑149 (MQ=255)
                                                                            cGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGGCATTa                                                                    >  1:1023275/1‑149 (MQ=255)
                                                                             ggAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTAc                                                                   <  2:539585/149‑1 (MQ=255)
                                                                                    ccTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCGTCACATCTTTGCACTTACGGGAATc                                                            >  1:687415/1‑149 (MQ=255)
                                                                                     cTTCGATCTAACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCt                                                           <  2:612799/149‑1 (MQ=255)
                                                                                      ttCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt                                                          <  1:585779/149‑1 (MQ=255)
                                                                                        cGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTa                                                        <  1:245628/149‑1 (MQ=255)
                                                                                          aTCTTATATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTAtt                                                      <  2:777430/149‑1 (MQ=255)
                                                                                            cTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTACCCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTGATTcc                                                    >  2:615928/1‑149 (MQ=255)
                                                                                            cTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTcc                                                    <  2:208741/149‑1 (MQ=255)
                                                                                                    ttCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGc                                                                         <  1:891635/120‑1 (MQ=255)
                                                                                                    ttCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGc                                                                         >  2:891635/1‑120 (MQ=255)
                                                                                                     tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGccc                                           <  2:532943/149‑1 (MQ=255)
                                                                                                      cAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTAGAGGACTGGTCACCTGGTAACCAACATCATCTCCAGAAACAACTTCAAATCTTTGCCATTACTGGAATCTTTATTCCCCTCggacaa                                          >  1:147496/1‑145 (MQ=255)
                                                                                                      cAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCa                                          <  1:705824/149‑1 (MQ=255)
                                                                                                        gAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGt                                        >  1:91149/1‑149 (MQ=255)
                                                                                                               atgatgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAg                                 <  1:442207/149‑1 (MQ=255)
                                                                                                                 gatgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGc                               <  2:823627/149‑1 (MQ=255)
                                                                                                                      cgccACGAACGCCTGATGCTACTCTTTACGCGTCATATCAGTCCTACAGGACTGGTCACCTGGTAAACTACAACGTCTCCAGCAACTTCTTCACATCTTAGCCATTACGGGAATCTGTATTCCTCTCGGCCCAGTCATACAGGCGAcct                          >  2:234479/1‑147 (MQ=255)
                                                                                                                         cACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTc                       <  2:556437/149‑1 (MQ=255)
                                                                                                                          aCGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCg                      <  2:532349/149‑1 (MQ=255)
                                                                                                                          aCGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTGTCg                      >  2:485771/1‑149 (MQ=255)
                                                                                                                               ggCCTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTAGATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTGTCGCTTCg                 >  1:47670/2‑149 (MQ=255)
                                                                                                                                   tGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCt             <  1:370475/149‑1 (MQ=255)
                                                                                                                                     ataagACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCgtgt           >  1:283966/5‑149 (MQ=255)
                                                                                                                                          aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCGCGGCCCAGTCATTCAGGCGCCGTTCCGCTGCGTGCTGTAGATg      >  1:80066/1‑149 (MQ=255)
                                                                                                                                             gtgtACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCCGGCGACGTTTCGCTTCGTCCTGTAGATGttt   <  1:199269/146‑1 (MQ=255)
                                                                                                                                              tttACGCGTCTTATCAGGCCTACAGGACTGGTCACCAGCTAACCTACATCTTCTCCAGCAACATCTTCACATCTTTGCCATTACGGGAAACTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTg  >  1:161132/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTG  >  NC_000913/2139574‑2139864

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