Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,139,726 T→G 7.3% intergenic (+241/+33) 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,7260TG 7.3% 163.2 / 4.1 69intergenic (+241/+33)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 T (24/39);  new base G (3/2);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

GGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAAC  >  NC_000913/2139579‑2139873
                                                                                                                                                   |                                                                                                                                                   
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)
                                                                                                                                              gcgTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAg      >  1:486841/1‑149 (MQ=255)
                                                                                                                                                  cTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTCCGGGAATCTTTAATCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAAc  <  2:97258/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTGCGCAGTAAC  >  NC_000913/2139579‑2139873

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: