Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,268 T→G 8.0% I27S (ATC→AGC) ‡ yegP → UPF0339 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,2680TG 8.0% 288.9 / 14.0 100I27S (ATC→AGC) ‡yegPUPF0339 family protein
Reads supporting (aligned to +/- strand):  ref base T (55/37);  new base G (5/3);  total (60/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

ACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGCCGCTAATCATCAA  >  NC_000913/2165123‑2165416
                                                                                                                                                 |                                                                                                                                                    
aCTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTc                                                                                                                                                  >  2:854463/1‑150 (MQ=255)
aCTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTc                                                                                                                                                  <  2:387809/150‑1 (MQ=255)
  tttAATCGGTAAACAGCACTTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCAc                                                                                                                                                <  2:739474/150‑1 (MQ=255)
  tttAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCAc                                                                                                                                                <  1:1238247/150‑1 (MQ=255)
      aTCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTTATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGGATAGTCTCACCatt                                                                                                                                            <  2:295893/150‑3 (MQ=255)
       tCGGAAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGGATAGTCTCACCattg                                                                                                                                           <  2:28024/150‑4 (MQ=255)
       tCGGAAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGGATAGTCTCACCattg                                                                                                                                           <  2:537509/150‑4 (MQ=255)
       atcgaaaaCAGCATCTAAAGGTGGTGTTTTCGATACACAATTTCAATCAAGGAGTCATTTTGGCTGGTGGGTTTGAACTCAGTAAGAGCAGTGATAATCAGGTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTGTCCTCCCCcgcg                                                                                                                                           >  2:1049880/6‑150 (MQ=255)
             aaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGCGCAGGGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCttt                                                                                                                                     >  2:944858/1‑150 (MQ=255)
               cAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGTGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTtat                                                                                                                                   >  2:338168/1‑150 (MQ=255)
               cAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTtat                                                                                                                                   >  2:1514696/1‑150 (MQ=255)
                  cACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATAcc                                                                                                                                >  2:1308937/1‑150 (MQ=255)
                  cACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATAcc                                                                                                                                <  2:1673623/150‑1 (MQ=255)
                  cACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATAcc                                                                                                                                <  1:1200926/150‑1 (MQ=255)
                  cACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATAcc                                                                                                                                <  1:1662623/150‑1 (MQ=255)
                           aTGCTGTTTTCGATACACAATTTCAATCAACGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAAcc                                                                                                                       <  2:86964/150‑1 (MQ=255)
                             gCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATTCCTCAAAAACctg                                                                                                                     >  1:250142/1‑149 (MQ=255)
                               tGTTTTCGTTACAGATTTTCTATGTAGGAGTCATTATGGCTGGTTTGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCATTGGTGAGACTTTCCTCAGCAGCGAGCTTTATACCTCAAAAACCTCTg                                                                                                                   <  2:1385154/150‑1 (MQ=255)
                               tGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTg                                                                                                                   <  1:172429/150‑1 (MQ=255)
                                gTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGc                                                                                                                  >  2:1427004/1‑150 (MQ=255)
                                gTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCATTATACCTCAAAAACCTCTGc                                                                                                                  >  2:1509501/1‑150 (MQ=255)
                                gTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAAACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGc                                                                                                                  >  2:822351/1‑150 (MQ=255)
                                  tttCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCgg                                                                                                                <  2:1418479/150‑1 (MQ=255)
                                  tttCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCgg                                                                                                                <  2:1306053/150‑1 (MQ=255)
                                     cGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGaaa                                                                                                             >  2:1669484/1‑150 (MQ=255)
                                         acacAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAggg                                                                                                         <  1:149295/150‑1 (MQ=255)
                                          cacaATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGc                                                                                                        >  2:1533086/1‑150 (MQ=255)
                                          cacaATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGc                                                                                                        >  1:1608378/1‑150 (MQ=255)
                                          cacaATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGCGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGc                                                                                                        >  2:1646457/1‑150 (MQ=255)
                                            caATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCAt                                                                                                      >  2:291706/1‑150 (MQ=255)
                                            caATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCAt                                                                                                      >  1:556878/1‑150 (MQ=255)
                                              aTTTCAATCAAGGAGTCATTATGGGTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATcg                                                                                                    <  1:306473/150‑1 (MQ=255)
                                                ttCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATcgcg                                                                                                  <  1:1669484/150‑1 (MQ=255)
                                                  caatcaaGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTc                                                                                                >  2:676203/1‑150 (MQ=255)
                                                  caatcaaGGAGGCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTGTGTGCTAACAGCGGGCCATGGTGAGACTATCATCACCAGCGCGGTTTATACGTCACAACCCTCTGCGGAAAAGGGCATCGCGTc                                                                                                >  2:279821/1‑150 (MQ=255)
                                                    atcaaGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCgg                                                                                              >  2:1260999/1‑150 (MQ=255)
                                                     tcaaGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCGTCGCGTCGGt                                                                                             <  1:426603/150‑1 (MQ=255)
                                                      caaGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTg                                                                                            <  2:608941/150‑1 (MQ=255)
                                                        aGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCg                                                                                          <  1:62077/150‑1 (MQ=255)
                                                           aGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAg                                                                                       >  2:626161/1‑150 (MQ=255)
                                                           aGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAg                                                                                       <  1:504260/150‑1 (MQ=255)
                                                           aGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTAAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAg                                                                                       >  2:273996/1‑150 (MQ=255)
                                                           aGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGCCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAg                                                                                       >  1:1139999/1‑150 (MQ=255)
                                                              cATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCaa                                                                                    >  1:1683428/1‑150 (MQ=255)
                                                              cATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCaa                                                                                    >  1:732446/1‑150 (MQ=255)
                                                                        ggttggttTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCaa                                                                          <  1:273996/150‑1 (MQ=255)
                                                                         gttggttTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCaag                                                                         <  1:1473408/150‑1 (MQ=255)
                                                                          ttggtgTGCACTCAGTACGAGCAGTGATAATCAGTACCGGTTTGTGCCAAAAGCTGGCTATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGGCGGAGCGTAGCAACAGCCCGCaaga                                                                        <  1:96925/150‑1 (MQ=255)
                                                                            ggttTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCaagaag                                                                      >  1:1523841/1‑150 (MQ=255)
                                                                             gttTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAagaaga                                                                     >  2:589965/1‑150 (MQ=255)
                                                                             gttTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAagaaga                                                                     >  1:17386/1‑150 (MQ=255)
                                                                                 gAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGc                                                                 >  1:548551/1‑150 (MQ=255)
                                                                                 gAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGc                                                                 >  2:796485/1‑150 (MQ=255)
                                                                                  aaCTCAGTAAGAGCAGTGATAATCAGTTCCGCTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGGCGGGGCGTAGCAACAGCCCGCAAGAAGAACGCt                                                                <  2:320615/150‑1 (MQ=255)
                                                                                    cTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTAt                                                              >  1:97470/1‑150 (MQ=255)
                                                                                     tCAGTAAGAGCAGTGATAATCAGTTCCTGTTTGTGCTAAAAGCGGGCAATTGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTTCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATg                                                             <  1:721146/150‑1 (MQ=255)
                                                                                      cAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATga                                                            <  2:1506610/150‑1 (MQ=255)
                                                                                      cAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATga                                                            <  1:611873/150‑1 (MQ=255)
                                                                                      cAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATga                                                            <  2:1643585/150‑1 (MQ=255)
                                                                                        gTAAGAGCAGTGATAATCAGATCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATgaga                                                          >  2:1634791/1‑150 (MQ=255)
                                                                                         tAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGaa                                                         <  1:1373185/150‑1 (MQ=255)
                                                                                         tAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGaa                                                         >  1:658008/1‑150 (MQ=255)
                                                                                           agagCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGaaaa                                                       >  2:415701/1‑150 (MQ=255)
                                                                                              gCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGaaaaaaa                                                    >  1:1489103/1‑150 (MQ=255)
                                                                                               cAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAAc                                                   >  2:946347/1‑150 (MQ=255)
                                                                                               cAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAAc                                                   <  1:963750/150‑1 (MQ=255)
                                                                                                 gTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTAACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACgg                                                 >  1:389377/1‑150 (MQ=255)
                                                                                                  tGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGc                                                >  2:629616/1‑150 (MQ=255)
                                                                                                     tAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCTCCAGCGAGCTTAATACCTCAAAACCCTCTGCGGAAAAGGGAAACGAGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAg                                             >  2:732294/1‑150 (MQ=255)
                                                                                                     tAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGCGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGACCGCTATGAGAAAAAAACGGCAAg                                             >  2:725918/1‑150 (MQ=255)
                                                                                                       aTCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTa                                           >  2:332229/1‑150 (MQ=255)
                                                                                                         cAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAAc                                         <  1:1514696/150‑1 (MQ=255)
                                                                                                          aGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCTGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCAATGAGAAAAAAACGTCAAGTaaag                                        >  2:289276/1‑148 (MQ=255)
                                                                                                          aGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACg                                        >  1:567184/1‑150 (MQ=255)
                                                                                                            ttCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGc                                      <  1:1460319/150‑1 (MQ=255)
                                                                                                            gtCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGc                                      <  1:1479657/149‑1 (MQ=255)
                                                                                                                ggTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAAt                                  <  1:626161/150‑1 (MQ=255)
                                                                                                                    tgtgCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAACCCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTGTGAGAAAAAAACGGCAAGTAACGGCAAATTCTa                              >  1:1111992/1‑150 (MQ=255)
                                                                                                                    tgtgCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTa                              >  1:29003/1‑150 (MQ=255)
                                                                                                                     gtgCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTAt                             >  1:885543/1‑150 (MQ=255)
                                                                                                                      tgCTAAAAGCGGGCAATGTTGAGACTATCCTCTCCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTAtt                            <  1:807633/150‑1 (MQ=255)
                                                                                                                       gCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTAttt                           >  1:1662035/1‑150 (MQ=255)
                                                                                                                        cTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCATGTAACGGCAAATTCTAttta                          >  2:1114758/1‑149 (MQ=255)
                                                                                                                        cTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTAttta                          >  2:1386206/1‑149 (MQ=255)
                                                                                                                        cTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCACTGCGCAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGCCAAGTAACGGCAAATACTAttta                          >  2:1164283/1‑149 (MQ=255)
                                                                                                                          aaaaGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCaa                        >  1:334220/1‑150 (MQ=255)
                                                                                                                          aaaaGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTGTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTACCGGCAAATTCTATTTCaa                        >  2:1272898/1‑150 (MQ=255)
                                                                                                                                gctcgctGGTGAGGATAGTCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACCGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTAGTTCAATCTGaa                  >  1:1162135/7‑150 (MQ=255)
                                                                                                                                 ctcgctGGTGAGGATAGTCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGaaa                 >  1:1229643/6‑150 (MQ=255)
                                                                                                                                 ctcgctGGTGAGGATAGTCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCAGaaa                 >  1:41132/6‑150 (MQ=255)
                                                                                                                                  gCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAg                <  2:1489103/150‑1 (MQ=255)
                                                                                                                                  gCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAg                <  1:415701/150‑1 (MQ=255)
                                                                                                                                  gCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAg                <  1:338168/150‑1 (MQ=255)
                                                                                                                                  gCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAg                >  2:1531333/1‑150 (MQ=255)
                                                                                                                                   cAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGc               >  1:845362/1‑150 (MQ=255)
                                                                                                                                    gctGGTGAGGATAGTCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGcc              >  1:909572/3‑150 (MQ=255)
                                                                                                                                    aaTGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACCGCCCGCAAGAAGAACGCTAAGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGcc              >  2:1172266/1‑150 (MQ=255)
                                                                                                                                     ctGGTGAGGATAGTCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGCCg             >  1:1011197/2‑150 (MQ=255)
                                                                                                                                            gaCTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGCCGCTAatca      <  2:1052251/150‑1 (MQ=255)
                                                                                                                                                aTCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAACACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGCCGCTAATCTTCaa  >  2:1604511/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGCCGCTAATCATCAA  >  NC_000913/2165123‑2165416

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