Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,367,510 T→G 10.6% S135R (AGC→CGC)  pspF ← psp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,367,5100TG10.6% 198.7 / 12.7 66S135R (AGC→CGC) pspFpsp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base G (5/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGT  >  NC_000913/1367362‑1367657
                                                                                                                                                    |                                                                                                                                                   
cTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCt                                                                                                                                                     <  3:213422/149‑1 (MQ=255)
 tttCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAGAGCGAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGTCTTGATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACCTTCACCTGCAATCGTTGGCTg                                                                                                                                                    >  3:545462/1‑149 (MQ=255)
           gcgcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAAcgc                                                                                                                                          >  4:433194/1‑149 (MQ=255)
           gcgcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCGGCAGCCAAcca                                                                                                                                          >  4:795669/1‑147 (MQ=255)
            cgcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCGGCAGCCAAccat                                                                                                                                         <  2:369030/149‑4 (MQ=255)
             gcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACgcgc                                                                                                                                        >  2:137782/1‑149 (MQ=255)
               aGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCGGCAGCCAAccattgc                                                                                                                                      <  3:795669/149‑7 (MQ=255)
                    tggCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGtt                                                                                                                                 >  1:59615/1‑149 (MQ=255)
                      gCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCa                                                                                                                               >  4:76074/1‑149 (MQ=255)
                      gCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCa                                                                                                                               <  2:745947/149‑1 (MQ=255)
                         gTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCg                                                                                                                            <  4:819743/149‑1 (MQ=255)
                           tGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTa                                                                                                                          <  1:889251/149‑1 (MQ=255)
                           tGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTa                                                                                                                          <  2:195958/149‑1 (MQ=255)
                             tACAACAGCAAAAGCCCGTCGGTTGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACTATCGCCGGGAGATCGGCATTCGTCGAGCATACCAACCGCACATTCACCTGCAAGGGTTGGCTGCTGCCAACGCGCTCCAGTTCACCGTACt                                                                                                                        >  4:848456/1‑149 (MQ=255)
                                  caTCAAAAGACAGTCAGTCGAGCAGGTCAACGCGAAAAGTGCCTTCATTGACCATCGCCGGGAAATCGGCCTACGTCGGCCATACCAACCGCACATTCACCTGCAATGTTTGACTACCGCCAACGCGCCCCAGTTCACCGTACTCAATc                                                                                                                   >  4:269903/1‑149 (MQ=255)
                                    tCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAc                                                                                                                 >  4:171155/1‑149 (MQ=255)
                                      aaaaGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                               <  3:521064/149‑1 (MQ=255)
                                      aaaaGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                               <  3:691516/149‑1 (MQ=255)
                                       aaaGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgcg                                                                                                              >  2:432274/1‑149 (MQ=255)
                                          gCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCaa                                                                                                           <  2:123585/149‑1 (MQ=255)
                                               tcggtcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAAtt                                                                                                      <  2:754542/149‑1 (MQ=255)
                                                  gtcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTtct                                                                                                   <  3:66477/149‑1 (MQ=255)
                                                   tcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCCCCGTACTCACTCCCGCGCAATAATTcctc                                                                                                  <  4:459031/149‑1 (MQ=255)
                                                   tcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCACCGGGAGATCGGCATCCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTtctc                                                                                                  <  4:346236/149‑1 (MQ=255)
                                                     gAGCAGGTCAGTGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCt                                                                                                >  2:781348/1‑149 (MQ=255)
                                                               gcgcGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACcatcat                                                                                      <  2:2148/149‑1 (MQ=255)
                                                                cgcgAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCatcatc                                                                                     <  2:74841/149‑1 (MQ=255)
                                                                 gcgAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCg                                                                                    >  2:829238/1‑149 (MQ=255)
                                                                  cgAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCgg                                                                                   >  2:850384/1‑149 (MQ=255)
                                                                     aaaGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGc                                                                                <  1:833579/149‑1 (MQ=255)
                                                                      aaGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGcc                                                                               <  3:273238/149‑1 (MQ=255)
                                                                         tGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                               <  3:399009/114‑1 (MQ=255)
                                                                         tGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                               >  4:399009/1‑114 (MQ=255)
                                                                           ccTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGc                                                                          >  1:654391/1‑149 (MQ=255)
                                                                             ttCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCGATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGa                                                                        >  4:508915/1‑149 (MQ=255)
                                                                                aTTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGtt                                                                     >  4:98298/1‑149 (MQ=255)
                                                                                      cATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCaa                                                               >  4:183487/1‑149 (MQ=255)
                                                                                       aTCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAg                                                              >  4:813075/1‑149 (MQ=255)
                                                                                        tcgtcgGGAGATCGGCATTCGCCGCGCATACCAACCGCACATTCCCCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCAGCGGTGCCGTAGCGAGTTCATCAAGa                                                             <  4:241699/149‑1 (MQ=255)
                                                                                         cGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGaa                                                            <  4:707583/149‑1 (MQ=255)
                                                                                           ccGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAAt                                                          <  1:468543/149‑1 (MQ=255)
                                                                                              ggAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGc                                                       <  4:327239/149‑1 (MQ=255)
                                                                                                  aTCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTAc                                                   >  2:139891/1‑149 (MQ=255)
                                                                                                   tCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTAcc                                                  >  1:154250/1‑149 (MQ=255)
                                                                                                                     aCCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGc                                                                                >  3:715509/1‑101 (MQ=255)
                                                                                                                     aCCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGc                                                                                <  4:715509/101‑1 (MQ=255)
                                                                                                                     aCCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCaaa                                <  3:813075/149‑1 (MQ=255)
                                                                                                                      ccAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAt                               <  3:653029/149‑1 (MQ=255)
                                                                                                                      ccAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAt                               >  2:149799/1‑149 (MQ=255)
                                                                                                                       cAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAtc                              >  3:104969/1‑149 (MQ=255)
                                                                                                                       cAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAtc                              >  1:685357/1‑149 (MQ=255)
                                                                                                                            gCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCt                         <  4:104969/149‑1 (MQ=255)
                                                                                                                             cacaTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCt                         >  2:269748/1‑148 (MQ=255)
                                                                                                                             cacaTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCt                         <  1:269748/148‑1 (MQ=255)
                                                                                                                             cacaTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTg                        <  3:183487/149‑1 (MQ=255)
                                                                                                                                aTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGAt                     >  1:765336/1‑149 (MQ=255)
                                                                                                                                    aCCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACg                 >  4:640696/1‑149 (MQ=255)
                                                                                                                                      cTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGtt               <  3:98298/149‑1 (MQ=255)
                                                                                                                                         agcgcgTTGGCGGCAGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGttttt            >  3:391438/6‑149 (MQ=255)
                                                                                                                                          gcgcgTTGGCGGCAGCCCACGCGCTACAGTTCCCCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGATCAAATCTCCCTGGATGACGTTTTTg           >  2:34803/5‑149 (MQ=255)
                                                                                                                                          gcgcgTTGGCGGCAGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTg           >  4:77911/5‑149 (MQ=255)
                                                                                                                                          aaTGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTg           >  4:206862/1‑149 (MQ=255)
                                                                                                                                              gTTGGCGGCAGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCa       >  1:281978/1‑149 (MQ=255)
                                                                                                                                               ttGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCAc      >  4:109323/1‑149 (MQ=255)
                                                                                                                                                  gCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCTCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTcc                           >  3:576723/1‑125 (MQ=255)
                                                                                                                                                  gCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTcc                           <  4:576723/125‑1 (MQ=255)
                                                                                                                                                   cTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGt  <  1:137782/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTTCGCGCTCGCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGT  >  NC_000913/1367362‑1367657

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