Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 208,433 G→C intergenic (‑39/+49) atsA ← / ← PA14_RS00940 arylsulfatase AtsA/ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463208,4330GC99.8% 2174.5 / NA 560intergenic (‑39/+49)atsA/PA14_RS00940arylsulfatase AtsA/ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (263/296);  total (264/296)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 560 total aligned reads displayed.

ATGGCCAGGGCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGTCCGTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTGCAGCAGCGCTTC  >  NC_008463/208285‑208577
                                                                                                                                                    |                                                                                                                                                
aTGGCCAGGGCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTc                                                                                                                                                >  1:8883685/1‑151 (MQ=255)
aTGGCCAGGGCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTc                                                                                                                                                >  2:2010349/1‑151 (MQ=255)
 tGGCCAGGGCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATTTCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCTTGCGCGGGTCTCGTCTGGATGGGAGGCGAGGGTGGGGCCCTCg                                                                                                                                               >  1:1062074/1‑151 (MQ=255)
  ggCCAGGGCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGc                                                                                                                                              >  2:7983263/1‑151 (MQ=255)
     cAGGGCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGcccc                                                                                                                                           >  2:2979087/1‑151 (MQ=255)
       gggCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCg                                                                                                                                         >  2:4843699/1‑151 (MQ=255)
         gCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCGCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGGGGGGCCCTCGCCCCCGcc                                                                                                                                       <  2:5686873/151‑1 (MQ=255)
          cGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGccc                                                                                                                                      >  2:3762208/1‑151 (MQ=255)
           gTCGTGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCTCCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCGGGATGGGAGGGGAGGGTGGGGCCCTCGCCCCCGCCCt                                                                                                                                     >  1:1364093/1‑151 (MQ=255)
            tCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCtt                                                                                                                                    <  1:938223/151‑1 (MQ=255)
             cGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGCGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTc                                                                                                                                   >  2:8658540/1‑151 (MQ=255)
              gAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTcc                                                                                                                                  >  1:693188/1‑151 (MQ=255)
              gAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTcc                                                                                                                                  <  2:6114122/151‑1 (MQ=255)
               aGGTTCGGCGTTGCGATCTCGCCGCCGAAGGCGCCGATATCGTAGAAGCCCCGGTCGTCGGCGACGATCACCAGTAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGACTGGCTGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTccc                                                                                                                                 <  1:869618/151‑1 (MQ=255)
               aGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTccc                                                                                                                                 >  1:10976595/1‑151 (MQ=255)
               aGGTGCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCCGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTccc                                                                                                                                 <  2:6218389/151‑1 (MQ=255)
                ggTTCGGCGTGGCGATCTCGCCGCCGAAGGCTCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGCGGGTGGGGCCCTCGCCCCCGCCCTTCCCg                                                                                                                                <  2:160669/151‑1 (MQ=255)
                ggTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCg                                                                                                                                <  1:2833304/151‑1 (MQ=255)
                   tCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAg                                                                                                                             <  2:4413079/151‑1 (MQ=255)
                   tCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAACTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAg                                                                                                                             >  1:10693958/1‑151 (MQ=255)
                    cGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAgg                                                                                                                            >  2:7012161/1‑151 (MQ=255)
                    cGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAgg                                                                                                                            >  1:8456009/1‑151 (MQ=255)
                     ggCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAggg                                                                                                                           <  1:1083545/151‑1 (MQ=255)
                     ggCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAggg                                                                                                                           >  1:6533846/1‑151 (MQ=255)
                      gCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAgggg                                                                                                                          >  1:11113725/1‑151 (MQ=255)
                        gTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCCTGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGa                                                                                                                        >  1:3957513/1‑151 (MQ=255)
                        gTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGa                                                                                                                        >  2:4632748/1‑151 (MQ=255)
                          ggCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGaaa                                                                                                                      >  1:4883477/1‑151 (MQ=255)
                            cGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAgg                                                                                                                    <  2:11037588/151‑1 (MQ=255)
                              aTCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGAGGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGACAGGGc                                                                                                                  >  2:10682014/1‑151 (MQ=255)
                               tctcGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGcc                                                                                                                 <  2:3445566/151‑1 (MQ=255)
                                 tcGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAg                                                                                                               <  2:4902964/151‑1 (MQ=255)
                                    ccgccgAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGc                                                                                                            <  2:7432603/151‑1 (MQ=255)
                                       ccgAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAg                                                                                                         <  2:5017950/151‑1 (MQ=255)
                                        cgAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAgg                                                                                                        >  2:2944306/1‑151 (MQ=255)
                                         gAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAggg                                                                                                       >  2:6345419/1‑151 (MQ=255)
                                          aaGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAgggg                                                                                                      <  1:3792706/151‑1 (MQ=255)
                                          aaGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAgggg                                                                                                      >  1:9046309/1‑151 (MQ=255)
                                          aaGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAgggg                                                                                                      >  1:6038551/1‑151 (MQ=255)
                                           aGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCAGCCGGACGGTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCGGCAGGGGGAAAGGGCCAGGGCAAGGGGc                                                                                                     >  2:3995237/1‑151 (MQ=255)
                                                    tatCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCa                                                                                            >  1:7209870/1‑151 (MQ=255)
                                                    tatCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCa                                                                                            <  2:78089/151‑1 (MQ=255)
                                                        ggAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCg                                                                                        >  1:3101064/1‑151 (MQ=255)
                                                          agaAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAg                                                                                      <  2:8843606/151‑1 (MQ=255)
                                                           gaAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGTAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAgg                                                                                     <  1:1470273/151‑1 (MQ=255)
                                                           gaAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCGTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAgg                                                                                     <  2:3730731/151‑1 (MQ=255)
                                                             aGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGaa                                                                                   >  1:3202783/1‑151 (MQ=255)
                                                                 cAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGACTTCg                                                                               >  2:2781452/1‑151 (MQ=255)
                                                                    gtcgtcGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAg                                                                            <  2:9228109/151‑1 (MQ=255)
                                                                      cgtcgGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAggg                                                                          >  2:5893666/1‑151 (MQ=255)
                                                                         cgGCGACGATCACCAGGACGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGGGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCGGGGCAAGGGGCTTTTCAGCCAGCGAGGAATTCGATGGGCAg                                                                       >  1:4225333/1‑151 (MQ=255)
                                                                            cgacgaTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGtt                                                                    <  1:33617/151‑1 (MQ=255)
                                                                             gacgaTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGttt                                                                   >  1:8290456/1‑151 (MQ=255)
                                                                             gacgaTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGtgt                                                                   >  1:2230008/1‑149 (MQ=255)
                                                                              acgaTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGtttt                                                                  <  2:8301/151‑1 (MQ=255)
                                                                                gaTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTcc                                                                <  1:6210552/151‑1 (MQ=255)
                                                                                gaTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTcc                                                                >  2:9577909/1‑151 (MQ=255)
                                                                                 ctCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCa                                                               <  2:5998883/150‑1 (MQ=255)
                                                                                 aTCACCAGGACCTTGGCGCGTTAGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTATCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCa                                                               <  2:5330983/151‑1 (MQ=255)
                                                                                 aTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCa                                                               <  1:1553445/151‑1 (MQ=255)
                                                                                 aTCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCa                                                               <  1:5473435/151‑1 (MQ=255)
                                                                                   cACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGACTTCGAAGGGCAGGTTTTCCAcc                                                             >  2:6870841/1‑151 (MQ=255)
                                                                                   cACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCAcc                                                             >  2:1223811/1‑151 (MQ=255)
                                                                                   cACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCACCGCCCTTCCCGCAGGGGGAAAGGTCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCAcc                                                             >  1:1242763/1‑151 (MQ=255)
                                                                                    aCCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCAccc                                                            >  1:5578198/1‑151 (MQ=255)
                                                                                     ccAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCg                                                           >  1:7709710/1‑151 (MQ=255)
                                                                                      cAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCAAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGt                                                          >  1:2551948/1‑151 (MQ=255)
                                                                                       aGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTc                                                         >  1:5232083/1‑151 (MQ=255)
                                                                                       aGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTc                                                         <  1:5518846/151‑1 (MQ=255)
                                                                                        ggcaGTGGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTcg                                                        <  2:6729243/148‑1 (MQ=255)
                                                                                           aGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCgcgg                                                     <  1:452093/151‑1 (MQ=255)
                                                                                            gtcggggCGAGTGCTCCTGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCgcggc                                                    <  1:9632454/148‑1 (MQ=255)
                                                                                              tGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCgcggcgg                                                  >  2:8170396/1‑151 (MQ=255)
                                                                                                gggCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCgcggcggcg                                                <  1:3630043/151‑1 (MQ=255)
                                                                                                 ggCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCgcggcggcgg                                               <  1:1746014/151‑1 (MQ=255)
                                                                                                 ggCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCgcggcggcgg                                               <  2:26015/151‑1 (MQ=255)
                                                                                                  gCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCgcggcggcggc                                              >  2:7596356/1‑151 (MQ=255)
                                                                                                    gTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGgcgc                                            >  1:1007734/1‑151 (MQ=255)
                                                                                                     tttGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAATCGAGGAATTCGAAGGGCAGGTTTTCCACCCGACGCGGCGGCGGcgcg                                           >  1:3015340/1‑151 (MQ=255)
                                                                                                       tGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGAc                                         <  1:9208115/151‑1 (MQ=255)
                                                                                                         cTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCACTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGt                                       <  2:425194/151‑1 (MQ=255)
                                                                                                            aTGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACt                                    <  1:5023156/151‑1 (MQ=255)
                                                                                                                 cgGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCg                               <  1:10346015/151‑1 (MQ=255)
                                                                                                                 cgGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCg                               <  1:4633492/151‑1 (MQ=255)
                                                                                                                 cgGGTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCg                               <  2:1536293/151‑1 (MQ=255)
                                                                                                                  gggTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGt                              >  2:3175319/1‑151 (MQ=255)
                                                                                                                   ggTCTCGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTc                             <  2:234333/151‑1 (MQ=255)
                                                                                                                     tctcGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTcgg                           <  2:9382214/151‑1 (MQ=255)
                                                                                                                        cGTCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGACTTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTcggcgg                        >  2:9011241/1‑151 (MQ=255)
                                                                                                                          tCTGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTc                      >  2:1046908/1‑151 (MQ=255)
                                                                                                                            tGGATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAg                    >  1:2935586/1‑151 (MQ=255)
                                                                                                                             ggATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGt                   <  2:10269283/151‑1 (MQ=255)
                                                                                                                             ggATGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGt                   >  1:9685056/1‑151 (MQ=255)
                                                                                                                                tGGGAGGAGAGGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCg                >  2:1934255/1‑151 (MQ=255)
                                                                                                                                     ggcgagGGTGGGGCCCTCGCCCCCGCCCGTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTGTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTgcag           <  2:2688528/148‑1 (MQ=255)
                                                                                                                                       agagGGTGGGGCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTgcagca         <  2:7280895/151‑1 (MQ=255)
                                                                                                                                          gggtggggCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTGCAGCAgcg      <  2:2532662/151‑1 (MQ=255)
                                                                                                                                          gggtggggCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGATTTTCAGCAAGCGAGGAATTCGTAGTGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTGCATCAgcg      <  2:10038348/151‑1 (MQ=255)
                                                                                                                                           ggtggggCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTGCAGCAgcgc     >  1:7431758/1‑151 (MQ=255)
                                                                                                                                             tggggCCCTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTGCAGCAGCGCtt   <  1:6744493/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATGGCCAGGGCGTCGAGGTTCGGCGTGGCGATCTCGCCGCCGAAGGCGCCGATATCGGAGAAGCCCAGGTCGTCGGCGACGATCACCAGGAAGTTGGGGCGTTTGCTCATGCGCGGGTCTCGTCTGGATGGGAGGAGAGGGTGGGTCCGTCGCCCCCGCCCTTCCCGCAGGGGGAAAGGGCCAGGGCAAGGGGCTTTTCAGCAAGCGAGGAATTCGAAGGGCAGGTTTTCCACCCGTCGCGGCGGCGGCGCGACGTACTCGCCGTCGGCGGTCAGTTCGTGCAGCAGCGCTTC  >  NC_008463/208285‑208577

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