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,4640AT93.9% 135.1 / ‑0.8 49E54V (GAG→GTG) fimCtype 1 fimbriae periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (23/23);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGA  >  NC_000913/4544221‑4544710
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
gcggcATCGTCAATGCCCAGGCCTGGTTCTCGTTAACCTATCAGTAATTGGTCAGCAGCTAATGTTATAACAGGAACAGTACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACAGTCTGCTTGCTTGCAGGTATCCTGATGTGCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGAGCCTCGTAATTTATCCGGCAGGGCAAAAACAAGTGTAACtt                                                                                                                                                                                                                                                  <  2:427482/251‑1 (MQ=255)
   acgcgctaccgttctaGGCCTGGTTCTCTTTAACCTATCAGTAAGTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTAGCCTGATGTTCATGGCAATGATGGTTGCCGGACGCTCTGAAGCCGGAGTGGCCGTACGTGCGACTCGCGGAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGcc                                                                                                                                                                                                                                               <  2:101233/236‑1 (MQ=255)
     aTCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGt                                                                                                                                                                                                                                             >  1:446349/1‑251 (MQ=255)
          cAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCATGGCAAAAACAAGAGCAACTTGCCGTGACa                                                                                                                                                                                                                                         >  2:60551/1‑250 (MQ=255)
                 cAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAAGACAGTGAGTAATAAAAACGTCAAGGGAATGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAATGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGTAAAAACAAGTGCAACTTGCCGTGACGAATAatga                                                                                                                                                                                                                                 <  2:653836/251‑1 (MQ=255)
                  aGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGcc                                                                                                                                                                                                                                               >  1:101233/1‑236 (MQ=255)
                      cTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGaa                                                                                                                                                                                                                             >  1:419754/1‑250 (MQ=255)
                                 tAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTAt                                                                                                                                                                                                                 >  1:506598/1‑251 (MQ=255)
                                       aTCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAAt                                                                                                                                                                                                            <  1:340611/250‑1 (MQ=255)
                                            tAATTGTTCAGTAGATAATGTGATAACAGGAATATTACAGTGATTAATATATATGTCTATTTCAGGTAATCGTATTAAATCACTTTCTGCGTGTTTGCTGGTATCTTGATGTTCATGGCAATGATTGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCTTTTATGTTAAAAACATTTTCATCTTGCCGTGACTAATAATGATGAAAATAGTACCTATTTAATTCAAt                                                                                                                                                                                                       <  2:104847/250‑1 (MQ=255)
                                                tGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATg                                                                                                                                                                                                   <  1:290096/250‑1 (MQ=255)
                                                 gTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATTATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATggg                                                                                                                                                                                                 <  1:652879/251‑1 (MQ=255)
                                                            aaTGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGcc                                                                                                                                                                                      >  1:271691/1‑251 (MQ=255)
                                                                 gATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGATTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGTCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATg                                                                                                                                                                                  <  1:420788/250‑1 (MQ=255)
                                                                 gATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACCCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATgg                                                                                                                                                                                 >  1:630700/1‑251 (MQ=255)
                                                                         gaacaggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAgg                                                                                                                                                                         >  1:616218/1‑251 (MQ=255)
                                                                             aggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATgg                                                                                                                                                                     <  1:55607/251‑1 (MQ=255)
                                                                                   gTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACtt                                                                                                                                                                                                                                                  >  2:186782/1‑168 (MQ=255)
                                                                                   gTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACtt                                                                                                                                                                                                                                                  <  1:186782/168‑1 (MQ=255)
                                                                                                          aaGGAAATCGCAGGAAATAACATTCTGCTTTCTGGCAGGTATCCTGATGTTCATGGCAATTATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGCCGCCTCCTCTGttt                                                                                                                                         <  2:10021/250‑1 (MQ=255)
                                                                                                               aaTCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGCTGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCCGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATATTACCTATTTAATTCa                                                                                                                                                                                                         >  1:402261/1‑181 (MQ=255)
                                                                                                               aaTCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATTGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAACTAGGACCTATTTAATTCAATCATGGGTGGAAAATGACGATAGTGTAAAGGATGGTCTTTTTATCGTTACGCTCCCTATGTTTGCGATg                                                                                                                                   >  2:296/1‑251 (MQ=255)
                                                                                                               aaTCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATg                                                                                                                                   >  1:121825/1‑251 (MQ=255)
                                                                                                               aaTCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAAT‑‑ATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCCTGGGTGGAAAATGCCAATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGa                                                                                                                                  >  2:336322/1‑250 (MQ=255)
                                                                                                               aaTCGCAGGAAATAACATGATGCTTGCTGGCAGGGAGCCTGATGTGCATGGGAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCGTAGGTGCGACTTGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAAGAGTACCTATTTAAGTCa                                                                                                                                                                                                         <  2:402261/181‑1 (MQ=255)
                                                                                                                  cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAg                                                                                                                                <  1:566338/251‑1 (MQ=255)
                                                                                                                  cGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAg                                                                                                                                <  1:163136/251‑1 (MQ=255)
                                                                                                                      ggAAATAACATTCTGCTTACTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGAGGAAAATGCCGATGGTGTAAAGTATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGa                                                                                                                             >  1:188247/1‑250 (MQ=255)
                                                                                                                        aaaTAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGTTGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGActcc                                                                                                                                                   >  2:467138/1‑223 (MQ=255)
                                                                                                                        aaaTAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGcc                                                                                                                                                   <  1:467138/226‑1 (MQ=255)
                                                                                                                        aaaTAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaaaa                                                                                                                          <  1:527394/251‑1 (MQ=255)
                                                                                                                             aCATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaaaaaa                                                                                                                        >  1:133642/1‑248 (MQ=255)
                                                                                                                             aCATTCTGCTTGCTGGCAGGTATCCGGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGaaaaaa                                                                                                                        <  2:133642/248‑1 (MQ=255)
                                                                                                                                                    ccTGATGTTCATGTCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATg                                                                                              <  1:373342/251‑1 (MQ=255)
                                                                                                                                                              aTGGCAATGATGTTTGCCGTACTCTCTGAAGCGGGAGTAGACTTAGGTGCGACTCGCGTAATTTATACGGCAGGGAAAAAACAAGTGCAACATGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAAGACCTTACGTATTCTTGATGCAataaataa                                                                                    <  2:544694/251‑1 (MQ=255)
                                                                                                                                                                ggCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAAcc                                                                                  <  2:446349/251‑1 (MQ=255)
                                                                                                                                                                  cAATGATGGTTGCCGGACGCGCGTAAGAGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTAACCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTGCCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGGATTCTTTATGCAACAAATAACCaa                                                                                <  2:470697/251‑1 (MQ=255)
                                                                                                                                                                         ggTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCca                                                                          <  1:226897/250‑1 (MQ=255)
                                                                                                                                                                                             cGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTc                                                     >  1:156374/1‑251 (MQ=255)
                                                                                                                                                                                                  gTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGAt                                               >  1:219437/1‑251 (MQ=255)
                                                                                                                                                                                                      ccTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAAc                                           >  1:601156/1‑251 (MQ=255)
                                                                                                                                                                                                               cgccTCGCGTAATTTATCCGGAAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGGAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGGTTGCGATGAAGGGGAAAAAAGAGAATACCTTACGTATTCTTGATGCATCAACTAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGAGAGAACGTTAAAGCg                                  >  2:21224/4‑251 (MQ=255)
                                                                                                                                                                                                                      cgTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGt                           >  1:114871/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaTTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGCCAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACGGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCa                         >  2:545382/1‑250 (MQ=255)
                                                                                                                                                                                                                          aTTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAAt                       >  1:383845/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATc              >  1:492215/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ggCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCtt                                                                                                  >  2:600000/1‑164 (MQ=255)
                                                                                                                                                                                                                                       ggCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCGGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCtt                                                                                                  <  1:600000/164‑1 (MQ=255)
                                                                                                                                                                                                                                               cAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCAATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTga  <  1:336322/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
GCGGCATCGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTGA  >  NC_000913/4544221‑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: