Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,338 T→A 11.1% intergenic (+140/‑74) pgm → / → ybfP phosphoglucomutase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,3380TA11.1% 185.5 / 3.0 76intergenic (+140/‑74)pgm/ybfPphosphoglucomutase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (37/28);  new base A (2/7);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 7.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGT  >  NC_000913/715100‑715566
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
tataAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGAt                                                                                                                                                                                                                                          >  2:784457/1‑239 (MQ=255)
tataAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGAt                                                                                                                                                                                                                                          >  2:770838/1‑239 (MQ=255)
tataAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGAt                                                                                                                                                                                                                                          >  2:469075/1‑239 (MQ=255)
        cTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGAATGACAAACACAGAATTGCCTCATGCGATACGCTTCTCAGGCCTACCAGGATGGGGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGCCgg                                                                                                                                                                                                                                  >  2:925557/1‑237 (MQ=255)
         tACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCATTATATTGAATTTAAGCGATTTTGTAGGCCGGTTAAGGCGTTc                                                                                                                                                                                                                                >  2:1077503/1‑239 (MQ=255)
         tACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATAAGGCCTAGGAGGATGGTGCAAAATATTGAATTTAAACGCTTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                                                >  2:587405/1‑239 (MQ=255)
             gCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                                                                                                            >  2:1203223/1‑238 (MQ=255)
             gCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTAAGAGGATGGGGCAATATATTGAATTTAAGCGATTTTGTAAGCCGGAAAAAGCGTTAAcg                                                                                                                                                                                                                            >  1:262909/1‑238 (MQ=255)
                      tCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTAATGCGCTACGCTTATCAGGCCTACGGGGATGGTGCAATATATTGAATTTAAGCGGTTTTTTAGGCCGGATAAGGGGTTACAGCCGCAACCgg                                                                                                                                                                                                                 >  1:873205/1‑239 (MQ=255)
                               aagaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAaaca                                                                                                                                                                                                         >  2:972791/1‑238 (MQ=255)
                               aagaaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCAACGAGGATGGGGCAAACTATTGAATTTAAGCGATTTTGAAGGCCGGAAAAGGGGTTAACGCCCCAACCGGCaaaaaaaa                                                                                                                                                                                                        >  1:202903/1‑239 (MQ=38)
                                  aaCATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGGTAAGGCGTTCCAGCCCCTTCCGGCAAAAACAACGa                                                                                                                                                                                                     >  2:32526/1‑239 (MQ=255)
                                       cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACt                                                                                                                                                                                                >  2:776751/1‑239 (MQ=255)
                                       cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACt                                                                                                                                                                                                >  2:766089/1‑239 (MQ=255)
                                       cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACt                                                                                                                                                                                                >  2:904495/1‑239 (MQ=255)
                                       cGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACAAACACt                                                                                                                                                                                                >  2:100318/1‑239 (MQ=255)
                                             caaaTTGAGAAAGAAGCGGTTGGAATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTca                                                                                                                                                                                          <  2:873205/236‑1 (MQ=255)
                                                 atGAAAAAAAAGCGGTTGAAATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTca                                                                                                                                                                                          <  2:262909/234‑1 (MQ=255)
                                                    agaaagaaGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTcaacaa                                                                                                                                                                                      >  2:528176/1‑236 (MQ=255)
                                                       atgaagCGGTTGGAATTTTTGGGTGAGTTCGGAAAACCGAGAAAACACATTTAATAAAAAAAGGGCGGTCGGAAGACCCCCCTTTTTTACGTATCACAAACAAAGAAGTGCTTGATGCGCTACGCTTATCAGGCATACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                                <  1:1077503/237‑1 (MQ=25)
                                                       aagaggttgTTTAGATTGTTAGCGAATTTCTGAAAACCGCGTAAACACATTTAAGAAAAAAAGGGCGTTCGCAAGATCGTCCTTGTTTACGTATGACGAACACAGAATTGCCTGTTGCGCTACACTTATCAGGCCTACGAGGATGGTGCAATATATGGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                                <  1:231351/231‑1 (MQ=37)
                                                       aagaagGGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACAAATTTAATAAAAAAAGGGCGGGCGCAAGATCGCCCTTTGTTACGTATGACAAACACAGAATTGCCTGAGGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                                <  2:721196/239‑1 (MQ=255)
                                                       aagaagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAg                                                                                                                                                                                <  2:202903/239‑1 (MQ=255)
                                                          aagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTa                                                                                                                                                                              >  2:1243461/1‑238 (MQ=255)
                                                          aagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTa                                                                                                                                                                              >  2:1272962/1‑238 (MQ=255)
                                                          aagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTGATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGAGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAg                                                                                                                                                                             <  2:1089835/239‑1 (MQ=255)
                                                          aagCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGAAGGGCTACCCTTCTCAGGCCTACGAGTGTGGTGCAAAATCTTGAATTTAAGCGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGCCAAAAAACTGAGAAg                                                                                                                                                                             >  1:231796/1‑239 (MQ=38)
                                                              ggTTGAAATTTTTAGGGAATTTATGAAAACCGCTTAAACCCATTTAAAAAAAAAAGGGCTGGCGCAAGATCGTCATTTTTTACGTATTACAAACACAAAATTGCCTTATGCGCTACGTTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTc                                                                                                                                                                          <  1:972791/238‑1 (MQ=25)
                                                                  ttgtttagtAGCGATGTTCTGAAAAACCGTGAAACACATTGAATAAAAAAAGGGCGGTCCCAATGTCGCCCTTTTTTACGTATGCCAAACACAGAATTGCGTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATAGTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGa                                                                                                                                                                     <  1:775116/231‑1 (MQ=37)
                                                                  tattttGTTAGGGAAGTTCTAAAAAACGGGTAAACAAATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGa                                                                                                                                                                     <  1:766089/235‑1 (MQ=255)
                                                                          tAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCaa                                                                                                                                                             >  2:433460/1‑239 (MQ=255)
                                                                            gCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTGTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATg                                                                                                                                                           <  2:1278592/239‑1 (MQ=255)
                                                                              gAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAAAACAAAATTGCCTGATGCGCTAGGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTCAGCGCTTTTGCAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAAACACAACGAACCCTTTTTCACAAAACTGAGAAGCTAAAGGAAAACCCaaaaaa                                                                                                                                                         >  1:655306/1‑235 (MQ=37)
                                                                               aaGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAAACCGAACACTTTGTCAACAAAATGAAAAGCCCAAGGAAAACCAAATaaaa                                                                                                                                                        >  1:79903/1‑238 (MQ=37)
                                                                                aGTTCTAAAAAACCCTTAAAACAATTTAAAAAAAAAAGGGCCGCTCCAAAGCGCCCGTTTTTCACTGATGAAAAACCCAGAATTGCCTGATGGTCTACGCTTATCAGGCCTACGGGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAg                                                                                                                                                        <  1:382494/238‑1 (MQ=21)
                                                                                   tCTGAAAAACGCGTAAACACATTTAAAAAAAAAAGGGCGGTCGAAAGATCGCCCTTTTTTACGTATGAAAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGAc                                                                                                                                                      <  2:952522/237‑1 (MQ=255)
                                                                                    cTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCaa                                                                                                                                                   <  2:1206318/239‑1 (MQ=255)
                                                                                               ttAAAAACTTTTAATAAAAAAAGGGGGGGTGCAAAGTCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTTGAGCGCTACGCTTAACAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTgg                                                                                                                                        <  1:1243461/238‑1 (MQ=37)
                                                                                                       caggtggttataaggagaaggTCGCAAGATCGCCCTTTTTTACGAAGGAAAAACACGAAATTGCCTGAGCCGCTACGCTTACCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCa                                                                                                                                  <  1:1256/219‑1 (MQ=37)
                                                                                                        tgtaaagaaaaaaaGGGCGGTCGAAAGATCGGCGTTTTTTACGTTGGACAAACACAGAATTGCTTGATGCGCTACGATTATCAGGCTTACGAGGATGGTGCAATATATGGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATa                                                                                                                                <  2:884862/232‑1 (MQ=37)
                                                                                                           taaaaaAAAAAGGGCGGTCCAAAGGTCCCCCTTTTTTACGTATTGCAATCACAGAATTGCCTGATGCGCTACGTTTACCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAgtg                                                                                                                             <  2:110913/239‑1 (MQ=37)
                                                                                                              aaaaaaagGGGGGGGCGCAATGTCGCCTTTTTTTACGTATGACAAACAAAGAATTGCATGATGCGTGACGCTTAGCAGGCCTACGAGGATGGTGCAATATATGGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCa                                                                                                                          <  2:370403/238‑1 (MQ=37)
                                                                                                                          ggtccaaatttCGCTCTTTTTTACGTAGACAAAAAAAAGAATTTCTTGTGGCTCTACGCTTACCAGCCCTGAGAGGAGGGTGAAATTTTTTGAATTTAAGAGTTTTTGGAGCCCGAATAAGGCGTCCACGCCGCACCCGGCAAAAACAACGAACACTTTGTAAACAAACTGAGTAGCTAAAGAAACTCCCAATGAAGCCAAACAGGAGTGTGGTCGTCATAGTGTCACTAATTACCGc                                                                                                               <  1:1166881/228‑1 (MQ=255)
                                                                                                                                      gCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCCCt                                                                                                                        >  2:1256/1‑217 (MQ=255)
                                                                                                                                              ttttCGTTTTACAAACCCAAAATTTCCTTTTGGGCTGCGGTTTTCACGTCTTAGGGGGTGGGGGAATTTATTTAAATTAACGGTTTTTGGAAGGCCGGTAAGGCTGTCACCGCCCCATCCGGCAAAAACAACAAACACTTTCTCAACAACCTGAGTGGCTAAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTTTCCCTAATTACCGCCACATTACTGCTGACTGCa                                                                                           <  1:1203223/239‑1 (MQ=255)
                                                                                                                                               ttACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAAACCAACAGGAGACTGGTCGTCCTAGGGTCACTAAATACCCCAACAATACTGCTGACTGCATg                                                                                         >  1:1222227/1‑239 (MQ=37)
                                                                                                                                                       gACAACAACAGAATTGCCTGAGGCGCTACGCTTAGCAGGCCTACGAGGATGGTGCAATTTATTGAATTTAAGCGATTTGGTAGGCCGAAAAAGGCTTTCACGCGGCATCCGGCAAAAACAACGAACACTTTGACAACAAACTGATTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAAc                                                                                 <  2:894352/239‑1 (MQ=255)
                                                                                                                                                             acacAGAATTGCTTGATGCGCTACGTTTTACAGGCCAACGAGGAGGGTGCAATATTTGAAATTAAAGGGTTTTTGTGGGCCGAATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACCTTTTGTCAACAACCTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGaa                                                                            <  2:79903/239‑1 (MQ=255)
                                                                                                                                                                        ccTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGATCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCACTGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACAGCCACATTACTGCTGACTCCATGCGCCCAACCTGAACAATCATCAtt                                                                 >  2:881699/1‑239 (MQ=255)
                                                                                                                                                                        ccTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGACAAGGCGTTCACGCCGCATCCGGCACAAACAACGAACCCTTTGTCAACAAACTGAGCAGCCCAAGGAAACACCAATGAAACCCACCAGCAGCATGGCCGCCATCGGGTCAGTAATTACCCCCACATTAATGCTGACTGCATCCGCCCAACCTGAACACTCATAAtt                                                                 >  1:846413/1‑239 (MQ=25)
                                                                                                                                                                             tGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGGGTCACTAATTACCGCCACATTActgct                                                                                                  >  2:504515/1‑201 (MQ=255)
                                                                                                                                                                               cgcTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                                                                                                            <  2:813322/79‑1 (MQ=255)
                                                                                                                                                                               cgcTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                                                                                                            >  1:813322/1‑79 (MQ=255)
                                                                                                                                                                                         aTCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGGGTAGCTCAAGGAAATCCCAATGCAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGGCGCGCTGATGCCCCCAACTTGAACAAtcctcacaacaagcagtcgggtaaa                                                >  1:468955/1‑220 (MQ=255)
                                                                                                                                                                                            tggataaagggaaggttgaaaTAAATTGATTTTAAGGATTTTTGAGGGCGGGGATAGGCGTCCGCCACGCTACCGGAAAAAAAAAACAAAAATTTTGCAACAAACTGAGTAGATAAAGGAACTCCCAATTAAGACCAACATGAGTGTGGCCGTCATGGTGTCCATAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGa                                             <  1:107924/220‑1 (MQ=255)
                                                                                                                                                                                              gCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTActg                                                                                                    >  2:925020/1‑182 (MQ=255)
                                                                                                                                                                                                 tacaagaaagttttatcatatTGAGTTTAGAGTATTTGGTGGGCGGGATAGGGGGTAGACTCCGAGTCCGGAAAAAAAAACGAACCCTTTGACAACAACCTGGGTAGTCTAGGGAACTCCCAATGAAAAACCAGAGGGGTTTGGACGTCAACGTGTCACAAATTGCCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGTCCCCt                                        <  1:466951/222‑1 (MQ=255)
                                                                                                                                                                                                 tACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAac                                                                                                                                                                                                   <  1:1067734/85‑1 (MQ=255)
                                                                                                                                                                                                 tACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAac                                                                                                                                                                                                   >  2:1067734/1‑85 (MQ=255)
                                                                                                                                                                                                   cGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAATGTCACTAATTACCCCCACATTACTGAGGACTGCATGCGCCCAACCCGAAAAATAATCATTACCTCGCCACTGGTTATTGACCCCTaa                                      >  1:146080/1‑239 (MQ=38)
                                                                                                                                                                                                         tGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTTTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGCCCCCTAAAGATaa                                >  2:959282/1‑239 (MQ=255)
                                                                                                                                                                                                                 ttattggAATTTAAGCATTTTTGGAGCCGGGAAAAGGGTTCCACGCCGCACCCGGCAAAAACAACGAACCTTTTGTCAACAAACTGAGTAGCTAAAGGAAATCCCAATGAAGCCCAACAGGAGTCTGGTCGTCATAGTGTCATTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCagag                        <  2:231796/233‑1 (MQ=255)
                                                                                                                                                                                                                      tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGCGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCACCATTAGCTGGCGACTGGTTACTGACCCCTAAAGAtctgcctcttataacc                   >  2:107924/1‑224 (MQ=255)
                                                                                                                                                                                                                      tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAGGGAAATCCAAATGAAACTCAACAGAAGTCTGGCCGTCATAATGTCACTAATTACCCCAACATTACGGATAACTGCATGCGCCCAACCTGAACAATCATCAATAGCTAGCGAATGGCTAATGACCCCAAAAGaaaaaaaaaaagaataa                   >  1:631881/1‑228 (MQ=255)
                                                                                                                                                                                                                      tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTCCTGACCCCTAAAGATAAAACCAGAGGATTa                   >  2:246146/1‑239 (MQ=255)
                                                                                                                                                                                                                      tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa                   >  2:1164908/1‑239 (MQ=255)
                                                                                                                                                                                                                      tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTAACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATaaca                              >  2:466951/1‑226 (MQ=255)
                                                                                                                                                                                                                      tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGGCTGGTCGGCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGAATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa                   >  2:43607/1‑239 (MQ=255)
                                                                                                                                                                                                                      tGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGCCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTa                   >  2:32403/1‑239 (MQ=255)
                                                                                                                                                                                                                        aaTTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACACGAGTCTGGTCGTCCTAGAGTCACTAATTACCGCCACATTCCTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGGCCCCTAAAGAAAAAAcaaagggagaaac                 >  1:881067/1‑227 (MQ=255)
                                                                                                                                                                                                                           taaaGCGTTTTTGGAGGCCGGAAAAGGCTTCCACGCCGCAACCGGAAAAAACAACGAACCCTTTGTCAACAAACTGAGTAGCTAAAGGAAACCCCAATGAAGCCAAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACggg              <  2:1222227/237‑1 (MQ=255)
                                                                                                                                                                                                                                 gttttgGTAGGCCGTATAGGGCTGCCACGCCGCATCCGGAAAAAACAACGAACCTTTTGTCAACAAACGGAGTAGCTAAAGGAAATCCCAATGAAGCCCAACAGGATGCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCa         <  2:354621/236‑1 (MQ=255)
                                                                                                                                                                                                                                 gggtaaaagagacaggAAAAGGCGTCCCCGCCGCACCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAACCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCa         <  2:48964/225‑1 (MQ=255)
                                                                                                                                                                                                                                 gATTTGGTAGGCCGAATAAGGCGTTCCCGCCGCATCCGGAAAAAACAACGAACCTTTTGTAAACAAACGGAGTAGCTCAAGGAAATCCCAAGGAAGACCAACAGGAGTCTGGTCGTAATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCa         <  2:591775/238‑1 (MQ=255)
                                                                                                                                                                                                                                  aattgttaaGGCTGGAAAAGGGTTTCCCGCCCCATGCGGACAAAACAACGAAATTTTTGACAAAAAACGGAGTAGCTCAGGGAAATCCCAATGACCACCAACAGGAGTCTGGTCGTTAAAGTGTCACTAATTCCCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCa         <  2:727406/230‑1 (MQ=255)
                                                                                                                                                                                                                                  aTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAAAGGATTACCGGGTAGCAtt       >  2:1012880/1‑239 (MQ=255)
                                                                                                                                                                                                                                   ttttttgcggtcGGATAGGGGGTTCCGCGCGCAACGGGAAAAACAAACGAACAATGTTGCAACAAACGGGATAGCTAGAAGAAACCCAAGGGAAAACCAGAGGGAGTCTGGTCGTCATAGGATCACAAATAACCGCCACATTACTGCTGACTGCATGCGCCCAACATGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCAtt       <  1:194085/228‑1 (MQ=255)
                                                                                                                                                                                                                                        gTAGGCCGGAAAGGGCTTTCACCCCGACTCCGGCAAAAACAACGAACCCTTTGTCAACAAACTGAGTAGCACAAGGAAACCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTTTCCCTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGt  <  2:1102400/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gTAGGCCGGAAAAGGTCTCCACGCCGCATCCGGAAAAAAAAACGAACCTTTTGTAAACAAGCGAAGTAGCTAAAGGAAATCCCAATGAAGACCAACAGGAGTTTGGGCTTCATAGTTTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGt  <  2:655306/239‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
TATAAGATCTACTGCGAAAGCTTCCTCGGTGAAGAACATCGCAAGCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAACACTTTGTCAACAAACTGAGTAGCTCAAGGAAATCCCAATGAAGACCAACAGGAGTCTGGTCGTCATAGTGTCACTAATTACCGCCACATTACTGCTGACTGCATGCGCCCAACCTGAACAATCATCATTAGCTGGCGACTGGTTACTGACCCCTAAAGATAAAACCAGAGGATTAACGGGTAGCATTGCGGT  >  NC_000913/715100‑715566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: