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

AACATTTGACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCT  >  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)
                                   aaTTCTGACATGAGCAGGCGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGTAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCACTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGGTGATGGGGTGGATGCCCTTCAGATGCCGGGAGGTTTCACTCCAGCCCTGATGGTTTGGGat                                                                                                                                                                                                    >  2:687764/1‑239 (MQ=255)
                                   aaTTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGATGCCCTTTGGTTGGTGGGGGGTTTTACCCCAACCCTGATGATAtgtggg                                                                                                                                                                                                    >  2:1035337/1‑237 (MQ=255)
                                               ggcAGGGGGCAATGGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACCTGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGcc                                                                                                                                                                                        <  1:1058437/237‑1 (MQ=255)
                                               cgcAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGGTATCGCGATGcc                                                                                                                                                                                        >  1:580977/1‑238 (MQ=255)
                                                    gTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGCTGGGGGGGATGCCCTTAGGTTGCTGGGAGGTTTCCCTCCGGACCTGATGATATGTGAGATCGCGATCCCAGGAAt                                                                                                                                                                                  >  1:370732/1‑239 (MQ=255)
                                                    gTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAAt                                                                                                                                                                                  >  1:897885/1‑239 (MQ=255)
                                                        cAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGaa                                                                                                                                                                               <  2:897885/238‑1 (MQ=255)
                                                        cAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGaa                                                                                                                                                                               <  1:797174/238‑1 (MQ=255)
                                                        cAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGaa                                                                                                                                                                               <  1:665341/238‑1 (MQ=255)
                                                        cAAGAGCATGGCACTATTGAATAAAGCCAAGCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGGAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGaa                                                                                                                                                                               <  1:136977/238‑1 (MQ=255)
                                                             gCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCt                                                                                                                                                                         <  1:173589/239‑1 (MQ=255)
                                                               aTGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATAGGGTGGATGCCCTTGGGTTGCTGGGGGGGTTCCCTCCAGACCTGTTGTATTGTGTATACGCGGTGCCACGAAAGAACGGGGTTa                                                                                                                                                                       >  2:867321/1‑239 (MQ=255)
                                                                 gCCAATATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTaaa                                                                                                                                                                     <  1:756007/239‑1 (MQ=255)
                                                                     cTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAActgc                                                                                                                                                                 <  1:1035337/239‑1 (MQ=255)
                                                                      tATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAActgct                                                                                                                                                                <  1:108175/239‑1 (MQ=255)
                                                                       aTTGAGTAAAGGCAGTTAGGGGAGAGAACGTGGCGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAActgctg                                                                                                                                                               <  1:58811/239‑1 (MQ=255)
                                                                           aGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTTGAGc                                                                                                                                                           >  1:172339/1‑239 (MQ=255)
                                                                            gTAAAGCCAGTAAGGGGAGAGAACAAGGCGCAGCCATTGGTCTGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGc                                                                                                                                                           <  1:367133/238‑1 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGGCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGttt                                                                                                                                                                                                                              <  2:162326/172‑1 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTTAGTTGCTGGGAGGTTTTACTCCAGACCTGATGATATTGTGTATCGGGATGGCCCGACTGAACGGGCTTAAACTGCTGGGGCat                                                                                                                                                         >  1:368721/1‑239 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGttt                                                                                                                                                                                                                              >  1:162326/1‑172 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTCAACTGCTGGGGCat                                                                                                                                                         >  1:887941/1‑239 (MQ=255)
                                                                             tAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTACCTCCAAGCCTGATGATATCGGGTATTCCGATGCCCCGCATGAACGGGCTTAAACTGATGGAGcct                                                                                                                                                         >  1:612712/1‑237 (MQ=255)
                                                                              aaaGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACATATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCat                                                                                                                                                         <  1:986216/238‑1 (MQ=255)
                                                                              aaaGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCat                                                                                                                                                         <  1:25581/238‑1 (MQ=255)
                                                                              aaaGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTACGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCat                                                                                                                                                         <  1:300777/238‑1 (MQ=255)
                                                                                aGCCCGTCAGGGGAGAGAACTTGTCGCAGCTATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCAtata                                                                                                                                                      <  2:484774/239‑1 (MQ=255)
                                                                                aGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGGTGCCCTTGAGTTGCTGGGGGGGTTCACTCCAGACCTGATGTTATGTGATATCGCGATGCCACGACTGAATGGGCATAAAATGGTGGGGCaaaaa                                                                                                                                                      >  2:543705/1‑235 (MQ=255)
                                                                                aGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTTGGAGGGTTCACTCCAGACCTTGTGGAATTGGAGATTGCGGTGCCCCGAATGAAAGGGCTTAAACTGCTgggggaaaaa                                                                                                                                                      >  2:108175/1‑231 (MQ=255)
                                                                                aGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGCGGTTTCACTCCAGACCTGGTGAGATGTGATATCCCGATTCCCCCAATGAAAGGGCTTAAACAACTGGGGCataaa                                                                                                                                                      >  2:58811/1‑237 (MQ=255)
                                                                                     gTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGTGTTGCTGGGAGGGTTTCCTCCCGGCCTGGTGGTATGTGATATGGTGGTGCCTGGCAGTAAAGGGCTTAAAACGCTGGGGCATAACCGTaa                                                                                                                                                 >  2:620359/1‑239 (MQ=255)
                                                                                          gggAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCCCTCCAGACCTGATGATATGTGCTATCGCGATGCCCCGAATGAACGGGCTTAAACTGCTGGGGCATAAACGGAACagag                                                                                                                                            >  2:796696/1‑239 (MQ=255)
                                                                                                    aTGACGCAGCCATTGGTCGGCAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGGGGATGCCCTTGGGATGCGGGGAGGGTTCACTCCCAACCTGGTGAGCGGTGATATCGCGATGCCAGGGAGTAAAGGGCTTAAACCGCCGGAGCCAAAAAGTAAAAGAAGCGACAAGa                                                                                                                                   >  2:564132/1‑238 (MQ=255)
                                                                                                      gACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGCTTGGGTGGATTCCCTTTAGTTGCGGGGAGGTTTCACTCCCGACCTGATGATATGTTATATTCGGGTGCCACGAATGAACGGGATTAAACTTCCGGGGAATAAACCGAACAGAGGCggaccaaac                                                                                                                                 >  2:716749/1‑230 (MQ=255)
                                                                                                           aGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGGGGGTTCACTCCAGACCTGCTGATATGTGATATCGGGATGCCACGAATGAACGGGCTTAAACTGCTGGGGCCTATCAGAAACAAAGGCGGCCAAACCCCAgg                                                                                                                            >  2:817466/1‑237 (MQ=255)
                                                                                                           aGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATTGCGATGCCACGAATGAACGGGCTTAAACTGCTGGGTCAAATACGGAAAAAGAGCGGCCAAACCCCaagt                                                                                                                           >  1:462688/1‑236 (MQ=255)
                                                                                                                       gAACACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGGGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCGCTTGCGTTGTTGGGAGGTTTCGCTCCCAGCCTCATGGTATGTGATATCGCGATGCCACCACTGAACAGGCTTCAAATGCTGGCGCATATACGTAACAAAGGCGACCCACCCCCCGTTCTgtgggac                                                                                                                  >  2:850055/1‑230 (MQ=255)
                                                                                                                             aGATTCCCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACt                                                                                                         <  1:817466/239‑1 (MQ=255)
                                                                                                                                            aaGATGGGCAGGTATTTCGTTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCCTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGat                                                                                          <  1:1059735/239‑1 (MQ=255)
                                                                                                                                              gATGGGGAGGTTTTTCGCTCGCTTCTGTATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGata                                                                                         <  2:580977/238‑1 (MQ=255)
                                                                                                                                                   gCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGTTGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAAAGGCGACCAGAc                                                                                                                                  >  2:702765/1‑192 (MQ=255)
                                                                                                                                                   gCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGAc                                                                                                                                  <  1:702765/192‑1 (MQ=255)
                                                                                                                                                    cAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCa                                                                                   <  1:796696/238‑1 (MQ=255)
                                                                                                                                                     cggTTTTTCGCTCGCTTTTGGGTTCATGGTTTTCCCTGTTGGGGGGGACAACGGTACGGGCGGGTGATGGGGTGGATGCCCTGGAGTTGCTGGGAGGTTTCACTCCAGACCTGAGGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCaa                                                                                  <  2:612712/237‑1 (MQ=255)
                                                                                                                                                     aGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCaa                                                                                  <  1:185147/238‑1 (MQ=255)
                                                                                                                                                         aTTTCGCTCGCTTCTGGATTCAGGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGt                                                                             <  1:602759/239‑1 (MQ=255)
                                                                                                                                                                 cgctTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTggg                                                                                                                                                                                                                                    <  2:379980/82‑1 (MQ=255)
                                                                                                                                                                 cgctTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTggg                                                                                                                                                                                                                                    >  1:379980/1‑82 (MQ=255)
                                                                                                                                                                              aTGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTgg                                                                    <  2:94420/227‑1 (MQ=255)
                                                                                                                                                                              aTGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTgg                                                                    >  1:94420/1‑227 (MQ=255)
                                                                                                                                                                              aTGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCCGGAGCATAAACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGGAAATATGGCAGAAATTGACAAACCGTTACGTCTGGCGTTTGatgatg                                                        >  1:911968/1‑239 (MQ=255)
                                                                                                                                                                                gtttttCTTAATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGtt                                                      <  2:627647/237‑1 (MQ=255)
                                                                                                                                                                                  ttttCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACGCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGaa                                                                                                                                                                               <  1:988007/116‑1 (MQ=255)
                                                                                                                                                                                   tttCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCtga                                                                                                                                                                                                                >  1:725599/1‑84 (MQ=255)
                                                                                                                                                                                   tttCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCtga                                                                                                                                                                                                                <  2:725599/84‑1 (MQ=255)
                                                                                                                                                                                           ttGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAAcc                                           <  2:172339/239‑1 (MQ=255)
                                                                                                                                                                                                    aCAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGGAATCTGCCCCTGAAAATATGGCAGAAAATGCCCAAGCGTTACGTTTGGGCGTTGGAAAGTTTTTGTTGAAAACAGTTaaaaa                                  >  2:798919/1‑237 (MQ=255)
                                                                                                                                                                                                                  ggcTGTTGGGGTGGATTGCCTTTAGTTGCTGGGAGGTTTCCCTCCAGGACTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCtgg                     <  2:368721/238‑3 (MQ=255)
                                                                                                                                                                                                                  ggcTGATGGGGTGTATGCCCTTGGGTTGTGGGGAGGTTTCCCTCCAGACCTGATTATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGGGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTg                     <  2:913933/238‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGGTTTGCTGGGAGGTTTCCCTCCAGGCCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCt  <  2:887941/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGAGTTGCTGGGGGGTTTCACTCCAGGCGTTATGATAAGTGATGTCGCGGTGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCATCCAGAACCAAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCt  <  2:244841/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCt  <  1:639464/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ccTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGGGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCgaga              <  1:1064342/227‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
AACATTTGACCAGAATTTTTATCTACACTTAAGTTAATTCTGACAGGCGCAGGTGGCAATAGCATGCCACTATTGAGTAAAGCCAGTCAGGGGAGAGAACATGACGCAGCCATTGGTCGGAAAACAGATTCTCATTGTTGAAGATGAGCAGGTATTTCGCTCGCTTCTGGATTCATGGTTTTCCTCATTGGGAGCGACAACGGTACTGGCGGCTGATGGGGTGGATGCCCTTGAGTTGCTGGGAGGTTTCACTCCAGACCTGATGATATGTGATATCGCGATGCCACGAATGAACGGGCTTAAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCT  >  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: