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

GGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGA  >  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)
                                        aTCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAAtt                                                                                                                                                                                                          <  1:320760/251‑1 (MQ=255)
                                          cAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGACGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGAACCTATTTAATTc                                                                                                                                                                                                         >  1:421187/1‑250 (MQ=255)
                                          aaGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAAAAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTc                                                                                                                                                                                                         >  2:348805/2‑250 (MQ=255)
                                                 tGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACAGTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATg                                                                                                                                                                                                  <  1:587376/250‑1 (MQ=255)
                                                  gTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTGCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAATTTGCCGTGACAAATAATGATGAAAAGAGTACCGATTTAATTCAATCATggg                                                                                                                                                                                                <  2:27778/251‑1 (MQ=255)
                                                          gATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATGCTGCTTGCTGGCAGGTAGCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCATAATTTATCCGGCAGGGCAAAAACAAGTGCAACGTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAAt                                                                                                                                                                                        <  2:392484/251‑1 (MQ=255)
                                                            aaaTGTGATAACAGGAACCGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCCGGAAATAACGTTCTGCTTGCTAGCAGGTATCCTGATGTTCATGGCAATGCTGGTTGCCGGACGCGCTGAAGCGGGAGTGTCCTTAGGTGCGACTCGCGTAATTGATCCGGACGGGCAAAAACACGTTCAACTTGCCGTGACAAATAATGATGACAATAGTACCTATTTAATTCAATCATGGGTTGAAAATGc                                                                                                                                                                                      >  2:445552/2‑251 (MQ=255)
                                                                gtgATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGAt                                                                                                                                                                                  >  1:482091/1‑251 (MQ=255)
                                                                 tgATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAACCAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGGGGAAAATGCCGATg                                                                                                                                                                                 >  2:497077/1‑251 (MQ=255)
                                                                           aacaggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGAGGAAAATGCCGATGGTGTAAAggg                                                                                                                                                                       >  1:358208/1‑250 (MQ=255)
                                                                              aggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATgg                                                                                                                                                                    <  2:288733/251‑1 (MQ=255)
                                                                                     tGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGtttt                                                                                                                                                             >  2:285868/1‑251 (MQ=255)
                                                                                          aataaAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGt                                                                                                                                                        >  1:282132/1‑251 (MQ=255)
                                                                                                aaCGTCACTGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGa                                                                                                                                                                                                                                           >  2:244357/1‑163 (MQ=255)
                                                                                                aaCGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGa                                                                                                                                                                                                                                           <  1:244357/163‑1 (MQ=255)
                                                                                                  cGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGACTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTAGTTTTATCGTGACGcctc                                                                                                                                                >  1:355400/1‑251 (MQ=255)
                                                                                                           aaGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTg                                                                                                                                       <  2:192162/251‑1 (MQ=255)
                                                                                                              gAAATAGCAGGAAATAACATCCTGCTTGATGGCAGGAATCATGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGACTTAGGTGCGACTCGCGTACTTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATACTGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATAGTGACGACTCCTCTGTTTGCg                                                                                                                                     >  2:623638/1‑250 (MQ=255)
                                                                                                               aaaTCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGAt                                                                                                                                   >  2:176896/1‑251 (MQ=255)
                                                                                                                   cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGaa                                                                                                                                <  1:111320/250‑1 (MQ=255)
                                                                                                                   cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAAGCATGGGTGGAAAATGCCGATGGTGGAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGaa                                                                                                                                <  1:281018/250‑1 (MQ=255)
                                                                                                                   cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCGGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTACCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAAGCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCGCCTCTGTTTGCGATGAAg                                                                                                                               <  2:61229/251‑1 (MQ=255)
                                                                                                                       ggAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGAGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGa                                                                                                                            >  2:198837/1‑250 (MQ=255)
                                                                                                                       gaaaaTAACATGCTGCTGGCCGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaa                                                                                                                           <  2:3224/249‑1 (MQ=255)
                                                                                                                              aCATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaaaaaa                                                                                                                       >  1:48060/1‑248 (MQ=255)
                                                                                                                              aCATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAaga                                                                                                                    <  2:427243/251‑1 (MQ=255)
                                                                                                                                   cTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAt                                                                                                                >  2:586487/1‑250 (MQ=255)
                                                                                                                                     gcttgctGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGCGCAACTTTCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTAGCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAt                                                                                                                <  1:45080/248‑1 (MQ=255)
                                                                                                                                            ggCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACg                                                                                                        <  1:283918/249‑1 (MQ=255)
                                                                                                                                                     ccagATGTTCATGGCAAAGATGGTTGCCGTACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAGCAAGTGCAACTTGCCGAGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATg                                                                                             >  2:223630/4‑251 (MQ=255)
                                                                                                                                                     ccTGATGTTCAGGGTAGTGATGGTTGTCTGAGGCGCTGAAGCGGTAGTGGCTATAGGGGCTAATCGTGGAATTTAGCCGGCAGGGCAAAAACAAAGGCAACTTGCCGTGACAAATAATTATGAAAATAGTACCTATGTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAATGATGGTCGTTTCATCGTGTCGCCTCCTCTTTGTGCGATGTAGGGAAATAAAgagagtcctctgcgtatgcttgatg                                                                                             <  2:55142/251‑22 (MQ=255)
                                                                                                                                                         aTGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGcaac                                                                                         >  2:542934/1‑251 (MQ=255)
                                                                                                                                                           gTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGcaacaa                                                                                       <  1:113485/251‑1 (MQ=255)
                                                                                                                                                           gTTCATGGCAATGATGGGTGCCGGACGCGCTGAAGCGGGAGTGGCATTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGaaaca                                                                                        >  1:173415/1‑250 (MQ=255)
                                                                                                                                                                 ggCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAAc                                                                                  <  1:87791/250‑1 (MQ=255)
                                                                                                                                                                                 ggACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCctct                                                                                                                                            <  1:302487/176‑1 (MQ=255)
                                                                                                                                                                                 ggACGAGCTGAAGCGGGCGTGGCCTTAGGTGCGCCTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGCCAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCctct                                                                                                                                            >  2:302487/1‑176 (MQ=255)
                                                                                                                                                                                  gACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCgg                                                                <  2:434349/251‑1 (MQ=255)
                                                                                                                                                                                    cgcgCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATGCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGa                                                              <  1:604945/251‑1 (MQ=255)
                                                                                                                                                                                                      gCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAGACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGaa                                            <  1:319712/251‑1 (MQ=255)
                                                                                                                                                                                                       ccTTAGGTGCGACTCTCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATAAAc                                           <  2:466183/251‑1 (MQ=255)
                                                                                                                                                                                                           aGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                                        <  1:457526/250‑1 (MQ=255)
                                                                                                                                                                                                              tGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGcc                                                                                                                                                  <  2:529787/141‑1 (MQ=255)
                                                                                                                                                                                                              tGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGcc                                                                                                                                                  >  1:529787/1‑141 (MQ=255)
                                                                                                                                                                                                                                          aaaaaaCAAGTGCAACTTGCCGTGACAAAGAATGATGAAAATAGTACCGATTTAACTTCATTACGGGTGGAAAATGTCGATGGTGTAAAGGACGGTCGCTTTACCGTGCCGCATCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAATAAATAACCAATTGCCACAGGATCTGGACAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAAtt        <  1:445552/250‑1 (MQ=255)
                                                                                                                                                                                                                                                cAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTga  <  1:2712/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
GGCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGA  >  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: