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.8% 1658.6 / NA 424P426A (CCC→GCC) phaCclass II poly(R)‑hydroxyalkanoic acid synthase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base G (233/189);  minor base T (0/1);  total (233/191)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 428 total aligned reads displayed.

GAACTACTGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCCCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCCTCT  >  NC_008463/5967877‑5968172
                                                                                                                                                    |                                                                                                                                                   
  actactGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCt                                                                                                                                                 >  1:3799378/1‑151 (MQ=255)
  actactGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGACGAACTGGTCAAGAGCAACCCGCTGAAACGCCCCGGCGCCCt                                                                                                                                                 >  1:9852569/1‑151 (MQ=255)
   ctactGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTg                                                                                                                                                <  1:2353090/151‑1 (MQ=255)
    tactGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTgg                                                                                                                                               >  1:9469474/1‑151 (MQ=255)
     actGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGa                                                                                                                                              <  1:6303712/151‑1 (MQ=255)
        gggTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGt                                                                                                                                           >  2:2287447/1‑151 (MQ=255)
         ggTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGtc                                                                                                                                          <  2:8176680/151‑1 (MQ=255)
         ggTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGtc                                                                                                                                          <  1:11361009/151‑1 (MQ=255)
             aacaacTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCg                                                                                                                                      <  2:10217427/151‑1 (MQ=255)
             aacaacTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCg                                                                                                                                      >  1:2627237/1‑151 (MQ=255)
              acaacTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCgg                                                                                                                                     <  1:8779323/151‑1 (MQ=255)
              acaacTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGAACGGCGAGTTGGACGAACTGTTCAAGAGCACCCCTCAGACCCTCCCCCGCGCCCTCGAGGTCTCCgg                                                                                                                                     >  2:10100147/1‑151 (MQ=18)
                 acTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCAc                                                                                                                                  <  1:5627703/151‑1 (MQ=255)
                       tGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCAt                                                                                                                            <  1:7391746/151‑1 (MQ=255)
                        gCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATc                                                                                                                           <  2:5524791/151‑1 (MQ=255)
                              cAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTg                                                                                                                     <  1:10428043/151‑1 (MQ=255)
                                 ccAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGTCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAg                                                                                                                  <  2:6756786/151‑1 (MQ=255)
                                  cAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGc                                                                                                                 <  1:8166867/151‑1 (MQ=255)
                                   aGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCa                                                                                                                <  2:577920/151‑1 (MQ=255)
                                     ccgccgGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAgg                                                                                                              >  2:1627917/1‑151 (MQ=255)
                                      cgccgGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGt                                                                                                             >  2:7811107/1‑151 (MQ=255)
                                      cgccgGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGt                                                                                                             >  2:6822835/1‑151 (MQ=255)
                                        ccgGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGa                                                                                                           <  1:10751530/151‑1 (MQ=255)
                                        ccgGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGa                                                                                                           <  1:4099933/151‑1 (MQ=255)
                                          ggCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACt                                                                                                         >  2:5662876/1‑151 (MQ=255)
                                          ggCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACt                                                                                                         >  2:8911817/1‑151 (MQ=255)
                                           gCGGTCGCCCTCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACtt                                                                                                        <  2:6127513/151‑1 (MQ=255)
                                             gTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGc                                                                                                      <  1:3058102/151‑1 (MQ=255)
                                             gTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGc                                                                                                      <  1:9360490/151‑1 (MQ=255)
                                             gTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGc                                                                                                      >  1:4591547/1‑151 (MQ=255)
                                              ttCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCg                                                                                                     <  2:1725207/151‑1 (MQ=255)
                                              ttCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCg                                                                                                     <  1:9534358/151‑1 (MQ=255)
                                                cGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCgcc                                                                                                   >  2:5429571/1‑149 (MQ=255)
                                                cGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGAc                                                                                                   >  1:8665370/1‑151 (MQ=255)
                                                cGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGAc                                                                                                   >  1:3431777/1‑151 (MQ=255)
                                                   cATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTc                                                                                                >  2:59935/1‑151 (MQ=255)
                                                    aTCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCt                                                                                               >  2:3280027/1‑151 (MQ=255)
                                                        tctACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtg                                                                                           <  2:10770051/151‑1 (MQ=255)
                                                        tctACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtg                                                                                           >  2:8865436/1‑151 (MQ=255)
                                                         ctACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgt                                                                                          >  2:192806/1‑151 (MQ=255)
                                                         ctACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgt                                                                                          >  1:1206584/1‑151 (MQ=255)
                                                          tACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgtg                                                                                         >  2:9399772/1‑151 (MQ=255)
                                                          tACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACtgtg                                                                                         <  1:10045491/151‑1 (MQ=255)
                                                              ggAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGc                                                                                     <  2:7334590/151‑1 (MQ=255)
                                                              ggAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGc                                                                                     >  2:6660040/1‑151 (MQ=255)
                                                               gAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGcc                                                                                    <  2:8670360/151‑1 (MQ=255)
                                                                aacaacGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGACGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCg                                                                                   >  1:5469234/1‑151 (MQ=255)
                                                                  caacGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGt                                                                                 >  1:7736112/1‑151 (MQ=255)
                                                                   aacGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTc                                                                                <  2:852821/151‑1 (MQ=255)
                                                                      gACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGa                                                                             <  1:7120529/151‑1 (MQ=255)
                                                                       acacCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGaa                                                                            >  1:4842152/1‑151 (MQ=255)
                                                                         accacGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACg                                                                          >  1:8043066/1‑151 (MQ=255)
                                                                         accacGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACg                                                                          <  1:10585189/151‑1 (MQ=255)
                                                                         accacGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACg                                                                          >  2:11556803/1‑151 (MQ=255)
                                                                          ccacGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGa                                                                         >  2:9417226/1‑151 (MQ=255)
                                                                             cgcgCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACca                                                                      <  1:1627105/151‑1 (MQ=255)
                                                                              gcgcCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACcac                                                                     <  1:6703017/151‑1 (MQ=255)
                                                                              gcgcCTGCCCGCCGCGCTGCACGGCGAGTTCGACGAACTGTTCCAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCTGTTGACTTGAGACTTGAAATGTGTCGCCGGTCTGACCGACcac                                                                     >  1:2010761/1‑151 (MQ=255)
                                                                                   tGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAc                                                                >  2:9746262/1‑151 (MQ=255)
                                                                                     cccGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAccc                                                              >  2:1327408/1‑151 (MQ=255)
                                                                                      ccgccgCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAcccc                                                             >  1:2785616/1‑151 (MQ=255)
                                                                                      ccgccgCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAcccc                                                             >  2:2932667/1‑151 (MQ=255)
                                                                                       cgccgcGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCAccccc                                                            <  1:4268481/151‑1 (MQ=255)
                                                                                        gccgcGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCt                                                           <  2:2406703/151‑1 (MQ=255)
                                                                                         ccgcGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTg                                                          <  1:4520909/151‑1 (MQ=255)
                                                                                           gcgcTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTggg                                                        <  2:10478098/151‑1 (MQ=255)
                                                                                             gcTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAg                                                      >  1:783648/1‑151 (MQ=255)
                                                                                                 cACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGt                                                  <  2:11323573/151‑1 (MQ=255)
                                                                                                  aCGGCGAGTTCGTCGACCTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTg                                                 <  2:10605342/151‑1 (MQ=255)
                                                                                                  aCGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTg                                                 >  1:2906335/1‑151 (MQ=255)
                                                                                                        aGTTCGTCGAACTGTTCAAGAGCAACCCGCTGACCCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACaa                                           <  2:5210532/151‑1 (MQ=255)
                                                                                                            cgtcgAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCg                                       <  1:2261267/151‑1 (MQ=255)
                                                                                                               cgAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGcc                                    >  1:3352749/1‑151 (MQ=255)
                                                                                                                 aaCTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAg                                  <  2:4642396/151‑1 (MQ=255)
                                                                                                                   cTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCTCCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGc                                >  1:10996625/1‑151 (MQ=255)
                                                                                                                      ttCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTAc                             <  2:2215994/151‑1 (MQ=255)
                                                                                                                      ttCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCGGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTAc                             <  2:4242017/151‑1 (MQ=255)
                                                                                                                       tCAAGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACt                            <  1:7015478/151‑1 (MQ=255)
                                                                                                                         aaGAGCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTgg                          <  2:9019828/151‑1 (MQ=255)
                                                                                                                           gagCAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGACGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACAGggc                        >  2:1458211/1‑151 (MQ=255)
                                                                                                                           gagCAACCCGCTGAACCGCCCCGGAGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCAGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGggc                        >  2:8346882/1‑151 (MQ=255)
                                                                                                                              cAACCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGggcggc                     >  2:11188610/1‑151 (MQ=255)
                                                                                                                               aaCCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCa                    <  2:2627237/151‑1 (MQ=255)
                                                                                                                               aaCCCGCTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATAACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCa                    >  1:5805041/1‑151 (MQ=255)
                                                                                                                                   cGCAGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTg                >  1:1410609/1‑151 (MQ=255)
                                                                                                                                     cTGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTCCTGGGCGGCAAGTGCg              >  2:7157341/1‑151 (MQ=255)
                                                                                                                                      tGAACCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGa             >  1:1807203/1‑151 (MQ=255)
                                                                                                                                        aaCCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCCCCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGt           >  2:7006410/1‑151 (MQ=255)
                                                                                                                                        aaCCGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGt           <  1:6188723/151‑1 (MQ=255)
                                                                                                                                          ccGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTc         <  1:5290922/151‑1 (MQ=255)
                                                                                                                                           cGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTCCTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCa        <  2:5031886/151‑1 (MQ=255)
                                                                                                                                           cGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCa        <  1:8584640/151‑1 (MQ=255)
                                                                                                                                           cGCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCa        >  2:8960370/1‑151 (MQ=255)
                                                                                                                                            gCCCCGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCAt       <  2:6443327/151‑1 (MQ=255)
                                                                                                                                              cccGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATcc     >  2:10962340/1‑151 (MQ=255)
                                                                                                                                               ccGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCTTCct    >  2:6500971/1‑151 (MQ=255)
                                                                                                                                               ccGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCct    >  2:2818340/1‑151 (MQ=255)
                                                                                                                                               ccGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCct    >  2:10143021/1‑151 (MQ=255)
                                                                                                                                                cGGCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCctc   <  1:1127239/151‑1 (MQ=255)
                                                                                                                                                 ggCGCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCctct  <  2:6336190/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAACTACTGGGTCAACAACTACCTGCTCGGCAACCAGCCGCCGGCGTTCGACATCCTCTACTGGAACAACGACACCACGCGCCTGCCCGCCGCGCTGCACGGCGAGTTCGTCGAACTGTTCAAGAGCAACCCGCTGAACCGCCCCGGCCCCCTGGAGGTCTCCGGCACGCCCATCGACCTGAAGCAGGTGACTTGCGACTTCTACTGTGTCGCCGGTCTGAACGACCACATCACCCCCTGGGAGTCGTGCTACAAGTCGGCCAGGCTACTGGGCGGCAAGTGCGAGTTCATCCTCT  >  NC_008463/5967877‑5968172

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: ‑