Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,140,400 G→A 6.4% R405W (CGG→TGG)  asmA ← 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,140,4000GA 6.4% 344.8 / 6.9 110R405W (CGG→TGG) asmAsuppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (56/46);  new base A (4/3);  total (61/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 111 total aligned reads displayed.

GACATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTC  >  NC_000913/2140251‑2140548
                                                                                                                                                     |                                                                                                                                                    
gACATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCg                                                                                                                                                      >  2:1439923/1‑150 (MQ=255)
  cATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTg                                                                                                                                                    >  1:411691/1‑150 (MQ=255)
   ctgtgCGTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTgg                                                                                                                                                   >  1:1147528/2‑150 (MQ=255)
      ggCCTGTCCTTGCCAGTTGTGGCGGAATGCGGCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTg                                                                                                                                                <  2:1098266/149‑1 (MQ=255)
       gCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTgg                                                                                                                                               <  2:155159/150‑1 (MQ=255)
             ccTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGttt                                                                                                                                         <  2:1566545/150‑1 (MQ=255)
              cTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTa                                                                                                                                        <  1:1597969/150‑1 (MQ=255)
              cTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTGCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTGGaaaacgt                                                                                                                                        <  1:1706084/150‑6 (MQ=255)
                tGCCAGTTGTGGCGGAATGCGGCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCACGGCTGGaaaacgttg                                                                                                                                      <  2:1421879/150‑8 (MQ=255)
                   cAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCcgc                                                                                                                                   >  1:1632545/1‑150 (MQ=255)
                   cAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCcgc                                                                                                                                   >  1:1231607/1‑150 (MQ=255)
                    aGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCcgcg                                                                                                                                  >  2:706901/1‑150 (MQ=255)
                        gtgGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGAtt                                                                                                                              <  2:346511/150‑1 (MQ=255)
                           gCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTACGGTCAACTAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAt                                                                                                                           >  2:1223062/1‑150 (MQ=255)
                              gAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGa                                                                                                                        >  2:170660/1‑150 (MQ=255)
                               aaTGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGAt                                                                                                                       <  1:1695099/150‑1 (MQ=255)
                                aTGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATg                                                                                                                      >  1:346628/1‑150 (MQ=255)
                                aTGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATg                                                                                                                      <  1:584228/150‑1 (MQ=255)
                                    gTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcgc                                                                                                                  >  1:26169/1‑150 (MQ=255)
                                          aTCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGTTAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAg                                                                                                            <  1:521257/150‑1 (MQ=255)
                                           tCTATGTCAGCCCCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGc                                                                                                           >  1:1045263/1‑150 (MQ=255)
                                            cTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCCGCCGTGGCGGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCg                                                                                                          >  1:1055869/1‑150 (MQ=255)
                                                gTCAGCACCGGAGAAATCACCAGCCAGAGACATATTTCCGGTCAACGCAATCCGAGAGTTAAACGCCTTCAGAATGGTACCACTCTCAACGATTTCCAGCCTTGGCAGAAAGTTTAACCCCTCATTTATTGATGTCGCGTCCAGCGAGcc                                                                                                      >  2:1610495/1‑150 (MQ=255)
                                                 tCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGccc                                                                                                     <  1:539574/150‑1 (MQ=255)
                                                  cAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCg                                                                                                    >  1:1639996/1‑150 (MQ=255)
                                                  cAGCACCGGAGAAATCACCAACCAGTGACATTTTTCAGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCg                                                                                                    >  1:1483956/1‑150 (MQ=255)
                                                   cgCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTACACGCCTTCCGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCgg                                                                                                   <  2:849794/149‑1 (MQ=255)
                                                    gCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGACATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGc                                                                                                  <  1:1039093/150‑1 (MQ=255)
                                                    gCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTTAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGc                                                                                                  <  2:1420406/150‑1 (MQ=255)
                                                             aaaTCATCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAAcc                                                                                         <  1:1009025/150‑1 (MQ=255)
                                                                    cagccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAc                                                                                  <  2:1125865/150‑1 (MQ=255)
                                                                    cagccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAc                                                                                  >  1:356497/1‑150 (MQ=255)
                                                                    cagccagTGACATTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGGGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAcc                                                                                 <  2:1128056/150‑1 (MQ=255)
                                                                     agccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAcc                                                                                 <  1:663445/150‑1 (MQ=255)
                                                                      gccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCg                                                                                <  1:1434745/150‑1 (MQ=255)
                                                                        cagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGtt                                                                              >  2:586150/1‑150 (MQ=255)
                                                                           tGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAg                                                                           <  1:1717795/150‑1 (MQ=255)
                                                                           tGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAg                                                                           <  1:502068/150‑1 (MQ=255)
                                                                           tGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAg                                                                           <  2:1231607/150‑1 (MQ=255)
                                                                           tGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGAGTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAg                                                                           >  2:629471/1‑150 (MQ=255)
                                                                                      cGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCa                                                                >  1:1196760/1‑150 (MQ=255)
                                                                                      cGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCa                                                                >  2:30167/1‑150 (MQ=255)
                                                                                       ggTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAg                                                               >  1:1292468/1‑150 (MQ=255)
                                                                                          cAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGTAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCCTGCCCGGCAGTGAAACCTGTCCACCCTTAAGTTTGGCCTGCAGTTg                                                            >  1:103635/1‑150 (MQ=255)
                                                                                          cAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTg                                                            >  1:1413614/1‑150 (MQ=255)
                                                                                          cAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTg                                                            <  2:219305/150‑1 (MQ=255)
                                                                                            aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAg                                                          >  2:1053327/1‑150 (MQ=255)
                                                                                            aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAg                                                          <  2:26169/150‑1 (MQ=255)
                                                                                              gAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCGTGCAGTTGAGTa                                                        >  2:752985/1‑150 (MQ=255)
                                                                                                aaTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAAt                                                      <  1:170660/150‑1 (MQ=255)
                                                                                                aaTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAAt                                                      <  1:590263/150‑1 (MQ=255)
                                                                                                aaTCGGATACTTAAACGACTTCAGAATTGTACCAATCGCAACGATTTCCAGCCGTGGATGGAAGTTGATCCTCGGATTTATTTATGTCGCGTCCATCGTGCCCGTCAGTGAAACCTGGCCACCGTTAAGGTTGCCCTGCAGTTGAGTAAt                                                      <  1:908258/150‑1 (MQ=255)
                                                                                                 aTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAAtt                                                     >  2:1400402/1‑150 (MQ=255)
                                                                                                    ggATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTcc                                                  >  2:1608333/1‑150 (MQ=255)
                                                                                                      aTAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAg                                                <  1:1439923/150‑1 (MQ=255)
                                                                                                       tAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGc                                               >  1:279286/1‑150 (MQ=255)
                                                                                                        aGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCa                                              >  2:399652/1‑150 (MQ=255)
                                                                                                         gTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCaa                                             >  2:1703527/1‑150 (MQ=255)
                                                                                                         gTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCaa                                             <  2:1529270/150‑1 (MQ=255)
                                                                                                         gTTAAACGCCTTCAGAATGGTACCAATCTAAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAAGTTCCCGCaa                                             >  1:407560/1‑150 (MQ=255)
                                                                                                           tAAACGCCTTCAGAATGGTCCCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACAGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAc                                           >  1:560200/1‑150 (MQ=255)
                                                                                                           tAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGTGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAc                                           >  2:21386/1‑150 (MQ=255)
                                                                                                             aaCGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAccc                                         >  1:1586614/1‑150 (MQ=255)
                                                                                                             aaCGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATATATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAccc                                         >  1:242749/1‑150 (MQ=255)
                                                                                                               cGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAACCCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCgc                                       >  2:1176922/1‑149 (MQ=255)
                                                                                                                 ccTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACt                                     <  2:1639996/150‑1 (MQ=255)
                                                                                                                 ccTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACt                                     <  2:1222552/150‑1 (MQ=255)
                                                                                                                  cTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGAAtt                                    >  2:28345/1‑150 (MQ=255)
                                                                                                                    tCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgt                                  >  1:1005012/1‑150 (MQ=255)
                                                                                                                     cAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                 <  2:1196760/150‑1 (MQ=255)
                                                                                                                         aTGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTc                             <  1:30167/150‑1 (MQ=255)
                                                                                                                           ggTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAt                           >  1:508255/1‑150 (MQ=255)
                                                                                                                           ggTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAt                           >  2:1644592/1‑150 (MQ=255)
                                                                                                                            gTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAtt                          >  1:478474/1‑150 (MQ=255)
                                                                                                                            gTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTATCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCAtt                          <  1:1344620/150‑1 (MQ=255)
                                                                                                                             tACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCGCTGCAGTTTAGTATTTTCCAGCAAAACCGACTTGTTGGTCAttt                         >  2:837598/1‑150 (MQ=255)
                                                                                                                              aCCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTg                        >  1:632550/1‑150 (MQ=255)
                                                                                                                                cAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCg                      <  1:1560323/150‑1 (MQ=255)
                                                                                                                                cAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCg                      <  2:1147528/150‑1 (MQ=255)
                                                                                                                                  aTCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTg                    <  1:706901/150‑1 (MQ=255)
                                                                                                                                  aTCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTg                    <  1:806246/150‑1 (MQ=255)
                                                                                                                                  aTCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTg                    <  2:1045263/150‑1 (MQ=255)
                                                                                                                                   tctcAACGTTTTCCATCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTgg                   <  2:1487072/150‑1 (MQ=255)
                                                                                                                                    ctcAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGc                  <  1:1223062/150‑1 (MQ=255)
                                                                                                                                     tcAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                 >  1:1029295/1‑150 (MQ=255)
                                                                                                                                     tcAACGTTTTCCAGCCGGGGCTGGAAGTGTATCCGCGGATTTCTTGCTGTCGCGGCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCTGCAAACCCGACTTGTTGGTCATTTGCGTGGct                 <  1:1610495/150‑2 (MQ=255)
                                                                                                                                        aCGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGGTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACa              >  1:624032/1‑150 (MQ=255)
                                                                                                                                         aaacttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACAt             >  2:526335/5‑150 (MQ=255)
                                                                                                                                          ttgttCCCGCCGGGGCTGGAAGTTTATCCGGGGATTTATTGTTGTAGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACAGTTACGTTTGCCCTGCAGTTGAGTAATTTCCAGCTCACCCGACTTGTTGGTCATTTGCGGGGCAACTTc            <  2:103635/147‑1 (MQ=255)
                                                                                                                                          aacttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATc            >  1:747683/4‑150 (MQ=255)
                                                                                                                                            cttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTg          >  1:441233/2‑150 (MQ=255)
                                                                                                                                            cttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTg          >  1:63973/2‑150 (MQ=255)
                                                                                                                                              tCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTa        >  2:1029799/1‑150 (MQ=255)
                                                                                                                                               ccAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTCAAACCTGTCAACCGTTAAGTTTGCCCTGCAGTTGAGTAAATTCCAGCAAACCCGAATTGTTGGTCAATTGCGTGGCAACATCTGTaa       >  1:867024/1‑150 (MQ=255)
                                                                                                                                               ccAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGAGCCCGGCAGTGAACCCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGCCTTGTTGGTCATTTGCGTGGCTACATCTGTaa       >  2:739914/1‑150 (MQ=255)
                                                                                                                                                 aGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTaaaa     >  2:1382007/1‑150 (MQ=255)
                                                                                                                                                  gCCGTGGCTGGAAGTTTATCCGCTGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAat    <  2:485013/150‑1 (MQ=255)
                                                                                                                                                  gCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGCAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCACCATCTGTAAAat    >  1:504284/1‑150 (MQ=255)
                                                                                                                                                   ccGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAatt   <  2:1488394/150‑1 (MQ=255)
                                                                                                                                                    cGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAattc  >  1:1162153/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACATGTGCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTC  >  NC_000913/2140251‑2140548

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