Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,003,134 T→G 23.0% C82W (TGT→TGG ycbU → putative fimbriae protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,003,1340TG23.0% 113.9 / 2.1 63C82W (TGT→TGGycbUputative fimbriae protein
Reads supporting (aligned to +/- strand):  ref base T (24/20);  new base G (2/16);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGATATTTCAGTGCGTTATCTTGTTCTTTAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAA  >  NC_000913/1002899‑1003361
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
aaaCGGTTTTTCAGTGCGTTACTTTGTTTTTTAGCTTTTTTAACAGCAATGGCGGGATGGCTGGGCCTGAACAAGTTAGTAGGAATATTTATGGGAATGTGGTTGTTCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGGGcc                                                                                                                                                                                                                                   <  2:796323/239‑1 (MQ=255)
aaaCGATATTTCAGTGGCTTCTCTTGTTCCTTAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGcc                                                                                                                                                                                                                                   <  1:1320805/239‑1 (MQ=255)
      tatTTCAGTGCGTTATCTTGTTCTTTAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCGGCTg                                                                                                                                                                                                                               >  1:802046/1‑237 (MQ=255)
                          ttCTTTAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAGTGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATTTCGATACTACCGGTTTCGCCGCTGGTATTATCGGCGTGGAAGCCCt                                                                                                                                                                                                         >  2:515254/1‑239 (MQ=255)
                          ttCTTTAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGTCCGCGAATACCGTAAGCAATTCGTCTATTTTGAATATCGATATCCCCGGTTGGGCCCCGGGTATAACTGGGGCGGACGACCt                                                                                                                                                                                                         >  2:845767/1‑239 (MQ=255)
                          ttCTTTAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGGGCCGCTGGGATTACTGGCGCGGGCGTCCt                                                                                                                                                                                                         >  1:87761/1‑239 (MQ=255)
                               tAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCACAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGATGTGCACCGGGTATTCCTGGCGCGGACGCCCtttttc                                                                                                                                                                                                    >  2:1047517/1‑238 (MQ=255)
                                    ttCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCCCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTGGTGCCGGGGGATTTACTGGGGCGAAGGGCCTTTTGAgcggg                                                                                                                                                                                               >  2:1319212/1‑238 (MQ=255)
                                                        ttGGTTGTGCCTGAACAAGTTTGTACGAATATTTATGGTAATGTGGTCGATCTGGGCTTTGATGGCGCCACCAAAAGTGCATTCAAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTgcgc                                                                                                                                                                            <  2:657619/237‑1 (MQ=255)
                                                             tGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGaca                                                                                                                                                                       <  1:205930/238‑1 (MQ=255)
                                                                   tGAACAATTTGTTTTGTATATTTTTGGGAATGGGGTCGATCAGGGTGGTGATGTGCCCCCAAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGGGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCa                                                                                                                                                                <  2:1100263/239‑1 (MQ=255)
                                                                     aaCAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACAAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGAGGTCCTGTTTAGCGGCGAGGCTGCAACCCTTGGGCCGacaacgctacaa                                                                                                                                                              >  1:321828/1‑230 (MQ=255)
                                                                          gTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGTCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGcc                                                                                                                                                                                                                                   >  2:603367/1‑165 (MQ=255)
                                                                          gTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGcc                                                                                                                                                                                                                                   <  1:603367/165‑1 (MQ=255)
                                                                            tAGTATGCATATTTATGGGAATGGGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAaactaac                                                                                                                                                         <  1:1142230/237‑1 (MQ=255)
                                                                               tATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGaca                                                                                                                                                   >  1:641458/1‑239 (MQ=255)
                                                                                    atatTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTCCTGGCGCGGACGTCCTTTTTAGCGGCGAGGCGGAAACCCTTGCGCCTGCACTGCTCAACCTAACTTACACAGGc                                                                                                                                              >  1:154606/1‑239 (MQ=255)
                                                                                    atatTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAgcg                                                                                                                                                                                                 <  1:901792/189‑1 (MQ=255)
                                                                                    atatTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAgag                                                                                                                                                                                                 >  2:901792/1‑187 (MQ=255)
                                                                                       tttATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGATGCTTGCCCCCTTGCGCCCCCACTGCTCAAACTAACCGTCACAGTCGGa                                                                                                                                           >  1:1222309/1‑239 (MQ=255)
                                                                                         tATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTACCTGACACAGGCGGAAg                                                                                                                                         >  1:1093621/1‑239 (MQ=255)
                                                                                         tATGGGAATGGGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATAGTCATATTGGTGATTTTAATAGCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATAGCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAg                                                                                                                                         <  1:412142/239‑1 (MQ=255)
                                                                                          aTGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGc                                                                                                                                        <  1:515254/239‑1 (MQ=255)
                                                                                              gAATGTGGTCGATCAGGGCTGTGATGGCGCCACCAAAAGGTCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCggt                                                                                                                                     <  1:1239570/238‑1 (MQ=255)
                                                                                               aaTGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGACTATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAaac                                                                                                                                                             >  1:1044372/1‑214 (MQ=255)
                                                                                                 tgtgGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCCCAAACTAACTGACACAGGCGGAAGCgggggga                                                                                                                                 >  1:1257641/1‑237 (MQ=255)
                                                                                                 tggggTCGATCAGGGCTGGAATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTa                                                                                                                                 <  1:1021895/236‑1 (MQ=255)
                                                                                                      tCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGGCACCCTTGCGACGACACTGCTCAAACTAAACGCAACAGGCGGAAGCGGGGGGATGGc                                                                                                                             >  1:461626/1‑238 (MQ=255)
                                                                                                                   tgATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGcgg                                                                                                               >  1:944118/1‑237 (MQ=255)
                                                                                                                   tgATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCCGACACCCTTTCGCCGGCAATGGTCCAACTAACTGGACAAGGAGGAAGCCGGGGGATGGCAAAgggggtgggcg                                                                                                               >  2:440534/1‑232 (MQ=255)
                                                                                                                        tCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGTGCCGACACTGCTCAAACTAACTGACACAGGCGGAAACGGTGGTATGGCCACGGGGAGTGCCGTGCaa                                                                                                          >  1:428412/1‑239 (MQ=255)
                                                                                                                             accaAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATGCCGTAAGCACTGCTGCTGATTTGAATATCGATATCCCCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACGCCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGc                                                                                                                 <  1:861309/227‑1 (MQ=255)
                                                                                                                               caAAAGTGCATTACAAAATATTCATATGGGGGGTTTTAATATCAGTGATTTTCAGGCCGCGAAGGCCGTAAGCACTGCTGCTGATTTGAATATGGATATCACCGGTTTGGCCGCTGGTATTACCGGCGCGGACGGCCTTTTTAGCGGCGAGGCTGACACCCGTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTg                                                                                                   <  2:428412/239‑1 (MQ=255)
                                                                                                                                      tcATTACAAATTATTCTTATTGTTTATTTTAATATCAGTGATTTTCAGGCCGCAAATTCCGTCAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCa                                                                                            <  2:87761/238‑1 (MQ=255)
                                                                                                                                                 tatTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAACCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGTGCAAAGTCAGCAAg                                                                                 >  1:183376/1‑239 (MQ=255)
                                                                                                                                                         ttGGTTTTTTTTATTTTAGTGATTTTAAGGCCCGGAATTCCGTTAGCACGGGTGCGGGTTGGAATATCGATATTCCCGGTTGTGCCGCTGGGATTTCTGGCGCGGACGTCCTTTTTAGCGTCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCg                                                                         <  2:1222309/239‑1 (MQ=255)
                                                                                                                                                           ggTGATTTTAATATCATTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCaaa                                                                                                                                                              <  2:1031436/153‑1 (MQ=255)
                                                                                                                                                           ggTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCaaa                                                                                                                                                              >  1:1031436/1‑153 (MQ=255)
                                                                                                                                                             tGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACCGACCCAGGCGGAAGCGGGGGGCTGGCCACGGGGACTGCCGTGAAAACTCTTGAGGCGCAAAGTCAGAAAAAAAAccacccca                                                                     >  2:629242/1‑233 (MQ=255)
                                                                                                                                                                ttttATTATTTGTGGTTTTCGGGCCGAGAATACCGTAAACACTGGTGCTGGTGTGAAATTGTATAGCCCCGGGTGGTGCGACGGTTTTACGGGCGCGGACGTCCTTGGTAGCGGCGAGGCGGACACCCTTGGGCCGACACTGCTCAACCGAGCTGATACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATc                                                                  <  2:802046/239‑1 (MQ=255)
                                                                                                                                                                       atCAGTGATTTTCAGGCCGCGAATACCGTAAGCCCTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCa                                                           <  1:461816/239‑1 (MQ=255)
                                                                                                                                                                        tCAGTGATTTTCAGGCCTCGAATACCGGAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAg                                                          <  1:531978/239‑1 (MQ=255)
                                                                                                                                                                                 ttCGGGCCGCGAATACCGTAAGCCCTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACg                                                 <  1:440534/239‑1 (MQ=255)
                                                                                                                                                                                 ttCAGGCCCGGAAGACCCTAAGCACTGCTGCGGATTTGAATATCGATTTCACCGGGTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACg                                                 <  1:1248864/239‑1 (MQ=255)
                                                                                                                                                                                  tCAGGCCGCGAATTCCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGc                                                <  1:318624/239‑1 (MQ=255)
                                                                                                                                                                                          cgcaTACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCCACACTGCTCACACTTACTGACACAGGCGGGAGCGGGGGGCTGGTCACGGGGGATGCAGTTGCACTTCTTCGAGCGCAAAAACCGCCAGAAAACCCGCTCAATGAAGTCCAGCCCCTTAACCCCATaaaa                                        >  2:7640/4‑239 (MQ=255)
                                                                                                                                                                                             atgccGTAATCGCTGCTGCTGATTTGAATGTCGATATCCCCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGcc                                     <  1:845767/236‑1 (MQ=255)
                                                                                                                                                                                                ccGTAAGCCCTGCTGCTGATTTGAATTGCGATATCCCCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTc                                                                                       <  2:269833/186‑1 (MQ=255)
                                                                                                                                                                                                        cctgctgctgATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAaca                           <  1:102975/237‑1 (MQ=255)
                                                                                                                                                                                                        aCTGCTGCTTATTTGAATTGCGATATCACCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTGTAGCGGCGAGGCTGAACCCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAaca                           <  1:277537/238‑1 (MQ=255)
                                                                                                                                                                                                        aCTGCTGCTGATTTGAATATCGATATCACCGGTTGGGCCGCTGGTATTACTGGCCCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAaca                           <  1:310535/238‑1 (MQ=255)
                                                                                                                                                                                                         ctgctgctgATTTGAATATCGATATCTCCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAacaca                         <  1:1319212/239‑1 (MQ=255)
                                                                                                                                                                                                          tgctgctgATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGGCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAacaca                         >  1:335400/1‑238 (MQ=255)
                                                                                                                                                                                                          tgctgctgATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGGGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAAt                                                                             >  1:269833/1‑186 (MQ=255)
                                                                                                                                                                                                           gctgttgATTGGAATCTGCATTTCACCGGTTGGTGCGCGGGTAGTACGGGGCGGGGCCTCCTTTTTAGCGGCGAGGATGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACt                       <  2:944118/239‑1 (MQ=255)
                                                                                                                                                                                                           gctgctgATTTGAATTTGGATATCCCCGGTTGTGCCGCTGGTATTACGGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACt                       <  1:1124718/239‑1 (MQ=255)
                                                                                                                                                                                                             tgcgggtttGAATATCGATATCACCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCa                     <  1:1047517/233‑1 (MQ=255)
                                                                                                                                                                                                               ctgATTTTAATACTGATAGCCCCGGTTGGGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCa                                                                                            <  2:690707/166‑1 (MQ=255)
                                                                                                                                                                                                                      taaTTTCGATATCGCCGGTTGGCCCGCTGGTATTGCGGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCTTTGCGCCGACACGGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAAc            <  2:154606/238‑1 (MQ=255)
                                                                                                                                                                                                                         tatCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGAAACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCCTGCAAATTCTTGATGCGCAAAGTCAGcca                                                                                  >  1:690707/1‑164 (MQ=255)
                                                                                                                                                                                                                             gATATCCCCGGTTGGGCCGCTGGTAATACTGGCCCTGACGTCCGTTTTAGCGGCGAGGCTGGCGCCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTta     <  1:608417/239‑1 (MQ=255)
                                                                                                                                                                                                                              atatCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGACAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACCCAAAAATCAACTTCGTtat    >  1:1134773/1‑239 (MQ=255)
                                                                                                                                                                                                                                ttcccaggTTGTGCCGTTGGTATTACTGGGGCGGTCTTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATaa  <  1:523062/233‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
AAACGATATTTCAGTGCGTTATCTTGTTCTTTAGCATTCTTAACATCCATGTCGGGATGGCTGGGCCTGAACAAGTTAGTATGCATATTTATGGGAATGTGGTCGATCAGGGCTGTGATGTCGCCACCAAAAGTGCATTACAAAATATTCATATTGGTGATTTTAATATCAGTGATTTTCAGGCCGCGAATACCGTAAGCACTGCTGCTGATTTGAATATCGATATCACCGGTTGTGCCGCTGGTATTACTGGCGCGGACGTCCTTTTTAGCGGCGAGGCTGACACCCTTGCGCCGACACTGCTCAAACTAACTGACACAGGCGGAAGCGGTGGTATGGCAACGGGGATTGCCGTGCAAATTCTTGATGCGCAAAGTCAGCAAGAAATCCCGCTCAATCAGGTCCAGCCTCTTACGCCCTTAAAAGCCGGGGATAACACACTCAAATATCAACTTCGTTATAA  >  NC_000913/1002899‑1003361

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