Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,139,716 T→G 6.3% intergenic (+231/+43) 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,7160TG 6.3% 312.0 / 10.8 111intergenic (+231/+43)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 (47/55);  new base G (5/2);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

Only 100 of 112 total aligned reads displayed.

AGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTG  >  NC_000913/2139574‑2139864
                                                                                                                                              |                                                                                                                                                    
aGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcgcg                                                                                                                                               <  2:1289600/150‑1 (MQ=255)
aGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcgcg                                                                                                                                               <  1:186958/150‑1 (MQ=255)
aGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAcgcg                                                                                                                                               <  1:1166520/150‑1 (MQ=255)
 gggTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGt                                                                                                                                              <  1:502762/150‑1 (MQ=255)
   gtggtgGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCGTAAAGCGTcg                                                                                                                                            <  1:1039209/150‑2 (MQ=255)
     ggtggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTa                                                                                                                                          >  1:74579/1‑150 (MQ=255)
       tggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                                        >  2:807144/1‑150 (MQ=255)
        gtagaTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCa                                                                                                                                       <  2:1577119/148‑1 (MQ=255)
        ggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCa                                                                                                                                       <  2:1654396/150‑1 (MQ=255)
            aTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGcc                                                                                                                                   <  2:1083215/150‑1 (MQ=255)
            aTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGcc                                                                                                                                   <  1:414760/150‑1 (MQ=255)
                 gCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACCGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAg                                                                                                                              <  2:862983/150‑1 (MQ=255)
                   cACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACTTATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGa                                                                                                                            >  2:1138075/1‑150 (MQ=255)
                       cTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTgg                                                                                                                        <  1:279602/150‑1 (MQ=255)
                           cACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACTCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCAc                                                                                                                    <  1:1533937/150‑1 (MQ=255)
                             cGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCt                                                                                                                  <  1:1439583/150‑1 (MQ=255)
                              gTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTg                                                                                                                 <  1:1057972/150‑1 (MQ=255)
                              gTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTg                                                                                                                 <  2:1527120/150‑1 (MQ=255)
                                tCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGt                                                                                                               <  2:470116/150‑1 (MQ=255)
                                    cGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAAcc                                                                                                           <  1:1621953/150‑1 (MQ=255)
                                    cGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAAcc                                                                                                           >  1:496850/1‑150 (MQ=255)
                                     ggCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTGCGATCGTACATATTCAGAATATAATGATGCCGCCACGCACGCCTGATGCGCCGCTTTACGCGTCTTATCCGGCCTACACGACTGGTCACCTGGTAAGCt                                                                                                          >  2:448184/1‑150 (MQ=255)
                                      gCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTa                                                                                                         <  1:517623/150‑1 (MQ=255)
                                        ttGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCAGGTAACCTACa                                                                                                       <  2:1184877/150‑1 (MQ=255)
                                         tGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGCTCTGACATATTCAGAATATAATGCTGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACag                                                                                                      >  1:234709/1‑149 (MQ=255)
                                            gCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTCACCTACAtctt                                                                                                   >  1:1621475/1‑150 (MQ=255)
                                            gCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAAACTACAtctt                                                                                                   >  1:1615075/1‑150 (MQ=255)
                                             cGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCAACAGGACTGGTCACCTGGTAACCTACAtcttc                                                                                                  <  2:1284433/150‑1 (MQ=255)
                                                 ccATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACGTCTTCTCCa                                                                                              >  1:242970/1‑150 (MQ=255)
                                                  cATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAg                                                                                             <  1:531044/150‑1 (MQ=255)
                                                     agtagCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCaa                                                                                          <  2:717083/150‑1 (MQ=255)
                                                     agtagCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTAATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAACTTCTCAAGCaa                                                                                          >  2:149455/1‑150 (MQ=255)
                                                     agtagCGAGGAACATCAAGTGAACGGAACAGCATTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCAGGTAACCTACATCATCTCCAGCaa                                                                                          >  2:1476565/1‑150 (MQ=255)
                                                       tagCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAct                                                                                        <  1:1718072/150‑1 (MQ=255)
                                                          cGAGGAACATCAAGTGAACGGAAGAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtct                                                                                     <  2:114715/150‑1 (MQ=255)
                                                          cGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtct                                                                                     <  2:437417/150‑1 (MQ=255)
                                                           gAGGAACATCAAGTGAACGGACCAGCCTTCGATCTTACATATTCAGAATATAATGATCCCGCCACGAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACAGGCCTGGTCACCTGGAAACCTACATCTGCTCCAGCAActtctt                                                                                    >  1:177982/1‑150 (MQ=39)
                                                           gAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtctt                                                                                    <  1:536262/150‑1 (MQ=255)
                                                            aGGAACATCAAGTGAACGGAACAGCCTTCGCTCTTACATTTTCAGAATATAATGATGCCGCCACGCACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACCGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACttcttc                                                                                   <  2:273600/150‑1 (MQ=255)
                                                            aGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGAAGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACttcttc                                                                                   <  1:326095/150‑1 (MQ=255)
                                                             ggAACATCAAGTGAACGGGACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTca                                                                                  >  2:92224/1‑150 (MQ=255)
                                                              gAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTcac                                                                                 <  2:1174247/150‑1 (MQ=255)
                                                                 cATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTGTACGCGTCTTATCAGGCCTACAGGACTGGTCCCCTGGTAACCTACATCTTCTCCAGCAACTACTTCACAta                                                                              >  2:351088/1‑149 (MQ=255)
                                                                       gTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGGAACTTCTTCACATCTTTGcc                                                                        <  2:512381/150‑1 (MQ=255)
                                                                       gTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGcc                                                                        >  1:692475/1‑150 (MQ=255)
                                                                       gTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGcc                                                                        <  1:1347245/150‑1 (MQ=255)
                                                                       gTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCATTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGcc                                                                        <  1:1321521/150‑1 (MQ=255)
                                                                        tGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCa                                                                       <  1:140935/150‑1 (MQ=255)
                                                                          aaCGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCAtt                                                                     <  2:479426/150‑1 (MQ=255)
                                                                          aaCGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCAtt                                                                     >  1:1624360/1‑150 (MQ=255)
                                                                              gAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACgg                                                                 <  1:1174720/150‑1 (MQ=255)
                                                                              gAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACgg                                                                 >  1:1346674/1‑150 (MQ=255)
                                                                                   gCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATc                                                            >  2:1597375/1‑150 (MQ=255)
                                                                                   gCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATc                                                            <  2:74579/150‑1 (MQ=255)
                                                                                     cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt                                                          <  1:616593/150‑1 (MQ=255)
                                                                                     cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt                                                          <  1:198448/150‑1 (MQ=255)
                                                                                     cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt                                                          <  1:239047/150‑1 (MQ=255)
                                                                                     cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTAATCAGGCATACAGGACTGGTCAACTGGTAACCTACATCTTCTACAGCAACTTCTTCAAATCTTTGCCATTACGGGACTCtt                                                          >  2:866087/1‑150 (MQ=255)
                                                                                      ttCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCttt                                                         <  2:535145/150‑1 (MQ=255)
                                                                                      ttCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCttt                                                         >  1:62934/1‑150 (MQ=255)
                                                                                       tCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTa                                                        >  1:113608/1‑150 (MQ=255)
                                                                                         gagcTTACATATTAAGAATATAATGATGCCGCCTCGAACGACGTATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTAtt                                                      <  2:869183/147‑1 (MQ=255)
                                                                                          aTCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTc                                                     <  2:1316491/150‑1 (MQ=255)
                                                                                              tACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCACt                                                 >  1:1592293/1‑150 (MQ=255)
                                                                                                   atTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGcc                                            <  2:412841/150‑1 (MQ=255)
                                                                                                    ttCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGccc                                           >  2:928938/1‑150 (MQ=255)
                                                                                                      cAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAg                                         >  2:1499996/1‑150 (MQ=255)
                                                                                                      cAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAg                                         <  1:1088534/150‑1 (MQ=255)
                                                                                                       aGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGt                                        >  1:381111/1‑150 (MQ=255)
                                                                                                        gAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGCCTGGTCACCTGGTACCCTACATCTTCTCCAGCACCTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTc                                       >  1:1241148/1‑150 (MQ=255)
                                                                                                        gAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTc                                       <  2:309749/150‑1 (MQ=255)
                                                                                                         aaTATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCa                                      <  2:113608/150‑1 (MQ=255)
                                                                                                          atatAATGATGCCGCCACGAACGCCTGGTGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCAt                                     >  1:212590/1‑150 (MQ=255)
                                                                                                            ataATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTc                                   <  1:807144/150‑1 (MQ=255)
                                                                                                             taATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCa                                  >  1:358366/1‑150 (MQ=255)
                                                                                                              aaTGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAg                                 >  2:706251/1‑150 (MQ=255)
                                                                                                              aaTGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAg                                 >  2:399663/1‑150 (MQ=255)
                                                                                                                tgatgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGc                               >  2:641333/1‑150 (MQ=255)
                                                                                                                 gatgGCGCCACGAAGGGCTGAGGGGACGCATTACGCGTCTTATCAGGCCTACAGGACTGGTCACGTGGTAACCGTCATCTTCTCCAGCATCTTCTTCACATCTTTGCCATTACGGGAATCTTTTTTCCGCTCGGCGCAGTCATTCAGGCg                              <  2:1579304/150‑1 (MQ=255)
                                                                                                                    gccgccACGAACGCCTGCTGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACCTCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGCCg                           >  1:300552/1‑150 (MQ=255)
                                                                                                                        ccACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTc                       <  1:1138075/150‑1 (MQ=255)
                                                                                                                              aCGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGTAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTCTCGCTTCg                 <  1:997125/150‑1 (MQ=255)
                                                                                                                              aCGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCg                 >  1:784287/1‑150 (MQ=255)
                                                                                                                               ggCCTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGt                >  2:540456/2‑150 (MQ=255)
                                                                                                                                gCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTc               >  1:1538431/1‑150 (MQ=255)
                                                                                                                                 ccTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTGCAGCAACTGCTTGACATCTTTGCCATTACGGGAATCTTTATTCAGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTcc              >  1:919506/1‑150 (MQ=255)
                                                                                                                                 ccTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTcc              >  2:1180857/1‑150 (MQ=255)
                                                                                                                                 ccTGATGCGACGCTTTACGCGTCTTATCAGGCATACAGGACTGGTCAGCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCATTGCCATTACGGGATTCTGTATTCCGCTCGGACCAGTCATTCAGGCGACGTTTCGcgaagtcc              >  1:264174/1‑143 (MQ=255)
                                                                                                                                 ccTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTcc              >  1:562721/1‑150 (MQ=255)
                                                                                                                                 ccTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGGTTCGCTTCGTcc              >  1:1319697/1‑150 (MQ=255)
                                                                                                                                  cTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGCAATCTTTATTCCACTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCt             >  1:215895/1‑150 (MQ=255)
                                                                                                                                   tGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTg            <  1:1021197/150‑1 (MQ=255)
                                                                                                                                     aTGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTa          <  1:1449483/150‑1 (MQ=255)
                                                                                                                                      taagACGCGTAAAGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACTTCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAg         >  2:678983/4‑150 (MQ=255)
                                                                                                                                      taagACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAg         >  2:624427/4‑150 (MQ=255)
                                                                                                                                      tGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAg         <  2:927678/150‑1 (MQ=255)
                                                                                                                                       gCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCATCAACTTCTTCACATCTTTGGCATTACGGGAATCTTTATTCAGCTCGGCCCAGTCATTCAGGCGAAGTTTCGCTTCGTCCTGTAGa        >  2:1100312/1‑150 (MQ=255)
                                                                                                                                           cGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGtt    <  2:1180961/150‑1 (MQ=255)
                                                                                                                                             tttttCGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCTTTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGTTGTTTg  <  1:149455/149‑1 (MQ=255)
                                                                                                                                             cTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTGTAGATGTTTg  <  1:1624615/150‑1 (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 < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: