Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,290,375 A→G 13.4% E45G (GAG→GGG)  rssB → PcnB‑degradosome interaction factor; response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,290,3750AG13.4% 186.5 / 2.8 78E45G (GAG→GGG) rssBPcnB‑degradosome interaction factor; response regulator
Reads supporting (aligned to +/- strand):  ref base A (24/41);  new base G (9/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 7.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACATTTGACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCT  >  NC_000913/1290142‑1290608
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
aaCATTTGACCAGAATTTTTATCTACCCTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGc                                                                                                                                                                                                                                        <  1:171140/239‑1 (MQ=255)
       gACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGGTTCCCTTGGGTTGGTgggggg                                                                                                                                                                                                                                 >  2:185147/1‑239 (MQ=255)
       gACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGAATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACTGTACTGGCGGCTGCTGGGGTGGGTGCCCTTTGGTTGGTGTGAgg                                                                                                                                                                                                                                 >  1:951998/1‑239 (MQ=255)
          caaaaTTTTTTTCTTCCATTAAGTTAATTCGGAAAGGGGGAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAAGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGtt                                                                                                                                                                                                                               <  2:820729/235‑1 (MQ=255)
          cAGAATTTTTATCTCCACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGttt                                                                                                                                                                                                                              <  1:766696/239‑1 (MQ=255)
            gAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGAAGCCCTTTGGTTGCTGGGAGGCTTc                                                                                                                                                                                                                             >  1:244841/1‑238 (MQ=255)
            gAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGAATTTCGCTCGCTTCTGGATTAATGGTTTTCCTCATTGGAGGCGCCAACGGTACTGGCGGCTGATTGGGGGGGTGCCCTTGGGTTGCAggagggatta                                                                                                                                                                                                                             >  2:215767/1‑230 (MQ=255)
             aaTTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGATGCCCTTGAGTTGCTGGGAGGTTTCa                                                                                                                                                                                                                            >  1:283202/1‑238 (MQ=255)
                 tttATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCa                                                                                                                                                                                                                       <  1:18129/239‑1 (MQ=255)
                    atttACCCTTAAGTTAATTTGGACCGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGAc                                                                                                                                                                                                                    <  2:684709/236‑1 (MQ=255)
                    aTCTACACTTAAGTTAATTCGGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGAc                                                                                                                                                                                                                    <  1:96093/239‑1 (MQ=255)
                    aTCTACACTTAAGCTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCCGGGGAGAGAACATGACGCAGCCACTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTAGCTTCTGGCTTCATGGCTTTCCTCATTGGGAGCGACAACGGTCCTGGCGGCTCAACGGGTGGATGCCCTTTAGTTGTCGGCAGGCTTCCCTCCCGa                                                                                                                                                                                                                     >  1:202650/1‑238 (MQ=255)
                       tACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGGTGCCCTTTGGTTGGTGGGAGGGTTTACCCCCGACAtg                                                                                                                                                                                                                 >  2:986216/1‑239 (MQ=255)
                         cacTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGTTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCtga                                                                                                                                                                                                                <  1:687764/238‑1 (MQ=255)
                         cacTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGAGTGGGTGGAGTCCCTTGAGTTGCTGGGGGGTTTCACTCCAGACCtgat                                                                                                                                                                                                               >  1:684709/1‑239 (MQ=255)
                                   aaTTCTGACATGAGCAGGCGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGTAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCACTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGGTGATGGGGTGGATGCCCTTCAGATGCCGGGAGGTTTCACTCCAGCCCTGATGGTT‑TGGGat                                                                                                                                                                                                    >  2:687764/1‑239 (MQ=255)
                                   aaTTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGATGCCCTTTGGTTGGTGGGGGGTTTTACCCCAACCCTGATGATA‑tgtggg                                                                                                                                                                                                    >  2:1035337/1‑237 (MQ=255)
                                               ggcAGGGGGCAATGGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACCTGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGcc                                                                                                                                                                                        <  1:1058437/237‑1 (MQ=255)
                                               cgcAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGGTAT‑CGCGATGcc                                                                                                                                                                                        >  1:580977/1‑238 (MQ=255)
                                                    gTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGCTGGGGGGGATGCCCTTAGGTTGCTGGGAGGTTTCCCTCCGGACCTGATGATA‑TGTGAGAT‑CGCGATCCCAGGAAt                                                                                                                                                                                  >  1:370732/1‑239 (MQ=255)
                                                    gTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAAt                                                                                                                                                                                  >  1:897885/1‑239 (MQ=255)
                                                        cAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGaa                                                                                                                                                                               <  2:897885/238‑1 (MQ=255)
                                                        cAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGaa                                                                                                                                                                               <  1:797174/238‑1 (MQ=255)
                                                        cAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGaa                                                                                                                                                                               <  1:665341/238‑1 (MQ=255)
                                                        cAAGAGCATGGCACTATTGAATAAAGCCAAGCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGaa                                                                                                                                                                               <  1:136977/238‑1 (MQ=255)
                                                             gCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCt                                                                                                                                                                         <  1:173589/239‑1 (MQ=255)
                                                               aTGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATAGGGTGGATGCCCTTGGGTTGCTGGGGGGGTTCCCTCCAGACCTGTTGTAT‑TGTG‑TATACGCGGTGCCACGAAAGAACGGGGTTa                                                                                                                                                                       >  2:867321/1‑239 (MQ=255)
                                                                 gCCAATATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTaaa                                                                                                                                                                     <  1:756007/239‑1 (MQ=255)
                                                                     cTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAActgc                                                                                                                                                                 <  1:1035337/239‑1 (MQ=255)
                                                                      tATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAActgct                                                                                                                                                                <  1:108175/239‑1 (MQ=255)
                                                                       aTTGAGTAAAGGCAGTTAGGGGAGAGAACGTGGCGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAActgctg                                                                                                                                                               <  1:58811/239‑1 (MQ=255)
                                                                           aGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTTGAGc                                                                                                                                                           >  1:172339/1‑239 (MQ=255)
                                                                            gTAAAGCCAGTAAGGGGAGAGAACAAGGCGCAGCCATTGGTCTGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGc                                                                                                                                                           <  1:367133/238‑1 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGGCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGttt                                                                                                                                                                                                                              <  2:162326/172‑1 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTTAGTTGCTGGGAGGTTTTACTCCAGACCTGATGATATTGTG‑TAT‑CGGGATGGCCCGACTGAACGGGCTTAAACTGCTGGGGCat                                                                                                                                                         >  1:368721/1‑239 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGttt                                                                                                                                                                                                                              >  1:162326/1‑172 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTCAACTGCTGGGGCat                                                                                                                                                         >  1:887941/1‑239 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTACCTCCAAGCCTGATGATA‑TCGGGTAT‑TCCGATGCCCCGCATGAACGGGCTTAAACTGATGGAGcct                                                                                                                                                         >  1:612712/1‑237 (MQ=255)
                                                                              aaaGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACATATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCat                                                                                                                                                         <  1:986216/238‑1 (MQ=255)
                                                                              aaaGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCat                                                                                                                                                         <  1:25581/238‑1 (MQ=255)
                                                                              aaaGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTACGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCat                                                                                                                                                         <  1:300777/238‑1 (MQ=255)
                                                                                aGCCCGTCAGGGGAGAGAACTTGTCGCAGCTATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCAtata                                                                                                                                                      <  2:484774/239‑1 (MQ=255)
                                                                                aGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGGTGCCCTTGAGTTGCTGGGGGGGTTCACTCCAGACCTGATGTTA‑TGTGATAT‑CGCGATGCCACGACTGAATGGGCATAAAATGGTGGGGCaaaaa                                                                                                                                                      >  2:543705/1‑235 (MQ=255)
                                                                                aGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTTGGAGGGTTCACTCCAGACCTTGTGGAA‑TTGGAGAT‑TGCGGTGCCCCGAATGAAAGGGCTTAAACTGCTgggggaaaaa                                                                                                                                                      >  2:108175/1‑231 (MQ=255)
                                                                                aGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGCGGTTTCACTCCAGACCTGGTGAGA‑TGTGATAT‑CCCGATTCCCCCAATGAAAGGGCTTAAACAACTGGGGCataaa                                                                                                                                                      >  2:58811/1‑237 (MQ=255)
                                                                                     gTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGTGTTGCTGGGAGGGTTTCCTCCCGGCCTGGTGGTA‑TGTGATAT‑GGTGGTGCCTGGCAGTAAAGGGCTTAAAACGCTGGGGCATAACCGTaa                                                                                                                                                 >  2:620359/1‑239 (MQ=255)
                                                                                          gggAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCCCTCCAGACCTGATGATA‑TGTGCTAT‑CGCGATGCCCCGAATGAACGGGCTTAAACTGCTGGGGCATAAACGGAACagag                                                                                                                                            >  2:796696/1‑239 (MQ=255)
                                                                                                    aTGACGCAGCCATTGGTCGGCAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGATGCCCTTGGGATGCGGGGAGGGTTCACTCCCAACCTGGTGAGC‑GGTGATAT‑CGCGATGCCAGGGAGTAAAGGGCTTAAACCGCCGGAGCCAAAAAGTAAAAGAAGCGACAAGa                                                                                                                                   >  2:564132/1‑238 (MQ=255)
                                                                                                      gACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGCTTGGGTGGATTCCCTTTAGTTGCGGGGAGGTTTCACTCCCGACCTGATGATA‑TGTTATAT‑TCGGGTGCCACGAATGAACGGGATTAAACTTCCGGGGAATAAACCGAACAGAGGCggaccaaac                                                                                                                                 >  2:716749/1‑230 (MQ=255)
                                                                                                           aGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGGGGGTTCACTCCAGACCTGCTGATA‑TGTGATAT‑CGGGATGCCACGAATGAACGGGCTTAAACTGCTGGGGCCTATCAGAAACAAAGGCGGCCAAACCCCAgg                                                                                                                            >  2:817466/1‑237 (MQ=255)
                                                                                                           aGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑TGCGATGCCACGAATGAACGGGCTTAAACTGCTGGGTCAAATACGGAAAAAGAGCGGCCAAACCCCaagt                                                                                                                           >  1:462688/1‑236 (MQ=255)
                                                                                                                       gAACACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGGGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCGCTTGCGTTGTTGGGAGGTTTCGCTCCCAGCCTCATGGTA‑TGTGATAT‑CGCGATGCCACCACTGAACAGGCTTCAAATGCTGGCGCATATACGTAACAAAGGCGACCCACCCCCCGTTCTgtgggac                                                                                                                  >  2:850055/1‑230 (MQ=255)
                                                                                                                             aGATTCCCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACt                                                                                                         <  1:817466/239‑1 (MQ=255)
                                                                                                                                            aaGATGGGCAGGTATTTCGTTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCCTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGat                                                                                          <  1:1059735/239‑1 (MQ=255)
                                                                                                                                              gATGGGGAGGTTTTTCGCTCGCTTCTGTATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGata                                                                                         <  2:580977/238‑1 (MQ=255)
                                                                                                                                                   gCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGTTGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAAAGGCGACCAGAc                                                                                                                                  >  2:702765/1‑192 (MQ=255)
                                                                                                                                                   gCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGAc                                                                                                                                  <  1:702765/192‑1 (MQ=255)
                                                                                                                                                    cAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCa                                                                                   <  1:796696/238‑1 (MQ=255)
                                                                                                                                                     cggTTTTTCGCTCGCTTTTGGGTTCATGGTTTTCCCTGTTGGGGGGGACAACGGTACGGGCGGGTGATGGGGTGGATGCCCTGGAGTTGCTGGGAGGTTTCACTCCAGACCTGAGGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCaa                                                                                  <  2:612712/237‑1 (MQ=255)
                                                                                                                                                     aGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCaa                                                                                  <  1:185147/238‑1 (MQ=255)
                                                                                                                                                         aTTTCGCTCGCTTCTGGATTCAGGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGt                                                                             <  1:602759/239‑1 (MQ=255)
                                                                                                                                                                 cgctTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTggg                                                                                                                                                                                                                                    <  2:379980/82‑1 (MQ=255)
                                                                                                                                                                 cgctTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTggg                                                                                                                                                                                                                                    >  1:379980/1‑82 (MQ=255)
                                                                                                                                                                              aTGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTgg                                                                    <  2:94420/227‑1 (MQ=255)
                                                                                                                                                                              aTGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTgg                                                                    >  1:94420/1‑227 (MQ=255)
                                                                                                                                                                              aTGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCCGGAGCATAAACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGGAAATATGGCAGAAATTGACAAACCGTTACGTCTGGCGTTTGatgatg                                                        >  1:911968/1‑239 (MQ=255)
                                                                                                                                                                                gtttttCTTAATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGtt                                                      <  2:627647/237‑1 (MQ=255)
                                                                                                                                                                                  ttttCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACGCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGaa                                                                                                                                                                               <  1:988007/116‑1 (MQ=255)
                                                                                                                                                                                   tttCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCtga                                                                                                                                                                                                                >  1:725599/1‑84 (MQ=255)
                                                                                                                                                                                   tttCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCtga                                                                                                                                                                                                                <  2:725599/84‑1 (MQ=255)
                                                                                                                                                                                           ttGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAAcc                                           <  2:172339/239‑1 (MQ=255)
                                                                                                                                                                                                    aCAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGGAATCTGCCCCTGAAAATATGGCAGAAAATGCCCAAGCGTTACGTTTGGGCGTTGGAAAGTTTTTGTTGAAAACAGTTaaaaa                                  >  2:798919/1‑237 (MQ=255)
                                                                                                                                                                                                                  ggcTGTTGGGGTGGATTGCCTTTAGTTGCTGGGAGGTTTCCCTCCAGGACTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCtgg                     <  2:368721/238‑3 (MQ=255)
                                                                                                                                                                                                                  ggcTGATGGGGTGTATGCCCTTGGGTTGTGGGGAGGTTTCCCTCCAGACCTGATTATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGGGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTg                     <  2:913933/238‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGGTTTGCTGGGAGGTTTCCCTCCAGGCCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCt  <  2:887941/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGAGTTGCTGGGGGGTTTCACTCCAGGCGTTATGATA‑AGTGATGT‑CGCGGTGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCATCCAGAACCAAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCt  <  2:244841/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCt  <  1:639464/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGGGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCgaga              <  1:1064342/227‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
AACATTTGACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATA‑TGTGATAT‑CGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCT  >  NC_000913/1290142‑1290608

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