Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,001,174 C→T 5.9% S653S (AGC→AGT elfC → putative outer membrane fimbrial subunit export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,001,1740CT 5.9% 210.7 / 2.5 67S653S (AGC→AGTelfCputative outer membrane fimbrial subunit export usher protein
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base T (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

CTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGG  >  NC_000913/1001030‑1001321
                                                                                                                                                |                                                                                                                                                   
cTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTGCTGGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAgcgc                                                                                                                                                   >  5:178106/1‑147 (MQ=255)
   gATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCa                                                                                                                                              >  4:93187/1‑149 (MQ=255)
     tGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATg                                                                                                                                            >  5:78392/1‑149 (MQ=255)
         ctctAACACCAGTATCGATAACGAAGGGCACACTACACTAAACCTGTGTTTAACGGATACATTACTCGATGAGGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGAtt                                                                                                                                        <  5:102336/149‑1 (MQ=255)
          tctAACACCAGTATCGATAACGAAGGGCACACTACCCTAACCCGGGGTGTCACGGAGACATTACTCGCGGACGGTGACCTGAGCGACAGCGTGGCACAGGGATATAACAGCGCGGGGAACACGGCGAATGGGAGCGCCAGCAGGGagga                                                                                                                                       >  6:231082/1‑146 (MQ=255)
          tctAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTa                                                                                                                                       >  6:27083/1‑149 (MQ=255)
           ctAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTa                                                                                                                                       >  7:198478/1‑148 (MQ=255)
            tAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTa                                                                                                                                       <  5:27083/147‑1 (MQ=255)
              acacCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACa                                                                                                                                     >  3:311626/1‑147 (MQ=255)
               cacCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAg                                                                                                                                  <  8:233102/149‑1 (MQ=255)
                acCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCACCAGGGATATAACAGCGAGGGGAATACGGCGAATGGTAGCGCCAGCAGGGATTACAAAgg                                                                                                                                 >  1:246746/1‑149 (MQ=255)
                 ccAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAggg                                                                                                                                <  1:181842/149‑1 (MQ=255)
                    gTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCTATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACCGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGAGTACAACGGGGCg                                                                                                                             >  3:17389/1‑149 (MQ=255)
                     tATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGt                                                                                                                            >  6:59603/1‑149 (MQ=255)
                       tCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCCGCATGGATTACAAAGGGGCGttt                                                                                                                          >  1:224042/1‑149 (MQ=255)
                         gATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGc                                                                                                                        <  7:35831/149‑1 (MQ=255)
                              cGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTg                                                                                                                         <  2:7431/143‑1 (MQ=255)
                              cGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTg                                                                                                                         >  1:7431/1‑143 (MQ=255)
                               gAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAAGGGCGTTTGCGGATGc                                                                                                                  >  1:195668/1‑149 (MQ=255)
                                 aGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGccc                                                                                                                >  6:23033/1‑149 (MQ=255)
                                    gCACACTACACAAAACCTGGGTTTAACGGAGACATTACTAGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAACGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAg                                                                                                             <  8:282753/149‑1 (MQ=255)
                                     cacaCTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGt                                                                                                            >  3:9359/1‑149 (MQ=255)
                                          tACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGTGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGc                                                                                                        >  8:271968/1‑148 (MQ=255)
                                            cacaAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTa                                                                                                      >  6:319590/1‑148 (MQ=255)
                                               aaaaCCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGTGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAAc                                                                                                  <  7:271968/149‑1 (MQ=255)
                                               aaaaCCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGACCGAGTGGGCTACAAc                                                                                                  >  2:251126/1‑149 (MQ=255)
                                                    cTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAg                                                                                             >  6:173989/1‑149 (MQ=255)
                                                       ggTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGc                                                                                            <  3:265795/147‑1 (MQ=255)
                                                            aaCGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGTGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACg                                                                                     >  8:80345/1‑149 (MQ=255)
                                                            aaCGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAAc                                                                                      <  3:93187/148‑1 (MQ=255)
                                                            aaCGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACg                                                                                     >  4:282321/1‑149 (MQ=255)
                                                             aCGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACg                                                                                     <  6:296092/148‑1 (MQ=255)
                                                             aCGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACg                                                                                     <  1:135716/148‑1 (MQ=255)
                                                                gagaCATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTTATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAg                                                                                 >  4:225786/1‑149 (MQ=255)
                                                                     aTTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTcaac                                                                            <  7:104969/149‑1 (MQ=255)
                                                                       tACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTcaacaa                                                                          >  4:63449/1‑149 (MQ=255)
                                                                        aCTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTg                                                                                                           >  8:2442/1‑115 (MQ=255)
                                                                          tCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTcaacaacaa                                                                       >  2:6843/1‑149 (MQ=255)
                                                                            gATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCaacaacaac                                                                      <  5:243524/148‑1 (MQ=255)
                                                                            gATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCaacaacaac                                                                      <  5:59603/148‑1 (MQ=255)
                                                                             aTGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTg                                                                    <  8:90882/149‑1 (MQ=255)
                                                                                 cGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAAc                                                                 >  6:5199/1‑148 (MQ=255)
                                                                                 cGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAAc                                                                 >  8:237460/1‑148 (MQ=255)
                                                                                  ggTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTa                                                               <  8:171831/149‑1 (MQ=255)
                                                                                          gAGCTACCGCGTGCAACAGGGATATAACCGCGAGGTGACAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCTGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCttt                                                       <  4:17389/149‑1 (MQ=255)
                                                                                             cTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAg                                                    <  4:311626/149‑1 (MQ=255)
                                                                                                   cGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGtt                                              <  6:158309/149‑1 (MQ=255)
                                                                                                       cAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGt                                          <  1:284144/149‑1 (MQ=255)
                                                                                                       cAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGt                                          <  1:251126/149‑1 (MQ=255)
                                                                                                           aGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGcc                                      <  4:9359/149‑1 (MQ=255)
                                                                                                             ggATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGc                                                                                            <  6:153053/93‑1 (MQ=255)
                                                                                                             ggATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGc                                                                                            >  5:153053/1‑93 (MQ=255)
                                                                                                               atatAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCAtt                                  >  5:153860/1‑149 (MQ=255)
                                                                                                               atatAACAGCGAGGGGAAAACGGCAAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCAtt                                  <  5:248417/149‑1 (MQ=255)
                                                                                                                tataACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTc                                 <  3:63449/149‑1 (MQ=255)
                                                                                                                  taACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGc                                       <  1:137342/141‑1 (MQ=255)
                                                                                                                  taACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGc                                       >  2:137342/1‑141 (MQ=255)
                                                                                                                    aCAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGTCCATTCACAg                             <  3:282321/149‑1 (MQ=255)
                                                                                                                     cAGCGAGGGGAAAACGGCTAATGGTAGTGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAgg                            <  7:80345/149‑1 (MQ=255)
                                                                                                                      aGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACACAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGta                                         <  5:226666/135‑2 (MQ=255)
                                                                                                                      aGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGtt                                         >  6:226666/1‑135 (MQ=255)
                                                                                                                      aGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTCAACTACGCTCTTTCAGACAGTTTAGTTGCCCATTCACAggg                           >  8:202548/1‑149 (MQ=255)
                                                                                                                                  aCGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTggg               <  8:198478/149‑1 (MQ=255)
                                                                                                                                  aCGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTggg               <  1:6843/149‑1 (MQ=255)
                                                                                                                                  aCGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTggg               <  8:43972/149‑1 (MQ=255)
                                                                                                                                       tAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAAt          >  4:264253/1‑149 (MQ=255)
                                                                                                                                            gTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTg     >  1:103909/1‑149 (MQ=255)
                                                                                                                                              aGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTggg   >  4:7556/1‑149 (MQ=255)
                                                                                                                                               gcgcCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGACCTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTgggg  >  6:288281/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGCGAGGGGAAAACGGCTAATGGTAGCGCCAGCATGGATTACAAAGGGGCGTTTGCGGATGCCCGAGTGGGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGG  >  NC_000913/1001030‑1001321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: