Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,544,464 A→T E54V (GAG→GTG)  fimC → type 1 fimbriae periplasmic chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,544,4640AT98.3% 176.4 / ‑6.7 60E54V (GAG→GTG) fimCtype 1 fimbriae periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (24/34);  minor base C (0/1);  total (24/36)
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.99e-01

GGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGA  >  NC_000913/4544220‑4544710
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
ggcggcATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCaa                                                                                                                                                                                                                                                     <  2:301195/249‑1 (MQ=255)
 gcggcATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACtt                                                                                                                                                                                                                                                  <  1:439036/251‑1 (MQ=255)
   ggcATCGCCAATGCCCAGCCCTGGTTATCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATCACAGGAACAGGACAGTGCGTAATAAACACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGGGCGACTCGCGTAATTTATCCGGCAGGGCACAAACAAGTGCAACTTGc                                                                                                                                                                                                                                                >  2:587376/1‑251 (MQ=255)
         gCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGAc                                                                                                                                                                                                                                          >  2:67420/1‑251 (MQ=255)
               gCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAaataa                                                                                                                                                                                                                                    <  1:136073/251‑1 (MQ=255)
                cccAGGCCTGGTTCTCTGTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAaataat                                                                                                                                                                                                                                   <  2:187967/251‑1 (MQ=255)
                 ccAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGGAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAatg                                                                                                                                                                                                                                  <  2:486347/251‑1 (MQ=255)
                   aGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTAGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAatgat                                                                                                                                                                                                                                <  2:574046/251‑1 (MQ=255)
                     gCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAAtgatga                                                                                                                                                                                                                              <  1:361569/251‑1 (MQ=255)
                        tGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGaaa                                                                                                                                                                                                                            >  2:283918/1‑250 (MQ=255)
                         ggTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGaaaa                                                                                                                                                                                                                           <  1:451650/250‑1 (MQ=255)
                          gTTCTCTTTAACCTATCAGTAATTGTTCATCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATa                                                                                                                                                                                                                         <  1:214672/251‑1 (MQ=255)
                                    aCCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTAttt                                                                                                                                                                                                               >  1:7939/1‑251 (MQ=255)
                                        aTCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACGAATAATGATGAAAATAGTACCTATTTaa                                                                                                                                                                                                             <  1:67420/249‑1 (MQ=255)
                                        aTCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑tt                                                                                                                                                                                                          <  1:320760/251‑1 (MQ=255)
                                          cAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGACGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGAACCTATTTAA‑TTc                                                                                                                                                                                                         >  1:421187/1‑250 (MQ=255)
                                          aaGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAAAAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTc                                                                                                                                                                                                         >  2:348805/2‑250 (MQ=255)
                                                 tGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACAGTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATg                                                                                                                                                                                                  <  1:587376/250‑1 (MQ=255)
                                                  gTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTGCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAATTTGCCGTGACAAATAATGATGAAAAGAGTACCGATTTAA‑TTCAATCATggg                                                                                                                                                                                                <  2:27778/251‑1 (MQ=255)
                                                          gATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATGCTGCTTGCTGGCAGGTAGCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCATAATTTATCCGGCAGGGCAAAAACAAGTGCAACGTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAAt                                                                                                                                                                                        <  2:392484/251‑1 (MQ=255)
                                                            aaaTGTGATAACAGGAACCGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCCGGAAATAACGTTCTGCTTGCTAGCAGGTATCCTGATGTTCATGGCAATGCTGGTTGCCGGACGCGCTGAAGCGGGAGTGTCCTTAGGTGCGACTCGCGTAATTGATCCGGACGGGCAAAAACACGTTCAACTTGCCGTGACAAATAATGATGACAATAGTACCTATTTAA‑TTCAATCATGGGTTGAAAATGc                                                                                                                                                                                      >  2:445552/2‑251 (MQ=255)
                                                                gtgATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGAt                                                                                                                                                                                  >  1:482091/1‑251 (MQ=255)
                                                                 tgATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAACCAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGGGGAAAATGCCGATg                                                                                                                                                                                 >  2:497077/1‑251 (MQ=255)
                                                                           aacaggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGAGGAAAATGCCGATGGTGTAAAggg                                                                                                                                                                       >  1:358208/1‑250 (MQ=255)
                                                                              aggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATgg                                                                                                                                                                    <  2:288733/251‑1 (MQ=255)
                                                                                     tGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGtttt                                                                                                                                                             >  2:285868/1‑251 (MQ=255)
                                                                                          aataaAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGt                                                                                                                                                        >  1:282132/1‑251 (MQ=255)
                                                                                                aaCGTCACTGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGa                                                                                                                                                                                                                                           >  2:244357/1‑163 (MQ=255)
                                                                                                aaCGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGa                                                                                                                                                                                                                                           <  1:244357/163‑1 (MQ=255)
                                                                                                  cGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGACTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTAGTTTTATCGTGACGcctc                                                                                                                                                >  1:355400/1‑251 (MQ=255)
                                                                                                           aaGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTg                                                                                                                                       <  2:192162/251‑1 (MQ=255)
                                                                                                              gAAATAGCAGGAAATAACATCCTGCTTGATGGCAGGAATCATGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGACTTAGGTGCGACTCGCGTACTTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATACTGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATAGTGACGACTCCTCTGTTTGCg                                                                                                                                     >  2:623638/1‑250 (MQ=255)
                                                                                                               aaaTCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGAt                                                                                                                                   >  2:176896/1‑251 (MQ=255)
                                                                                                                   cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGaa                                                                                                                                <  1:111320/250‑1 (MQ=255)
                                                                                                                   cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAAGCATGGGTGGAAAATGCCGATGGTGGAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGaa                                                                                                                                <  1:281018/250‑1 (MQ=255)
                                                                                                                   cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCGGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTACCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAAGCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCGCCTCTGTTTGCGATGAAg                                                                                                                               <  2:61229/251‑1 (MQ=255)
                                                                                                                       ggAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGAGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGa                                                                                                                            >  2:198837/1‑250 (MQ=255)
                                                                                                                       gaaaaTAACATGCTGCTGGCCGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaa                                                                                                                           <  2:3224/249‑1 (MQ=255)
                                                                                                                              aCATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaaaaaa                                                                                                                       >  1:48060/1‑248 (MQ=255)
                                                                                                                              aCATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAaga                                                                                                                    <  2:427243/251‑1 (MQ=255)
                                                                                                                                   cTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAt                                                                                                                >  2:586487/1‑250 (MQ=255)
                                                                                                                                     gcttgctGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGCGCAACTTTCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTAGCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAt                                                                                                                <  1:45080/248‑1 (MQ=255)
                                                                                                                                            ggCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACg                                                                                                        <  1:283918/249‑1 (MQ=255)
                                                                                                                                                     ccagATGTTCATGGCAAAGATGGTTGCCGTACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAGCAAGTGCAACTTGCCGAGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATg                                                                                             >  2:223630/4‑251 (MQ=255)
                                                                                                                                                     ccTGATGTTCAGGGTAGTGATGGTTGTCTGAGGCGCTGAAGCGGTAGTGGCTATAGGGGCTAATCGTGGAATTTAGCCGGCAGGGCAAAAACAAAGGCAACTTGCCGTGACAAATAATTATGAAAATAGTACCTATGTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAATGATGGTCGTTTCATCGTGTCGCCTCCTCTTTGTGCGATGTAGGGAAATAAAgagagtcctctgcgtatgcttgatg                                                                                             <  2:55142/251‑22 (MQ=255)
                                                                                                                                                         aTGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGcaac                                                                                         >  2:542934/1‑251 (MQ=255)
                                                                                                                                                           gTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGcaacaa                                                                                       <  1:113485/251‑1 (MQ=255)
                                                                                                                                                           gTTCATGGCAATGATGGGTGCCGGACGCGCTGAAGCGGGAGTGGCATTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGaaaca                                                                                        >  1:173415/1‑250 (MQ=255)
                                                                                                                                                                 ggCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAAc                                                                                  <  1:87791/250‑1 (MQ=255)
                                                                                                                                                                                 ggACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCctct                                                                                                                                            <  1:302487/176‑1 (MQ=255)
                                                                                                                                                                                 ggACGAGCTGAAGCGGGCGTGGCCTTAGGTGCGCCTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGCCAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCctct                                                                                                                                            >  2:302487/1‑176 (MQ=255)
                                                                                                                                                                                  gACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCgg                                                                <  2:434349/251‑1 (MQ=255)
                                                                                                                                                                                    cgcgCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TGCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                                              <  1:604945/251‑1 (MQ=255)
                                                                                                                                                                                                      gCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAGACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGaa                                            <  1:319712/251‑1 (MQ=255)
                                                                                                                                                                                                       ccTTAGGTGCGACTCTCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATAAAc                                           <  2:466183/251‑1 (MQ=255)
                                                                                                                                                                                                           aGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                                        <  1:457526/250‑1 (MQ=255)
                                                                                                                                                                                                              tGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGcc                                                                                                                                                  <  2:529787/141‑1 (MQ=255)
                                                                                                                                                                                                              tGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGcc                                                                                                                                                  >  1:529787/1‑141 (MQ=255)
                                                                                                                                                                                                                                          aaaaaaCAAGTGCAACTTGCCGTGACAAAGAATGATGAAAATAGTACCGATTTAACTTCA‑TTACGGGTGGAAAATGTCGATGGTGTAAAGGACGGTCGCTTTACCGTGCCGCATCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAATAAATAACCAATTGCCACAGGATCTGGACAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAAtt        <  1:445552/250‑1 (MQ=255)
                                                                                                                                                                                                                                                cAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTga  <  1:2712/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
GGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAA‑TTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGA  >  NC_000913/4544220‑4544710

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