Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,264 A→G 8.1% T26A (ACT→GCT) ‡ yegP → UPF0339 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,2640AG 8.1% 281.8 / 12.3 99T26A (ACT→GCT) ‡yegPUPF0339 family protein
Reads supporting (aligned to +/- strand):  ref base A (54/37);  new base G (5/3);  total (59/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 = 8.37e-01

CGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGCCGCTAATCA  >  NC_000913/2165116‑2165412
                                                                                                                                                    |                                                                                                                                                    
cGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGAc                                                                                                                                                     <  1:1498338/150‑1 (MQ=255)
       aCTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTc                                                                                                                                              <  2:387809/150‑1 (MQ=255)
       aCTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTc                                                                                                                                              >  2:854463/1‑150 (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:537509/150‑4 (MQ=255)
              tCGGAAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGGATAGTCTCACCattg                                                                                                                                       <  2:28024/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:1673623/150‑1 (MQ=255)
                         cACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATAcc                                                                                                                            >  2:1308937/1‑150 (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                                                                                                    >  1:1608378/1‑150 (MQ=255)
                                                 cacaATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGc                                                                                                    >  2:1533086/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                                                                                   <  1:504260/150‑1 (MQ=255)
                                                                  aGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAg                                                                                   >  2:626161/1‑150 (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                                                                 >  1:17386/1‑150 (MQ=255)
                                                                                    gttTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAagaaga                                                                 >  2:589965/1‑150 (MQ=255)
                                                                                        gAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGc                                                             >  2:796485/1‑150 (MQ=255)
                                                                                        gAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGc                                                             >  1:548551/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                                                        <  1:611873/150‑1 (MQ=255)
                                                                                             cAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATga                                                        <  2:1506610/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:1531333/1‑150 (MQ=255)
                                                                                                                                         gCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAg            <  1:338168/150‑1 (MQ=255)
                                                                                                                                         gCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAg            <  1:415701/150‑1 (MQ=255)
                                                                                                                                         gCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAg            <  2:1489103/150‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
CGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAAAACCTCTGCGGAAAAGGGCATCGCGTCGGTGCGTAGCAACAGCCCGCAAGAAGAACGCTATGAGAAAAAAACGGCAAGTAACGGCAAATTCTATTTCAATCTGAAAGCCGCTAATCA  >  NC_000913/2165116‑2165412

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: ‑