Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 397,917 C→T 12.3% L16L (CTG→TTG)  yaiW → microcin Bac7 uptake protein, outer membrane surface‑exposed lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913397,9170CT12.3% 153.5 / 3.6 69L16L (CTG→TTG) yaiWmicrocin Bac7 uptake protein, outer membrane surface‑exposed lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (25/35);  new base T (7/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 7.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCA  >  NC_000913/397683‑398132
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                             
cTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATCCCTTGAGCTAAAAAGGAGAAGTAATGTCGCCCGTCCATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTTGt                                                                                                                                                                                                                           >  2:506929/1‑239 (MQ=255)
cTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTTTggg                                                                                                                                                                                                                           >  2:204078/1‑238 (MQ=255)
       ttgtgaCGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTCACTGATTAACTCATGGAACACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg                                                                                                                                                                                                                    <  2:392735/236‑1 (MQ=255)
       tgttgaCGCAATTTACCACAGTTTTTGGGCAGGTTGGCGGTGCTTTCCAGTTCCTGAGTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg                                                                                                                                                                                                                    <  1:270625/236‑1 (MQ=255)
       tgattaCGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGGGCTTTCCAGTACCTGATTAACTCATGGACCACACTTGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGATCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg                                                                                                                                                                                                                    <  2:55406/239‑1 (MQ=255)
       tgatgaCGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg                                                                                                                                                                                                                    <  2:718342/239‑1 (MQ=255)
       tgatgaCGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg                                                                                                                                                                                                                    <  2:658190/239‑1 (MQ=255)
       tgaaggcggaGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCCTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCg                                                                                                                                                                                                                    <  2:107540/230‑1 (MQ=255)
         attcctattgttgccaccctttttGGTGAGTTTGTCGGTGTTTGACAGTCCCTGATTAACTCATGGACCACTCTGGATGGGGTGATGTCTATCTACAAACGTCCGCGCAGCTTTGAACATGAGCGGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAATAGGAGAAGTAATGCCGCGCGTCAATCCTTTATCTTCTCTGCCATTGTTGGCGGTTCTCGTTTTAgatgc                                                                                                                                                                                                                  <  1:597093/220‑5 (MQ=255)
         atgaCGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGGGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGATGTCATTGTTGGCGGTTCTGGGTTTTACggg                                                                                                                                                                                                                  >  2:369602/1‑239 (MQ=255)
          tgaCGCAGATTACCAACGTTTTTGGGCAGGTTCGCAGTGCTTTCCGGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCggg                                                                                                                                                                                                                  <  2:120872/238‑1 (MQ=255)
                    tgccATCGTTTTTGCTCAGGTTCGCGGTGCTTTCCTGTCCCTGATTCAATACGGGCACACAGCGGTTGAGTTGATGTCTACCTACAAGCGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGCAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTGTAGCGGGGTGTAGTAGTc                                                                                                                                                                                                       <  1:544028/237‑1 (MQ=255)
                     tccaccttttttGGTCAGTGTCGCGGGTCTTTCCAGTCCCTGGGTAACTCATGGGCCACAGTGGTTTAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCGTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCa                                                                                                                                                                                                      <  1:221482/232‑1 (MQ=255)
                            tttttGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGAACGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTGCCTGTCATTGTTGGCGGTGTCGGTTTTGAGGGGGgtgaggtttcacggccc                                                                                                                                                                                                 >  2:807431/1‑224 (MQ=255)
                             ttttGGGCAGGTTCGCGGTTCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGc                                                                                                                                                                                              <  1:651755/239‑1 (MQ=255)
                                       gTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAAAATGAGCTGGATGGTTACAAAATTCAGGAAGTAACCCATACCTTGGGCTAAACAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTTTGGTTTTGGCGGGGGGAGGAAGTGCGGCGCCGCAACCGTTGaa                                                                                                                                                                                 >  1:172606/1‑239 (MQ=255)
                                                 cTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGgtgt                                                                                                                                                                                                              <  2:372111/203‑1 (MQ=255)
                                                      cAGTTCCTGATTAACTCAGGGGCCCCGCTGGTTGAGTTGATGTCTATCTACAATCGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCGCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGc                                                                                                                                                                 <  2:781292/239‑1 (MQ=255)
                                                        ttccctttttAACTAATGGTCCACACGGGTTGAGTTGAAGTGTTTCTTAACACGGCTGCTTAGTTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTTAGCGAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGa                                                                                                                                                               <  1:883872/231‑1 (MQ=255)
                                                          tccTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGAt                                                                                                                                                              <  2:523394/237‑1 (MQ=255)
                                                          cccTTATTTACTCATGGTCCACTTTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTGGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATa                                                                                                                                                             <  1:84953/238‑1 (MQ=255)
                                                           ccTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTACAAAGGAGAAGTAATGTCGCGCGTCAATCCGCTATCTTCGCTGTCATTGTTGGCGGTTTTGGGTTTTGGGGGGGGTAGTGGTACGGGTCCGCACCCGTTGAAAAAAGGGGGAAAAGGGaaaa                                                                                                                                                            >  1:450694/1‑236 (MQ=255)
                                                           ccTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGGTTTGTAGGGGGGTGGGGGGCAGGCTCCGGAACCCTTGAAAAAGGAGGaaaaaaaggaag                                                                                                                                                            >  1:388695/1‑232 (MQ=255)
                                                              gATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTTTGGTTTTACAGGGGGGAAGGAGTCAGTCTCCCCTACCGTTTGAAAAGAGGGCAaaaaggaaaaaag                                                                                                                                                           >  1:143005/1‑230 (MQ=255)
                                                              gATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTTGCGGGGGTGAGGAGTCAGGGTCCGCAACCCTTGGAAAAAGGGGGaaaaggggtagaag                                                                                                                                                          >  2:129989/1‑229 (MQ=255)
                                                              gATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGCGGTAGTAGTCCGGCTCCGCAACCGTTGTAAAAAGGCGCAAAAGCGATAGATg                                                                                                                                                         >  2:800196/1‑239 (MQ=255)
                                                              gATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAACGGGGGGGAGTAGGTAAGCTCCCCAACATTTTGAAAAGGGGGGAAAACGcacgagag                                                                                                                                                          >  2:442260/1‑232 (MQ=255)
                                                                   aCTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATCGATGTagcg                                                                                                                                                    >  2:891605/1‑239 (MQ=255)
                                                                      cATGGACCACACTGGTTGAGCTGATGTCTATCTACAACCGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGACGTAATGTCGCGCGTACATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTACAAGGCAGGCTCCGCAACCGTTGGACAAGGGTGGAAAAGCGCTAGacagatcgag                                                                                                                                                  >  2:161224/1‑229 (MQ=255)
                                                                         ggACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCCTTGTTGGCGGTTCTGGTTTTAACGGGGTGTAGTAGTCAGGCCCCGCAACCGCTGCAAAAGGGCGAAAAAGCGATAGATGTAGCGAGcgtc                                                                                                                                              >  2:410449/1‑239 (MQ=255)
                                                                         ggACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGGTTTAGCGGGGGGTGGGAGTGCGGGCCCGCAACCGCTGGAAAAAGGCGaaaaaaagaaaaaagtagagaggggg                                                                                                                                              >  1:428408/1‑218 (MQ=255)
                                                                               cgtttgtttgtttgaGGTCTACATTCAAACGTCTGCGACGCGTTGAAAATGAGGTGGATTGTGACAAAATTCAGGAAGTAACCCAGCCCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCTGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCgg                                                                                                                                        <  1:92571/228‑1 (MQ=255)
                                                                               cacTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTTTGGTTTTAGCGGGGGGTAGAAGGCAGGCTCCGCCAACGTTTAAAAAGGGGTGAAAAGCGATCGATGTAGCGAGCGTCGTGCg                                                                                                                                         >  2:815180/1‑238 (MQ=255)
                                                                                          ttgaTGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGGAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATTTAGCGAGCGGCGTGCGGCAGAAGATg                                                                                                                               >  2:488332/1‑237 (MQ=255)
                                                                                              tGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTTGAGGGGGGTAAGAGTCAGGCCCTGCAACCCGTGACAAAGGGGGGAAAAACGCAAGATGTAACGAGAGCCGGGCGGCGGAAAAGGccgcgg                                                                                                                         >  1:553764/1‑235 (MQ=255)
                                                                                               ttttatcttcaaactgtttcccAGCTTTGAACAGAAGCTGGATGGGGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGCCGCGCGACAATCCTCTATCTTCGCTGTCATTGTGGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGa                                                                                                                        <  1:38935/218‑1 (MQ=255)
                                                                                               gTCTATCTACAAACGTCGTCGCAGCTTTGAACATGAGCTGGATGGTGACAATATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGa                                                                                                                        <  2:27461/239‑1 (MQ=255)
                                                                                                    ttaaaaaaGGTCTGCGAAGCTTTGAACAGGAGGTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTTAGGTAAAAACGAGAAGTAATGTCGCGCGTCAATCCTCTAACTTCGCTGTCAGTGTTGGCGGTTCTGGTGTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGAGAGATGTGGCGGGGGTCGTGCGGCAGAAGATGCCTGCGAGCGTg                                                                                                                   <  1:669992/234‑1 (MQ=255)
                                                                                                                tGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGGTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGGGAGCGTCGTGCGGCAGAAGAGGCCTGGGGGCGGGAAAAAcccggat                                                                                                       >  2:435238/1‑234 (MQ=255)
                                                                                                                tGCGCAGCTTTGAACATGAGCTGGATGGTAACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTGCTGGGTTTTGCGGGGGGGAGTAGTACGTCTCCGCCACCCTTTTAAAACGGCTCAAAAGCGGAAGAGATATCGACGGTAGTGCGGCCGAAAATAACTGCGAgagcaaaaaaacgggag                                                                                                       >  1:412487/1‑223 (MQ=255)
                                                                                                                   tacgCTTTTAACTTTGGCTGGATGGTGACCAAATCAAGGAAGTAACCCATACCTGGGTTTAATAAGGGACTGGAATGTCGGCGGTCAATCGTCTATCTTGGTGGTCATTGTTGGCGGTTCTTGTTTTAGCGGCGCGTAGTAGTCAGGCGCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCTAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGcc                                                                                                    <  1:789302/236‑1 (MQ=255)
                                                                                                                   gcAGCTTTGAACATGAGCTTGATGGTTACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGc                                                                                                     <  2:428408/238‑1 (MQ=255)
                                                                                                                   gcAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGcc                                                                                                    <  2:461619/239‑1 (MQ=255)
                                                                                                                   gcAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGGCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGcc                                                                                                    <  2:262701/239‑1 (MQ=255)
                                                                                                                       cTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGc                                                                                                     <  2:141341/234‑1 (MQ=255)
                                                                                                                               aTGAGCTGGATGTGTACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATATATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGAt                                                                                        <  2:886798/239‑1 (MQ=255)
                                                                                                                                 gAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTTGAAAAGGGCGAAAAAGCGATAGATGTAGCGAGGGCCGGGCGGCAGAAGATGCCCGCGAGCGGGAAAGACCCCGATGGCCGGGGGAAAAAtg                                                                                       >  1:141341/1‑237 (MQ=255)
                                                                                                                                 gAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCCCCGCAACCGTTGGAAAAGAGCGAAAAAGCGATAGATGTAGGGAGCGGCCTGCGGCAGAAGAGGCCTGGGAGCGGGAAAGAACCCGAGGCCCTGGCGaaaaaaa                                                                                       >  1:79771/1‑234 (MQ=255)
                                                                                                                                 gAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGGTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGGGTCGTGCGGCAGAAGATGCCTGCGAGGGGGAAAGAACGCGAGGCCCGGGGGAAAGaac                                                                                       >  1:169457/1‑236 (MQ=255)
                                                                                                                                 gAGCTGGATCGGGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGGGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGACGGTGAAAAACCGCGAGACCGGGGGGGAAGaaa                                                                                       >  1:165941/1‑236 (MQ=255)
                                                                                                                                         tGGTGAAAAAATTCAGGAAGGAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTGGCGCGTCAATCCTCTATCTGCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTGAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGccacc                                                                               <  1:506929/238‑1 (MQ=255)
                                                                                                                                            tGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGCAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCCGAAGATGCCTGCGGGCGGGAAAAACCGCGGTGCCTGGGCGAAAGAGCTCGCccccac                                                                             >  2:634269/1‑237 (MQ=255)
                                                                                                                                                           ggTAACCCTTTCCTTTGGCTAAAAGGGGTGAGTTATGTCGCTCGGCAATGCGTTATCTGCGCTTGCATGGTTGGCTGTTCTGGTTTTAGCGGGGTGTAGTAGGCAGGCTCCGCAACCGGTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAgggg                                                             <  1:221560/237‑1 (MQ=255)
                                                                                                                                                                     aCCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCCGAAGATGCCTGCGAGCGTGAAAGACACCGAGGCCCGGGCGAAAAAACTCCCCAACCCCTTTTAAAGTCAGGGGCTGGCACCGa                                                  >  2:97450/1‑238 (MQ=255)
                                                                                                                                                                      ccTTGAGCTAAAAAGGAGAAGGAATGGCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACg                                               <  2:169457/239‑1 (MQ=255)
                                                                                                                                                                          gAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTTTGGTTTTATCGGGGGGGAAGAGTCAGGCTCCGCCACCGGTGGAAAAAGGCGAAAAAGCGAAAGATGTAACGTGCGGCGGGGTGCAGAAGAGGCCCGCGCGCGTGAAAGACCGCGCTGCCCGGGCGAGAGGTCTCCCCACCACCCTTGAAAATCAGGGGCCTGCCACCgcccagg                                             >  1:607676/1‑233 (MQ=255)
                                                                                                                                                                          gAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGGGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCCAGCGGAGTGCGGCAGGAGATGCCTGAGAGCGGGAAAAACCGCGATGCCTGGGGGTAAGATACCGCCCCCCCCTTTGAAAGGCCGGGGGGGGGACCggacgtgg                                            >  1:147462/1‑232 (MQ=255)
                                                                                                                                                                           aGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGGGTAGTAGTTAGGCTCCGCCACCGTTGAAAAAGGGCGGAAAAGCGAAAGATTTACGGGGCGCGGGTGGGCAGAAGAACGCCGCGGGGGGGAAAAACCGCGTGGCCTGGGGGGAAGAGCATGCCCCCCACATTGTAAGCCAGGGGCAGGcccgcgcgcggga                                           >  1:175539/1‑227 (MQ=255)
                                                                                                                                                                           aGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGGTTTCGCGGGGTGTTGTAGTAAGGCTCCGCAACCGATGCAAAAAGGCGGAAAACCGATCACGGTAGCGAGCCGCGGGCGGACTCAAATGCCAGCACACAGGACAAACCGCGCTGCCCGGGCGACAAAACTCGCCAACACCTTGGAAAAGCAGGGGCggggaaccacaccgga                                           >  1:190148/1‑226 (MQ=255)
                                                                                                                                                                             cTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCAGTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGGTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCTCCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGaa                                         <  2:172606/238‑1 (MQ=255)
                                                                                                                                                                                       gaattaatttcggcgtgcattCCTCTATCTCCGCTGCTATGTTTGGCGGTTCGGGTTTTAGCGGGGTTTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCt                               <  1:488332/218‑1 (MQ=255)
                                                                                                                                                                                       gaAGTAATGTCTCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCt                               <  2:450694/238‑1 (MQ=255)
                                                                                                                                                                                        aaGTAATGGCGCGCGTCAATCCTCTATCGTCGCTGTCATTGGTGGCGGTTTTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTc                              <  2:538783/238‑1 (MQ=255)
                                                                                                                                                                                                  gcgcgTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTTTGGTTTTAGCGGGGTGTAGTCGGCAGGGCCCCCAACCGTTGAAAAAAGGAGAAAAAGCCATAGAAGTTGCGCGCGGCGTGCGGCGAGAGAGGCCTGCGAGCGTGACAGACACGTATGCCTTGGCTCACAGTATCCCCACCACCTTTGCAAGTCAGGGCGAGGCCCCGACCCCTGAAAAACgcatcgccgggccgcggg                    >  1:840069/1‑222 (MQ=255)
                                                                                                                                                                                                  gcgcgTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCgg                   <  2:143005/239‑1 (MQ=255)
                                                                                                                                                                                                           ttctctttgttcggttgcATTGTTGGGGGTTTTGGTTTTAGCGGGGTGGAGTAGTCAGGCTCCGCAACGGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCgcagc          <  1:204078/222‑1 (MQ=255)
                                                                                                                                                                                                                  tCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGGTAGATGTAGCGAGCGGCGTGCGGGCGAAGATGCCGGCGAGCGTGAAAAACCACGATGCCTGGGCGGAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCGGGCGCCGGCGCTGGAAAACGACTGCCCGGTGCCGGGGGGGGGgcagacggcagc   >  1:654294/1‑231 (MQ=255)
                                                                                                                                                                                                                   cTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCa  <  2:175539/239‑1 (MQ=255)
                                                                                                                                                                                                                   cTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCa  <  2:213126/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                             
CTCGGCCTGATGACGCAGATTACCAACGTTTTTGGTCAGGTTCGCGGTGCTTTCCAGTACCTGATTAACTCATGGACCACACTGGTTGAGTTGATGTCTATCTACAAACGTCTGCGCAGCTTTGAACATGAGCTGGATGGTGACAAAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCA  >  NC_000913/397683‑398132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: