Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 12,026 T→G intergenic (‑240/‑137) yaaI ← / → dnaK DUF2541 domain‑containing protein YaaI/chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,0260TG97.7% 131.9 / ‑4.8 43intergenic (‑240/‑137)yaaI/dnaKDUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (26/16);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATC  >  NC_000913/11779‑12262
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
gATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGACTGCGATTCATT‑CTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAA‑TTTTTGCATCTCCCCCTTGATGACGTGGTTTAc                                                                                                                                                                                                                                          >  1:16881/1‑251 (MQ=255)
  tttCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACg                                                                                                                                                                                                                                         >  1:198251/1‑251 (MQ=255)
      aTAAGGATTTTCTTAGTGGGAAGAGGTAGAGGGAAGAATACCCACTAGTTTACTGCGGATAAAGAGAATTGTCAGTTACGTAATCTTTTCTTTTTATTACAATTTTTTGAGGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGCCAATTTTACTTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCGCCCCCTTGATGACGTGGGTTACGAccc                                                                                                                                                                                                                                     <  2:301696/251‑1 (MQ=255)
                               gtatagggATTAATACCCACTAGTTTACTTCTGAAAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCagtgagt                                                                                                                                                                                                            <  1:405988/246‑1 (MQ=255)
                                  gagggaggaaacccccttagagtACTCCTGATAAATAGAAGGTTCTGTTACGTAACCTTTCCGTTTTATTACAATTTTTTGAGGAATGACTCGGCTGTGATTCCTT‑CTTTTTCTGAATAAAATTGTGGTCATTTTTCCGTATTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCTCCTTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGtgagtgag                                                                                                                                                                                                         <  2:141778/229‑1 (MQ=255)
                                      gATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTg                                                                                                                                                                                                     >  1:70318/1‑251 (MQ=255)
                                                               gataaagagaaggttagggtactgaaTCTTTTCTTTTTATTACAATTTTTTGATTAATGCTTCGGCTGTGCATCATTTCTTT‑TCTGAATAAAATTGCTGTCAATTTTACGCTTTGTCCGGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAAGTGTTGTATCTCCACCTTGATGCCGTGGGTTCCGACCCCATTTAGTAGTTAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTg                                                                                                                                                                             <  1:39028/226‑1 (MQ=255)
                                                                      agaaggTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCAtt                                                                                                                                                                                                                                 <  2:53886/191‑1 (MQ=255)
                                                                      agaaggTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCAtt                                                                                                                                                                                                                                 >  1:53886/1‑191 (MQ=255)
                                                                        aaggttCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACg                                                                                                                                                                   <  1:259878/247‑1 (MQ=255)
                                                                                                   ttattaCAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc                                                                                                                                        >  1:204541/1‑251 (MQ=255)
                                                                                                      ttaCAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcac                                                                                                                                      >  1:395683/1‑250 (MQ=255)
                                                                                                          aaTTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatgat                                                                                                                                 >  2:217678/1‑251 (MQ=255)
                                                                                                             tttttGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAGAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAACATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCCCAACCACatgatgcc                                                                                                                              >  1:66905/1‑249 (MQ=255)
                                                                                                                ttGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCgc                                                                                                                           >  1:332322/1‑250 (MQ=255)
                                                                                                                       gtGCCTTGGCGGCGATTCATT‑CTTTATATGAATACATTTGCTGTCAATTTTACGTCTTGTCCTGCCATTTCGAGACATTTCTGCGCAAAAGAACAAAAAATTTTTGCATCTCCCCCTTGCTGTCGTGGGTTACGACCCAATTTAGTAGTCAACCGCAGTGAGTGAGTCTGAAAAAACATGAAATTGGGCAGTTGAACCCAGACGTTTCGACCATATTACAAACTCAAAACCACATTATGACCGAatatat                                                                                                                     >  2:49330/2‑250 (MQ=255)
                                                                                                                          ccTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTg                                                                                                                 >  1:322606/1‑251 (MQ=255)
                                                                                                                               gcgtcGATTCATT‑CTTTATCTGAATAAAATTGCGGTCAAGTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTGTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGAc                                                                                                            <  2:292306/247‑1 (MQ=255)
                                                                                                                                        catt‑cTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGAt                                                                                                    <  1:293793/250‑1 (MQ=255)
                                                                                                                                         ccttctttttcggAATAAAATTGTTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAAGTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATg                                                                                                   <  1:253303/239‑1 (MQ=255)
                                                                                                                                                  tttGAATAAAATTGCGTTCAAGTTTCCGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAaat                                                                                           <  2:481609/248‑1 (MQ=255)
                                                                                                                                                         aaaaTTGCTGTCAATTTTACGTTGTTCCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTa                                                                                   <  2:70318/251‑1 (MQ=255)
                                                                                                                                                         aaaaTTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGCAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATTATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGt                                                                                    >  2:502702/1‑250 (MQ=255)
                                                                                                                                                                cTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCt                                                                            >  2:58575/1‑251 (MQ=255)
                                                                                                                                                                               cTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGCCGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGACTATATAGTGGAAACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAActc                                                             >  1:82162/1‑251 (MQ=255)
                                                                                                                                                                                         cATATCGCGAAATTTCTGCGCAAACGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCg                                                   >  1:200040/1‑251 (MQ=255)
                                                                                                                                                                                          atatCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCg                                                   >  2:156198/1‑250 (MQ=255)
                                                                                                                                                                                           tatCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAACTAATTGGTATCGACCTGAGTACTACCAACTCTTGTGTAGCGAt                                                 >  2:253051/1‑251 (MQ=255)
                                                                                                                                                                                                  aaaTTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCACATATAGTAGTCAACCGCAGTGAGGGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGCATATATAGTGGAGACGTTAAGATGGGTAAAATAATTGGTATAGACCTGGGTACTACAAACTCTTTTGTAGCGATAatggat                                          >  2:17325/1‑251 (MQ=255)
                                                                                                                                                                                                   aaTTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGAGGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTatggatg                                         >  1:466947/1‑251 (MQ=255)
                                                                                                                                                                                                          gcgcAAAAGCACAAAAAATTTATGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATAGAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGcacc                                    <  1:95912/249‑1 (MQ=255)
                                                                                                                                                                                                           cgcAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTAACAACTCTTGTGTAGCGATTATGGATGGcaccacc                                 >  1:229265/1‑250 (MQ=255)
                                                                                                                                                                                                             caacaacgcaAACAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTAAAGACTCACAACCACATGATGACCGAATATAAAGTGGAGCCGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACAAACTCTTGTGTCGCGATTATGGATGGCACCActcc                               >  2:418961/9‑251 (MQ=255)
                                                                                                                                                                                                             cAAAAGCCCAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCATTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTTCCAACTCTTGTGTAGCGATTATGGATGGcacctct                                 <  1:278196/249‑4 (MQ=255)
                                                                                                                                                                                                             cAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCActcc                               >  2:218885/1‑251 (MQ=255)
                                                                                                                                                                                                                aaGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTcg                            >  1:359977/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTcgtg                          <  1:371128/251‑3 (MQ=255)
                                                                                                                                                                                                                      aaaaaaTTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGttt                                                                                                        <  2:16230/169‑1 (MQ=255)
                                                                                                                                                                                                                      aaaaaaTTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGttt                                                                                                        >  1:16230/1‑169 (MQ=255)
                                                                                                                                                                                                                       aaaaaTTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTg                     >  2:118848/1‑251 (MQ=255)
                                                                                                                                                                                                                                tGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCACCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACCGACTCACAACCACATGCTGACCGAATATATAGTGGAGACGTTTAAATGGGTAAAATAATTGGTATCGACCTGGGTCCTACCCACTCTTGTGAAGCGATTATGGATGGCACCCCTCATCGTGTGCTGGGGAACGcc            >  2:308310/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aTCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGa          >  2:26822/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ccTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATc  <  1:58575/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
GATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATC  >  NC_000913/11779‑12262

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