Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,412,898 C→T 11.4% A190T (GCC→ACC)  acrS ← acrAB operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,412,8980CT11.4% 179.3 / 4.7 66A190T (GCC→ACC) acrSacrAB operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (46/12);  new base T (2/6);  total (48/18)
Fisher's exact test for biased strand distribution p-value = 3.92e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTTGATACTTTACTGATGCAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCATCTCATCATTAAATTCACATTTGTGATATA  >  NC_000913/3412669‑3413128
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       
cTGTTTTATTTTTTACTGATTCAGTGGATTGAACTGCTGGCAACAAAAACAGAAAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtata                                                                                                                                                                                                                                <  2:706008/239‑1 (MQ=255)
cTGTTTGATACTTTACTGATTCAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtata                                                                                                                                                                                                                                <  1:943684/239‑1 (MQ=255)
cTGGTTGATACTTTACTGATGCAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtata                                                                                                                                                                                                                                <  1:735381/239‑1 (MQ=255)
 tttttgtttattttCTGATGTAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGTTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATaa                                                                                                                                                                                                                               <  2:477933/229‑1 (MQ=255)
 tGTTTGATTCTTTTCTGATGCAGTGGATTTAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATaa                                                                                                                                                                                                                               <  2:5651/239‑1 (MQ=255)
     tGATACTTTACTGATGCAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGAt                                                                                                                                                                                                                           >  1:309194/1‑239 (MQ=255)
             tACTGATGCAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTAAATGAAACTTAATTCATCCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCTTGAACATTCTTAAAACGTTATCGCCCATCGCGGGGGTTTGTTTATAAGGATCATatac                                                                                                                                                                                                                    >  2:223052/1‑235 (MQ=255)
                  aTGCAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCat                                                                                                                                                                                                              >  1:500019/1‑239 (MQ=255)
                          gATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATaaa                                                                                                                                                                                                                              <  1:808454/215‑1 (MQ=255)
                             tGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCCTTAAcc                                                                                                                                                                                                   >  1:78516/1‑239 (MQ=255)
                                 ctgctgGAAAAAAAAAAAGAAAAAGGCCGAAATATGAACTCAATTAACCCAATAAACTTGTTTACTGTGGACAACCTAAATTTTCAGAACAAATTGTAATTTTACATGCCCCTAAATTCATTCGTTTGATAAATTAATTTTTTTATTTTTTCATCGGGAATAAACATTCTTAATACGTTACCGACAAGGGCGGGGGTTTGTTTAAAAAGATCATAACCCGCCATATTCATTAACCAGtt                                                                                                                                                                                               <  1:68269/239‑1 (MQ=255)
                                    ctgGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTg                                                                                                                                                                                            >  1:185447/1‑239 (MQ=255)
                                                     cAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAg                                                                                                                                                                            >  1:630455/1‑238 (MQ=255)
                                                      aaaTGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCCCTGAAgg                                                                                                                                                                           >  1:561464/1‑238 (MQ=255)
                                                            cgaaatagaaaCTCAGTAAACCCACAAAAATTGGTTAATGGTGACAATTTAAATTTTAAGGAAAGAATGTGAATTAACATGACTCTAAATTCATTCTTTTGATAAATTAATTTCTTTATTTTTTCATCTGGCATGAACATTTTTAATACTTTATCGCCCAGAGGGGGGGTTTGTTAATAAAGATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCAt                                                                                                                                                                     <  1:431266/229‑1 (MQ=255)
                                                              cAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACTTGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATcca                                                                                                                                                                  >  1:656812/1‑237 (MQ=255)
                                                              cAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGAAAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTATTACGTTTTCGACCAGATCGGGGGCTGGTTTATAAAGATCATAACCCGCAATATTCATTAACCAGTTTTGAACAAATCGGCTGAAGACACCACCaa                                                                                                                                                                  >  2:857164/1‑239 (MQ=255)
                                                              cAATATGAACTCAGTTAATCCACTAAAACTGGTTAACGGTGACGAACTGACTTTTCAGGCCAGATTGCGAATTTACGTGACACTTAATTAATTAGTTTGATGAATTAATTTCGTTATGTTTTCATATGGCCTGCACATTCTTAATACGTTATCGACCCGACCGGGGGTTTGTTTATAAAGATCTTAACCCGCCATTTTCATTAACCAGCTTTGACCCATACTCCTGCAGGCACCATcgc                                                                                                                                                                  >  1:149383/1‑237 (MQ=255)
                                                                         cAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACAAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATTTTCATTAACCAGTTTTGACCAATTCCGCTGAAGGCACCATCATTATTAaccatc                                                                                                                                                      >  1:1044157/1‑239 (MQ=255)
                                                                           gTTAACCCAAAAAATCTGTTAACTGGTGAGGAACTGATTTTGCGGGAGAGATTGTGTATGTAGAGGACATTAAATTCATTTGTTTGATAAATTATTTTGTTTATTTTTTCATCTGGCATGAACATTCTTATTATTTTACCGAGCAGAGGGGGGGTTTTTTTATAAAGATCAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATcac                                                                                                                                                    <  1:452138/239‑1 (MQ=255)
                                                                                    cTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCt                                                                                                                                                                                                                                        <  2:518168/147‑1 (MQ=255)
                                                                                    cTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCt                                                                                                                                                                                                                                        >  1:518168/1‑147 (MQ=255)
                                                                                    cTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTaa                                                                                                                                           >  1:143818/1‑239 (MQ=255)
                                                                                    cTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCAGCTGAAGGCACCATCAATAATAATCATAACAACATCTaa                                                                                                                                           >  1:764177/1‑239 (MQ=255)
                                                                                    cTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTACGTAATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACAAGTTTTGAACAATTCCGCTGAAGGAACCATCACTAATAATAATCacaacaccaca                                                                                                                                           >  1:702001/1‑234 (MQ=255)
                                                                                       aaaTTGTTTAACTGGAACAAATTGAATTTTAAGGACAGATTGTAATTTTACTTGCCTATTAATTCATTTGTTTGATGAATTAATTTGGTTATTTTTTAATCTGGCATGAACATTCTTAATCCGTTACCGCCCAGAGCGGGGGTTTGTTTAAAAAGATCAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATc                                                                                                                                        <  1:857164/239‑1 (MQ=255)
                                                                                          cTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAatcatccaaa                                                                                                                                                  >  1:940880/1‑222 (MQ=255)
                                                                                           tGGTTAATGGTGACAAACTGAATTGTCAGACAAGTATTTTATTTTACATGTCTTCTAATTTATTGTTTGTAATAATTATTTTGTTTATTTTTCAATCTGGCATAAACATTTTTAATCCGTTATCGACGAGGCGGGGTGTTTTTTTATAAAAATCAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAgg                                                                                                                                    <  1:724425/239‑1 (MQ=255)
                                                                                             gTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATaa                                                                                                                                                                                                                               >  2:471945/1‑147 (MQ=255)
                                                                                             gTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATaa                                                                                                                                                                                                                               <  1:471945/147‑1 (MQ=255)
                                                                                             gTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATATAGATCATAACCCGCCCTATTCATTAACCAGTTTTGAACAATTCCGATGAAGGCCCCATAAAAAATAATACAACCAACACCAAAATCGAGGta                                                                                                                                  >  2:299350/1‑238 (MQ=255)
                                                                                             gTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCTTATTCATTAACCAGTTTTTAACAATTCCGCTCAAGCCACCATCAATAAAATTCAACACACAATCTAAATAGGGGta                                                                                                                                  >  2:22017/1‑238 (MQ=255)
                                                                                              ttAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCTTATTCATTAACCAGTTTTGaa                                                                                                                                                                                          >  2:1092526/1‑183 (MQ=255)
                                                                                                  aaggggacagAAATGAATTTACAGGAAAAAATGTATATTTACATGACTCTTAATTCCTTGTTTTGATGAATTAATTTCTTTATTTTTTCATGTGGAAAGAACATTCTTAATACGTTATCGACCAGAGGGGGGGCTTGTTTATAAAGATCAAACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTAtt                                                                                                                               <  1:479695/230‑1 (MQ=255)
                                                                                                         aCGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtata                                                                                                                                                                                                                                <  2:532891/136‑1 (MQ=255)
                                                                                                         aCGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtata                                                                                                                                                                                                                                >  1:532891/1‑136 (MQ=255)
                                                                                                         aCGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGAACCATCAATAATAATCATCACACCATCAAAACCGAGGTTATTTGcacaa                                                                                                                        >  2:908825/1‑235 (MQ=255)
                                                                                                                 aaTTTTCAGGACAGAATGCGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAAAAATCCTCACAAAATCTAAATCAAGGTTATTTGATACaaaa                                                                                                                     >  2:479695/1‑234 (MQ=255)
                                                                                                                         ggACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCCCCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTAAACAACATTGTTACTCac                                                                                                        >  2:494027/1‑239 (MQ=255)
                                                                                                                                         tACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGt                                                                                                                                                                                                                                     >  1:770272/1‑99 (MQ=255)
                                                                                                                                         tACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGt                                                                                                                                                                                                                                     <  2:770272/99‑1 (MQ=255)
                                                                                                                                                          attcGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATACAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCATTAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACACACTTGTTGCTGACACCCCTGCAAATCTTCCCGGCGATTCTTCGGAt                                                                       >  2:1029153/1‑239 (MQ=255)
                                                                                                                                                           ttcGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGGGAGTCTGCGGGtt                                                                      >  1:931989/1‑239 (MQ=255)
                                                                                                                                                           ttcGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGAtt                                                                      >  1:879123/1‑239 (MQ=255)
                                                                                                                                                            tcGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtata                                                                                                                                                                                                                                >  2:661025/1‑85 (MQ=255)
                                                                                                                                                            tcGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtata                                                                                                                                                                                                                                <  1:661025/85‑1 (MQ=255)
                                                                                                                                                            tcGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCCTCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTTTTGCTGACACGCCTGCAATACTTCCCGGAAAATCTCCGGAtaa                                                                     >  2:523767/1‑237 (MQ=255)
                                                                                                                                                            tcGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTa                                                                     >  1:288552/1‑239 (MQ=255)
                                                                                                                                                                              ttCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGAAACGCCTGCAATACTTCGCGGCAAGTCTGCGGATTAAAGCCCATCTTTTcgcg                                                    >  1:370643/1‑239 (MQ=255)
                                                                                                                                                                                                      aaCATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCAtctc                            >  1:671709/1‑239 (MQ=255)
                                                                                                                                                                                                      aaCATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATAACTCCCTCGGCCAGCAtctc                            >  1:371788/1‑239 (MQ=255)
                                                                                                                                                                                                           tCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACCCAACCTTGTTGCTGCAACGCCTGCAATACTTCCCGGGGAGTCTGAGGATTACACCCCACCTTTTCGCGTTTCACTCCCCCGGCCGGCCACTcatcat                       >  2:317727/1‑239 (MQ=255)
                                                                                                                                                                                                           tCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCACCtcat                          >  1:232987/1‑236 (MQ=255)
                                                                                                                                                                                                                     ttATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATcaca                                                                                                                                                   >  1:445130/1‑105 (MQ=255)
                                                                                                                                                                                                                     ttATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATa                                                                                            >  1:225829/1‑160 (MQ=255)
                                                                                                                                                                                                                     ttATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAg                                                                  >  1:943593/1‑186 (MQ=255)
                                                                                                                                                                                                                     ttATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCATCTCATCATTAAATTcaca             >  1:184263/1‑239 (MQ=255)
                                                                                                                                                                                                                      tATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCATCTCATCATTAAATTCACAt            >  1:994642/1‑239 (MQ=255)
                                                                                                                                                                                                                       aTCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGTACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGCGAGTCTGCGGATTAAAGCCCCTCTTTTCGCGTATCACACCCTCTGCCAGCAGCTACTCATTAAATTACCAtt           >  2:1063100/1‑239 (MQ=255)
                                                                                                                                                                                                                       aTCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACACCCTTGTTGCTGACACGCCTGCAATCCTTCGCGGACAGTCTCCGGATTAAATCCCATCTTTTCGCGTATCACCCCCTCGGCACCCCTCTCATCATAAAATTCACAtt           >  2:906432/1‑239 (MQ=255)
                                                                                                                                                                                                                       aTCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATa                                                                                            >  1:741241/1‑158 (MQ=255)
                                                                                                                                                                                                                       aTCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCAGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTACGGATTAAAGCCCATCTTTTCCCGTATCACTCCTTCGGCCAGCATCTCATCATTAAATTCACCtt           >  1:513774/1‑239 (MQ=255)
                                                                                                                                                                                                                               gagCGGGGGCTTGGTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACACTTCCGCTGTCGGCACCATCACTAATAATCATCACAACATCTAAATCGAGGTTATTTCCTACACAACCTTGTTGCTGACACCCCTGCAATACTTCCCGGTGATCCTGCGAATTAAAGCCCATCTTTTCCAGTATCCCTCCCCCGCCACACATCTCATCATTAAATTCACAATTGTGaaaa   >  2:460031/1‑236 (MQ=255)
                                                                                                                                                                                                                                agCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCATCTCATCATTAAATTCACATTTGTGAtata  >  1:718209/1‑239 (MQ=255)
                                                                                                                                                                                                                                agCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTACGCGGAGAGTCTGCGGATTAAAACCCATCTTTTCGCGAATTACTCCCTCGGCCAGCATCTCATCATTAAATTCACATTtgtgtaata  >  1:502102/1‑234 (MQ=255)
                                                                                                                                                                                                                                 gggggggTTTTGTTAAAAAAAACATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCATCTCATCATTAAATTCACATTTGTGAtata  <  1:317727/236‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                       
CTGTTTGATACTTTACTGATGCAGTGGATTGAACTGCTGGCAACAGAAACAGACAAATGCCGCAATATGAACTCAGTTAATCCACTAAAACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAGTCTGCGGATTAAAGCCCATCTTTTCGCGTATCACTCCCTCGGCCAGCATCTCATCATTAAATTCACATTTGTGATATA  >  NC_000913/3412669‑3413128

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: