Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,342,193 C→T 5.3% E659E (GAG→GAA yfaL ← putative autotransporter adhesin YfaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,342,1930CT 5.3% 242.4 / 2.0 75E659E (GAG→GAAyfaLputative autotransporter adhesin YfaL
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base T (2/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CACCAACCCACAGTGACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAA  >  NC_000913/2342048‑2342341
                                                                                                                                                 |                                                                                                                                                    
cACCAACCCACAGTGACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCAt                                                                                                                                                  >  1:662950/1‑150 (MQ=255)
        cacaGTGACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGAtgc                                                                                                                                          <  2:968390/150‑1 (MQ=255)
        cacaGTGACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGAtgc                                                                                                                                          <  1:146237/150‑1 (MQ=255)
         acaGTGACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGAtgct                                                                                                                                         >  1:15767/1‑150 (MQ=255)
              gACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCg                                                                                                                                    >  1:559276/1‑150 (MQ=255)
              gACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAATCAATGTGCTACCCTCATCCTGATGCTGCCCg                                                                                                                                    >  2:336132/1‑150 (MQ=255)
                    gCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGctgctg                                                                                                                              >  1:73236/1‑150 (MQ=255)
                     cATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGt                                                                                                                             <  1:1174712/150‑1 (MQ=255)
                      aTAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTc                                                                                                                            >  2:537737/1‑150 (MQ=255)
                      aTAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTc                                                                                                                            >  2:1217221/1‑150 (MQ=255)
                          ggaaggaTATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGcc                                                                                                                        <  1:30534/150‑1 (MQ=255)
                           gaaggaTATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACACTCATCCTGATGCTGCCCCCTGCTGTCAGCCa                                                                                                                       >  2:329068/1‑150 (MQ=255)
                                   tCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTgcg                                                                                                               >  1:672630/1‑150 (MQ=255)
                                    cAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTgcgc                                                                                                              >  1:1398506/1‑150 (MQ=255)
                                      gCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGccc                                                                                                            >  2:506023/1‑150 (MQ=255)
                                       cAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGcccc                                                                                                           <  2:550449/150‑1 (MQ=255)
                                           aTCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTg                                                                                                       >  2:1163771/1‑150 (MQ=255)
                                            tCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGc                                                                                                      >  2:157855/1‑150 (MQ=255)
                                               ccTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCgtg                                                                                                   >  1:1087592/1‑150 (MQ=255)
                                               ccTTTCAGGGTGCCTTCTTCAACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCgtg                                                                                                   <  1:459216/150‑1 (MQ=255)
                                                   tCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTa                                                                                               <  1:873653/150‑1 (MQ=255)
                                                          gCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGc                                                                                        >  2:540751/1‑150 (MQ=255)
                                                          gCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGc                                                                                        >  2:425462/1‑150 (MQ=255)
                                                          gCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGc                                                                                        >  1:66376/1‑150 (MQ=255)
                                                           ccTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCa                                                                                       >  1:709206/1‑150 (MQ=255)
                                                             ttcttcGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATc                                                                                     >  2:942828/1‑150 (MQ=255)
                                                              tcttcGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCa                                                                                    <  2:29705/150‑1 (MQ=255)
                                                               cttcGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCCCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCaa                                                                                   >  1:555593/1‑150 (MQ=255)
                                                                 tcGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGctgctg                                                                                                                              >  2:759155/1‑105 (MQ=255)
                                                                 tcGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGctgctg                                                                                                                              <  1:759155/105‑1 (MQ=255)
                                                                 tcGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACt                                                                                 >  2:1299570/1‑150 (MQ=255)
                                                                  cGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTg                                                                                <  2:15767/150‑1 (MQ=255)
                                                                        gTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACtt                                                                          >  1:204075/1‑150 (MQ=255)
                                                                        gTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAATCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACtt                                                                          <  1:336132/150‑1 (MQ=255)
                                                                         tACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTc                                                                         >  1:151406/1‑150 (MQ=255)
                                                                          aCCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCa                                                                        >  1:748288/1‑150 (MQ=255)
                                                                           ccgccgACTGCGGTGTTCCGCTGGCGGTCAGCTGCAGTGTACCAGCCCCCGTTTTAGTCAATGGGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGGCATCATTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGACACTTCAc                                                                       <  2:1419182/150‑1 (MQ=255)
                                                                             gccgACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCg                                                                     <  1:1217221/150‑1 (MQ=255)
                                                                             gccgACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCg                                                                     >  2:592521/1‑150 (MQ=255)
                                                                             gccgACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGAGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCg                                                                     >  1:1338917/1‑150 (MQ=255)
                                                                                gACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCTTct                                                                  <  2:151406/150‑2 (MQ=255)
                                                                                gACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCg                                                                  <  2:559276/150‑1 (MQ=255)
                                                                                   tGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGgcg                                                               <  1:537737/150‑1 (MQ=255)
                                                                                     cGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGgcgcg                                                             >  1:1292721/1‑150 (MQ=255)
                                                                                      gTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCTTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGgcgcgc                                                            >  1:556337/1‑150 (MQ=255)
                                                                                      gTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCTTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGgcgcgc                                                            >  1:560906/1‑150 (MQ=255)
                                                                                           aCCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACGCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTTCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTc                                                       <  1:185382/150‑1 (MQ=255)
                                                                                           aCCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTc                                                       <  1:468331/150‑1 (MQ=255)
                                                                                            ccGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCa                                                      <  1:157855/150‑1 (MQ=255)
                                                                                              gCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAat                                                    >  2:199185/1‑150 (MQ=255)
                                                                                               cTGGCGGTCAGCTCCAGGGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGCTGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAata                                                   <  1:687262/150‑1 (MQ=255)
                                                                                                 ggCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATc                                                 >  1:305897/1‑150 (MQ=255)
                                                                                                  gCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCa                                                <  2:73236/150‑1 (MQ=255)
                                                                                                   cGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCAc                                               >  1:557955/1‑150 (MQ=255)
                                                                                                            ccAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGc                                      >  1:180124/1‑150 (MQ=255)
                                                                                                              aGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCt                                    >  1:617070/1‑150 (MQ=255)
                                                                                                                 gtACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGcc                                 >  2:1091869/1‑150 (MQ=255)
                                                                                                                 gtACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGcc                                 <  2:672630/150‑1 (MQ=255)
                                                                                                                  tACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCg                                <  1:506023/150‑1 (MQ=255)
                                                                                                                  tACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCg                                <  2:147861/150‑1 (MQ=255)
                                                                                                                      aGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATc                            <  2:66376/150‑1 (MQ=255)
                                                                                                                      aGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACAGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCaga                            >  2:1329707/1‑148 (MQ=255)
                                                                                                                        cccccGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTg                          <  2:1398506/150‑1 (MQ=255)
                                                                                                                           ccGTTTTAGTCAATGAGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGACGAATCTGaat                       <  2:1338917/150‑2 (MQ=255)
                                                                                                                              ttttAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAAcgcg                    <  1:199185/150‑1 (MQ=255)
                                                                                                                               tttAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCgt                   >  2:666826/1‑150 (MQ=255)
                                                                                                                                   gTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCa               <  1:879135/150‑1 (MQ=255)
                                                                                                                                   gTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCa               <  1:540751/150‑1 (MQ=255)
                                                                                                                                        tgtgCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTc          >  2:456610/1‑150 (MQ=255)
                                                                                                                                        tgtgCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTc          <  2:1292721/150‑1 (MQ=255)
                                                                                                                                          tgCTACCTTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCgg        <  2:556337/150‑1 (MQ=255)
                                                                                                                                          tgCTACCTTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCgg        <  2:560906/150‑1 (MQ=255)
                                                                                                                                          tgCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCgg        <  2:748288/150‑1 (MQ=255)
                                                                                                                                          tgCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCgg        >  1:475715/1‑150 (MQ=255)
                                                                                                                                           gttACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGa       >  2:1152455/3‑150 (MQ=255)
                                                                                                                                                ccTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTaa  <  1:1163771/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CACCAACCCACAGTGACGAAGCATAAGGAAGGATATCAGCAACATCACCTTTCAGGGTGCCTTCTTCGACACGTACCGCCGACTGCGTTGTACCGCTGGCGGTCAGCTCCAGTGTACCAGCCCCCGTTTTAGTCAATGTGCTACCCTCATCCTGATGCTGCCCGCTGCTGTCAGCCATCAGTGCGCCCCACTGCGTGTCTACCCCCGCATCAACTGCCACTTCACCGTCGGCGCGCATTTCAATATCACGTCCGTGGCCTGCCGCATCTGAACGCGTGTCAACGTCGGAGGTAA  >  NC_000913/2342048‑2342341

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