Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,968,025 C→G P426A (CCC→GCC)  phaC → class II poly(R)‑hydroxyalkanoic acid synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,968,0250CG99.6% 1762.5 / NA 461P426A (CCC→GCC) phaCclass II poly(R)‑hydroxyalkanoic acid synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (236/222);  minor base T (0/2);  total (237/224)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 463 total aligned reads displayed.

ACTACTGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCCCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCCTC  >  NC_008463/5967879‑5968171
                                                                                                                                                  |                                                                                                                                                  
actactGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGcgccc                                                                                                                                                 >  1:8413847/1‑150 (MQ=255)
    ctGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCTAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGa                                                                                                                                             <  2:4766601/150‑1 (MQ=255)
       ggTCAACAACTACCGGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTCATGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGt                                                                                                                                          <  2:4975832/150‑1 (MQ=255)
           aacaacTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTcc                                                                                                                                      >  2:3170316/1‑150 (MQ=255)
             caacTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCgg                                                                                                                                    <  1:2109704/150‑1 (MQ=255)
             caacTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCgg                                                                                                                                    <  2:9060844/150‑1 (MQ=255)
              aacGACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGc                                                                                                                                   <  2:2442343/150‑1 (MQ=255)
               acTCCCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCa                                                                                                                                  <  2:2836642/150‑1 (MQ=255)
               acTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCa                                                                                                                                  >  1:4219711/1‑150 (MQ=255)
                cTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCAc                                                                                                                                 >  1:6790128/1‑150 (MQ=255)
                   ccTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGcc                                                                                                                              <  2:7880644/150‑1 (MQ=255)
                      gCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCAt                                                                                                                           >  2:10135050/1‑150 (MQ=255)
                       cTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATc                                                                                                                          <  1:10960352/150‑1 (MQ=255)
                        tCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCg                                                                                                                         >  2:8786780/1‑150 (MQ=255)
                              accatgcgccgGAGTTCGACATCATCTACTGGAAGAACGACAACACGCTCATGCCCGGCGAGCTGCACGTCGAGTTCGTCGAAATGTTCAAGAGCAACCCGCTGTACCTCCACTTCTCCCTGGAGGTCTCCGGCACGCCCATCGACCTGa                                                                                                                   <  1:3094629/144‑1 (MQ=255)
                              aCCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGa                                                                                                                   >  1:9135156/1‑150 (MQ=255)
                                  gccgccGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCa                                                                                                               >  1:10491878/1‑150 (MQ=255)
                                  gccgccGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCa                                                                                                               >  1:7149981/1‑150 (MQ=255)
                                         gCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACt                                                                                                        >  1:8829160/1‑150 (MQ=255)
                                         gCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACt                                                                                                        <  1:3655676/150‑1 (MQ=255)
                                         gCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACt                                                                                                        >  1:504677/1‑150 (MQ=255)
                                         gCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCCGCCCGCCGCGCTGCACGTCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACt                                                                                                        <  2:4165910/150‑1 (MQ=255)
                                          cGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACtt                                                                                                       >  2:4223714/1‑150 (MQ=255)
                                              cGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGa                                                                                                   <  1:1398613/150‑1 (MQ=255)
                                                aCATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACt                                                                                                 <  1:9632671/150‑1 (MQ=255)
                                                  aTCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTc                                                                                               <  2:7674811/150‑1 (MQ=255)
                                                  aTCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTc                                                                                               >  2:4520631/1‑150 (MQ=255)
                                                     ctctACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGAGTTCTAc                                                                                            >  2:949880/1‑150 (MQ=255)
                                                     ctctACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTAc                                                                                            >  2:3754228/1‑150 (MQ=255)
                                                       ctACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtg                                                                                          <  2:1036078/150‑1 (MQ=255)
                                                        tACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgt                                                                                         <  1:5583304/150‑1 (MQ=255)
                                                        tACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgt                                                                                         <  1:9710785/150‑1 (MQ=255)
                                                         aCTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgtg                                                                                        >  1:7728902/1‑150 (MQ=255)
                                                         aCTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgtg                                                                                        <  1:3449959/150‑1 (MQ=255)
                                                          cTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTgtgt                                                                                       >  2:4444362/1‑150 (MQ=255)
                                                                 aacGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGt                                                                                <  2:10050559/150‑1 (MQ=255)
                                                                    gACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTg                                                                             >  2:7445588/1‑150 (MQ=255)
                                                                      caccacGCGCCGGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGaa                                                                           <  1:10785217/150‑1 (MQ=255)
                                                                        ccacGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACg                                                                         >  1:628070/1‑150 (MQ=255)
                                                                          acGCGCCTGCCCGCCGCGCTGCACGGCGAGGTCGTCGAACTGTTCAAGCGCAACCCGCTGAACCGCCCCGGCGCCCGGGAGGTCTGCGGCACGCCCAGCGTGCTGACGCCGTGGACTGGCGCCTGCGACTGTGTCGCCGGTCTGAGCGag                                                                       >  2:5668694/1‑149 (MQ=255)
                                                                          acGCGCCGGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGAc                                                                       <  2:4041110/150‑1 (MQ=255)
                                                                           cgcgCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGAcc                                                                      >  1:9287127/1‑150 (MQ=255)
                                                                               ccTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTTTTCAAGAGCAACACGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCTACCTGAAGCAGGTGACTTGCGACTTCAACTGTGTCGCCGGTCTGAACGACcacaa                                                                  >  1:3529327/1‑149 (MQ=255)
                                                                               ccTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACAt                                                                  >  2:1733827/1‑150 (MQ=255)
                                                                               ccTGCCCGCCGCGCTGCACGGCGAGTTCGGCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCGGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGAGACTTGCGACTACTACTGTGTCGCCGGTCTGAACGACCACAt                                                                  >  2:8240223/1‑150 (MQ=255)
                                                                                 tGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCa                                                                <  1:521347/150‑1 (MQ=255)
                                                                                 tGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCACCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCa                                                                >  2:6435721/1‑150 (MQ=255)
                                                                                  gCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAc                                                               <  2:7269119/150‑1 (MQ=255)
                                                                                  gCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAc                                                               <  1:1735148/150‑1 (MQ=255)
                                                                                  gCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTATACTGTGTCGCCGGTCTGAACGACCACATCAc                                                               <  2:5455193/150‑1 (MQ=255)
                                                                                     cgccgcGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACGGTGTCGCCGGTCTGAACGACCACATCAcccc                                                            <  2:9614086/150‑1 (MQ=255)
                                                                                      gccgcGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAccccc                                                           >  2:6757493/1‑150 (MQ=255)
                                                                                         gcgcTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTgg                                                        >  1:10289606/1‑150 (MQ=255)
                                                                                             tGCACGTCGAGTTAGTCGAACTGTTCAAGAGCAACCAGGTTATCCGCCACGGTGCCCTGGAGGTCTGCGGCTCTCCAATAGTCGTGAAGATGTTGTATTGCGACTTCTACTGTGTAGCAGGTCTGATCGAACTGATCACCCGCTGGgtgt                                                    <  2:2231729/150‑1 (MQ=255)
                                                                                             tGCACGGCGAGTTCGTCGAACTGATCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGt                                                    <  1:5413477/150‑1 (MQ=255)
                                                                                               cACGGCGAGTTCGTCGAACTGATCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCg                                                  <  1:226624/150‑1 (MQ=255)
                                                                                                    cGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTa                                             <  1:4855124/150‑1 (MQ=255)
                                                                                                      aGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACa                                           >  2:1298018/1‑150 (MQ=255)
                                                                                                      aGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACa                                           >  2:6659942/1‑150 (MQ=255)
                                                                                                       gTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACaa                                          >  1:4523857/1‑150 (MQ=255)
                                                                                                        ttCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAg                                         <  2:5223800/150‑1 (MQ=255)
                                                                                                         tcgtcgACCTGTTCAAGAGCAACCCGCTGCACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCGCCTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGt                                        <  2:6133816/150‑1 (MQ=255)
                                                                                                         tcgtcgAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGt                                        <  1:2405885/150‑1 (MQ=255)
                                                                                                         tcgtcgAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGt                                        <  2:7345058/150‑1 (MQ=255)
                                                                                                         tcgtcgAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGt                                        <  1:6197086/150‑1 (MQ=255)
                                                                                                           gtcgAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCg                                      >  1:934446/1‑150 (MQ=255)
                                                                                                             cgAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGc                                    >  1:7485760/1‑150 (MQ=255)
                                                                                                              gAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGcc                                   >  1:1215828/1‑150 (MQ=255)
                                                                                                                 cTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAgg                                >  2:10932551/1‑150 (MQ=255)
                                                                                                                  tGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGc                               <  2:10685892/150‑1 (MQ=255)
                                                                                                                  tGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGc                               <  2:8313942/150‑1 (MQ=255)
                                                                                                                   gTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCt                              <  2:10225632/150‑1 (MQ=255)
                                                                                                                    ttCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTa                             <  1:2920668/150‑1 (MQ=255)
                                                                                                                     tCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCGAc                            >  2:4895573/1‑150 (MQ=255)
                                                                                                                      cAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGACGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACt                           <  2:1489/150‑1 (MQ=255)
                                                                                                                       aaGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTg                          <  1:10055247/150‑1 (MQ=255)
                                                                                                                       aaGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTg                          <  1:3039741/150‑1 (MQ=255)
                                                                                                                       aaGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTg                          <  2:2063150/150‑1 (MQ=255)
                                                                                                                        agagCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTgg                         <  2:8596291/150‑1 (MQ=255)
                                                                                                                        agagCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTgg                         <  2:8914735/150‑1 (MQ=255)
                                                                                                                         gagCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTggg                        <  1:3230130/150‑1 (MQ=255)
                                                                                                                           gCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACTAGTCGGCCAGGCTACTGggcg                      >  1:5699995/1‑150 (MQ=255)
                                                                                                                           gCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGggcg                      <  1:398230/150‑1 (MQ=255)
                                                                                                                            cAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGggcgg                     >  1:6459146/1‑150 (MQ=255)
                                                                                                                             aaCCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGggcggc                    >  2:9405625/1‑150 (MQ=255)
                                                                                                                             aaCCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGggcggc                    >  2:8164849/1‑150 (MQ=255)
                                                                                                                             aaCCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGggcggc                    >  2:8040619/1‑150 (MQ=255)
                                                                                                                               cccGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCaa                  <  2:10292481/150‑1 (MQ=255)
                                                                                                                                  gCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTg               <  1:7316065/150‑1 (MQ=255)
                                                                                                                                   cTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGc              >  1:2663351/1‑150 (MQ=255)
                                                                                                                                    tGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCg             >  2:3286923/1‑150 (MQ=255)
                                                                                                                                    tGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCg             <  1:7961611/150‑1 (MQ=255)
                                                                                                                                     gAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGa            >  2:8289313/1‑150 (MQ=255)
                                                                                                                                     gAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGa            >  1:9033075/1‑150 (MQ=255)
                                                                                                                                      aaCCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAg           <  2:59517/150‑1 (MQ=255)
                                                                                                                                        ccGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGtt         <  1:7637416/150‑1 (MQ=255)
                                                                                                                                           ccccGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCAt      >  2:2984385/1‑150 (MQ=255)
                                                                                                                                             ccGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATcc    <  1:9914629/150‑1 (MQ=255)
                                                                                                                                               ggCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAGGTGCGAGTTCATCctc  >  1:4596537/1‑150 (MQ=255)
                                                                                                                                               ggCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCctc  <  2:1924270/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACTACTGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCCCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCCTC  >  NC_008463/5967879‑5968171

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