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.2% 127.0 / NA 47E54V (GAG→GTG) fimCtype 1 fimbriae periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base A (2/1);  major base T (16/25);  minor base G (2/1);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 5.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTG  >  NC_000913/4544228‑4544709
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
cGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTg                                                                                                                                                                                                                                          >  1:137050/1‑250 (MQ=255)
           aGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAatgat                                                                                                                                                                                                                              >  1:196753/1‑251 (MQ=255)
            ggCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAatgatg                                                                                                                                                                                                                             <  2:262476/251‑1 (MQ=255)
                 gggtctcTATAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATTGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGaaaa                                                                                                                                                                                                                         <  1:3664/247‑1 (MQ=255)
                    tctcTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAg                                                                                                                                                                                                                      <  1:428496/250‑1 (MQ=255)
                         ttAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGGGTGGCCTTAGGTGCGACTCGAGTAATTTATCCGGCAGGGCAAAACCAAGGGCAACTTTCCGTGCCAACAAATGATGAAAATAGTACCt                                                                                                                                                                                                                 >  2:213966/1‑250 (MQ=255)
                                  cAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGAGGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTc                                                                                                                                                                                                        <  1:53015/250‑1 (MQ=255)
                                      aaTTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTc                                                                                                                                                                                                        <  1:404216/246‑1 (MQ=255)
                                      aaTTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATc                                                                                                                                                                                                    >  2:48297/1‑250 (MQ=255)
                                      aaTTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGGGACAAATAATGCTGAAAATAGTACCTATTTAATTc                                                                                                                                                                                                        >  2:404216/1‑246 (MQ=255)
                                       aTTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTTTTTAATTCAATCAt                                                                                                                                                                                                  >  2:270169/1‑251 (MQ=255)
                                        ttGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGAAGGGCAAAAACAAGGGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTAAATcctg                                                                                                                                                                                                 >  2:288030/1‑248 (MQ=255)
                                                cagATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGGCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGaaa                                                                                                                                                                                         <  2:198560/251‑1 (MQ=255)
                                                    tAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGc                                                                                                                                                                                     >  1:490775/1‑251 (MQ=255)
                                                                  gaacaggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAgg                                                                                                                                                                       <  1:270169/251‑1 (MQ=255)
                                                                    acaggaCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGAt                                                                                                                                                                     <  1:149222/251‑1 (MQ=255)
                                                                         aCAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCg                                                                                                                                                                <  1:220623/251‑1 (MQ=255)
                                                                              gAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGtttt                                                                                                                                                            <  1:372796/250‑1 (MQ=255)
                                                                              gAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTa                                                                                                                                                           <  1:288030/251‑1 (MQ=255)
                                                                                    taaAAACGTAAATGTAAGGAAATCGCAGGACATCACATTCTGCTTGCTGGCAGGTCTCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAACATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTg                                                                                                                                                      >  2:461386/1‑250 (MQ=255)
                                                                                      aaaacttctagtgaatgAAATCGCAGGAAATAACAGTCTGCTGGATGGCAGGTATTCTGATGGTCATGGCAATGATGGTTGCCGGACGCGGTTAAGCGTGAGTGTCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGTCAAAAACAAGTGCAACTTGACGTGACAAATAATGATGAAAATAGTTCCTATTTAATTCATTCATGGGTGGAAAATGTCGATGGTGTAAAGGATGGTCGTTTTATCGTGACg                                                                                                                                                   <  1:302208/235‑1 (MQ=255)
                                                                                      aaaaCGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACg                                                                                                                                                   >  2:151342/1‑251 (MQ=255)
                                                                                                        aagcGCAAGAAATAACAGTCTGATTGCTGGCAGGTATCCCTATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTGAGGTGCGACTCGCGTAATTGATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGGCATTTTAGCGTGTCGCCACCTCTGTGTTCGATg                                                                                                                                 <  2:35335/248‑1 (MQ=255)
                                                                                                             cagtaaaTAACATTCCGCTTGCGGGCAGGTTTCCGGTGGTTCATGTCAAGGATGGTTGCCGGACGCGCTGAAGCGTGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGACAATATTCCCTATTTAATCTAAACATGTGTGGAAAATTCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCGGTTTGCGCTGAAggg                                                                                                                            <  2:359566/247‑1 (MQ=255)
                                                                                                              aGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGGCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGa                                                                                                                           <  1:99952/251‑1 (MQ=255)
                                                                                                                                 gcttgcaaggaTCCTGATGTTCATGTCAATGATGGTTGCCGTACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGTGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGAGGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATCCCGATGGTGTAAAGGATGTTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATTAAGGGAAAAAAAGAGAATACCTTAc                                                                                                        <  1:265439/241‑1 (MQ=255)
                                                                                                                                               tGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGc                                                                                           >  1:260217/1‑250 (MQ=255)
                                                                                                                                                       aTGGCAATGATGGTTGTCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCTTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATaa                                                                                  <  1:454637/251‑1 (MQ=255)
                                                                                                                                                          gCAATGATGGTTGCCGGCCGCGCTGAGGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAACCAAGTGCAAATTGCCGTAACAAATAATGCTGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGGTGGTGTAAAGGATGGTCGTTTTATCGTGACGACTCCTCTGGTTGCGATGAAGGGAAAAAACGAGAATACCTTACGGATTCTTGATTCAACACATAAcc                                                                                >  2:59055/1‑250 (MQ=255)
                                                                                                                                                                      gCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGGGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTTCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATACCCAATTGCCACAgg                                                                    <  2:77628/250‑1 (MQ=255)
                                                                                                                                                                       ccGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGCTGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGAc                                                                  >  1:65763/1‑251 (MQ=255)
                                                                                                                                                                        cGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAAGAACCAATTGCCACAGGAcc                                                                 <  1:72256/251‑1 (MQ=255)
                                                                                                                                                                          gACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAACACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTTTCGTGACGACTCCTCTGTTTGCGCTGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACTAATAACCACTTGCAACGGGACCgg                                                               >  2:233104/1‑251 (MQ=255)
                                                                                                                                                                                cTGACGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGGGACAAATAATGATGAAAATAGTACCTATTTAATTCAATAATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGCAGCCTCCTCTGTTTGCGATGAAGGGGAAAAAAGAGAATACCTTACGTATTCTTGTTGCAACAAATAACCAATTGCCACAGGACCGGGAAAGt                                                         >  2:368478/1‑251 (MQ=255)
                                                                                                                                                                                 tGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTaa                                                                                                                                                                          <  1:28898/137‑1 (MQ=255)
                                                                                                                                                                                 tGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTaa                                                                                                                                                                          >  2:28898/1‑137 (MQ=255)
                                                                                                                                                                                                   aGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGtt                                       <  1:398951/250‑1 (MQ=255)
                                                                                                                                                                                                                     ttATCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAGCCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATg                     <  1:461386/250‑1 (MQ=255)
                                                                                                                                                                                                                        tCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcac                                                                       >  1:163896/1‑197 (MQ=255)
                                                                                                                                                                                                                        tCCGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCGATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCcac                                                                       <  2:163896/197‑1 (MQ=255)
                                                                                                                                                                                                                        gccGGCAGGGCAAAAACAAATGCCACTTGCAGTGCCAAATAATGATGAAGATAGTACATCTTTAATTCAATCATGGGTGGAACATGCCGATGGTGTAAAGGCTGGTCGTTTTATCCTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGCATACATTAAGTATTATTGCTGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGTTGAAAGTTaaacagattcagtaaaaggga                  >  2:469027/2‑232 (MQ=255)
                                                                                                                                                                                                                         ccGGCAGGGCAAAAACAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATa                 <  1:48297/250‑1 (MQ=255)
                                                                                                                                                                                                                          cGGCAGGGCCAAAACAAGTGCCACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGAAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAACAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACCGGACCGGGAAAGTTTCTTCTGGATGAACGTAAAAGCGCTTCCGTCAATGGATaa                >  2:26423/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGgtg                                                                                                                                                                             >  2:112065/1‑86 (MQ=255)
                                                                                                                                                                                                                                 gCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGgtg                                                                                                                                                                             <  1:112065/86‑1 (MQ=255)
                                                                                                                                                                                                                                      aaCAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTCCCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGAGGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATTGATAAATCAAAATTGACt   <  1:311470/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aCAAGTGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGGCAATGGATAAATCAAAATGGACTg  <  2:445360/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
CGCCAATGCCCAGGCCTGGTTCTCTTTAACCTATCAGTAATTGTTCAGCAGATAATGTGATAACAGGAACAGGACAGTGAGTAATAAAAACGTCAATGTAAGGAAATCGCAGGAAATAACATTCTGCTTGCTGGCAGGTATCCTGATGTTCATGGCAATGATGGTTGCCGGACGCGCTGAAGCGGGAGTGGCCTTAGGTGCGACTCGCGTAATTTATCCGGCAGGGCAAAAACAAGAGCAACTTGCCGTGACAAATAATGATGAAAATAGTACCTATTTAATTCAATCATGGGTGGAAAATGCCGATGGTGTAAAGGATGGTCGTTTTATCGTGACGCCTCCTCTGTTTGCGATGAAGGGAAAAAAAGAGAATACCTTACGTATTCTTGATGCAACAAATAACCAATTGCCACAGGACCGGGAAAGTTTATTCTGGATGAACGTTAAAGCGATTCCGTCAATGGATAAATCAAAATTGACTG  >  NC_000913/4544228‑4544709

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