Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,231,349 A→C 100% Q227P (CAG→CCG)  ygjR → putative NAD(P)‑dependent dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,231,3490AC96.4% 644.0 / 3.4 162Q227P (CAG→CCG) ygjRputative NAD(P)‑dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base C (68/88);  total (70/92)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 165 total aligned reads displayed.

ATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCATAT  >  CP009273/3231205‑3231496
                                                                                                                                                |                                                                                                                                                   
 tatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCAc                                                                                                                                                 >  1:593120/1‑148 (MQ=255)
    cGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGGTTACGGTGATTTCAGCGTCACCTTTCCGCACTCa                                                                                                                                              >  1:1808288/1‑148 (MQ=255)
    cGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGGGGTGATGGATTACGGTGATTTCAGCGTCACCATGCCGCACTCa                                                                                                                                              >  1:1808512/1‑148 (MQ=255)
       tttttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaaaa                                                                                                                                           >  1:294833/1‑148 (MQ=255)
       tttttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaaaa                                                                                                                                           >  1:3059574/1‑148 (MQ=255)
       tttttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaaaa                                                                                                                                           >  2:900024/1‑148 (MQ=255)
               gTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGa                                                                                                                                   <  2:293504/148‑1 (MQ=255)
               gTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGa                                                                                                                                   <  1:1424844/148‑1 (MQ=255)
                tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGAt                                                                                                                                  <  2:497297/148‑1 (MQ=255)
                      gTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTc                                                                                                                            <  1:1124154/148‑1 (MQ=255)
                        cggcggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCt                                                                                                                          <  2:1099630/148‑1 (MQ=255)
                         ggcggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTg                                                                                                                         <  1:1889581/148‑1 (MQ=255)
                         ggcggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTg                                                                                                                         <  2:102788/148‑1 (MQ=255)
                              tGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcgag                                                                                                                    >  1:2966265/1‑148 (MQ=255)
                               ggCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcgagc                                                                                                                   <  2:3455494/148‑1 (MQ=255)
                                gCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcgagcg                                                                                                                  >  2:1474746/1‑148 (MQ=255)
                                 cGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcgagcga                                                                                                                 <  1:2446425/148‑1 (MQ=255)
                                  gTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGcgagcgag                                                                                                                <  1:1745650/148‑1 (MQ=255)
                                   ttATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGTCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCgaga                                                                                                               <  2:239796/148‑1 (MQ=255)
                                   ttATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCgaga                                                                                                               <  1:1638395/148‑1 (MQ=255)
                                   ttATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCgaga                                                                                                               >  2:226125/1‑148 (MQ=255)
                                    tATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGAt                                                                                                              <  2:177283/148‑1 (MQ=255)
                                    tATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGGGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGAt                                                                                                              >  2:1051512/1‑148 (MQ=255)
                                       ttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                           >  1:2654788/1‑148 (MQ=255)
                                       ttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                           <  1:735163/148‑1 (MQ=255)
                                       ttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                           <  1:735055/148‑1 (MQ=255)
                                       ttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                           >  2:1946943/1‑148 (MQ=255)
                                       ttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                           >  1:3202814/1‑148 (MQ=255)
                                         ggTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAgg                                                                                                         >  2:596243/1‑148 (MQ=255)
                                         ggTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAgg                                                                                                         >  2:596235/1‑148 (MQ=255)
                                          gTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAggg                                                                                                        <  2:3407777/148‑1 (MQ=255)
                                            gAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCg                                                                                                      >  2:1339657/1‑148 (MQ=255)
                                             aGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGGGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGa                                                                                                     >  2:1835099/1‑148 (MQ=255)
                                               ccGAAAAGTGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTACCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAg                                                                                                   <  1:2682501/148‑1 (MQ=255)
                                               ccGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAg                                                                                                   <  2:2889669/148‑1 (MQ=255)
                                                  aaaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAg                                                                                                <  1:2764479/148‑1 (MQ=255)
                                                   aaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                               <  2:1339732/148‑1 (MQ=255)
                                                   aaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                               >  2:1185617/1‑148 (MQ=255)
                                                   aaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                               >  1:461020/1‑148 (MQ=255)
                                                       cGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCg                                                                                           <  2:3202814/148‑1 (MQ=255)
                                                       cGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCg                                                                                           <  2:1637917/148‑1 (MQ=255)
                                                        gTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGc                                                                                          <  1:2492935/148‑1 (MQ=255)
                                                         tGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCt                                                                                         >  2:2957372/1‑148 (MQ=255)
                                                          gCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTg                                                                                        >  2:3019934/1‑148 (MQ=255)
                                                           cAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTgg                                                                                       <  2:844318/148‑1 (MQ=255)
                                                           cAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTgg                                                                                       <  2:1374308/148‑1 (MQ=255)
                                                             ggCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTg                                                                                     >  1:2800056/1‑148 (MQ=255)
                                                              gCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGa                                                                                    <  2:2240064/148‑1 (MQ=255)
                                                              gCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGa                                                                                    >  2:1216049/1‑148 (MQ=255)
                                                               cAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGAt                                                                                   <  1:3103729/148‑1 (MQ=255)
                                                               cAACCGCCAGTGTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGAt                                                                                   <  2:1905413/148‑1 (MQ=255)
                                                                 aCCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTg                                                                                 <  1:982539/148‑1 (MQ=255)
                                                                 aCCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTg                                                                                 >  1:154352/1‑148 (MQ=255)
                                                                  ccGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGa                                                                                <  1:2514380/148‑1 (MQ=255)
                                                                  ccGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGa                                                                                >  1:2339452/1‑148 (MQ=255)
                                                                        gTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACt                                                                          <  2:2997334/148‑1 (MQ=255)
                                                                          ttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGt                                                                        <  1:377942/148‑1 (MQ=255)
                                                                            gCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCt                                                                      <  1:2143422/148‑1 (MQ=255)
                                                                             cTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTg                                                                     >  2:1397708/1‑148 (MQ=255)
                                                                                gCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAAt                                                                  >  1:3523880/1‑148 (MQ=255)
                                                                                 cAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATg                                                                 <  1:435308/148‑1 (MQ=255)
                                                                                 cAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATg                                                                 <  1:435924/148‑1 (MQ=255)
                                                                                 cAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATg                                                                 <  2:153425/148‑1 (MQ=255)
                                                                                  aaGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGc                                                                <  2:3383297/148‑1 (MQ=255)
                                                                                   aGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGcc                                                               <  2:348107/148‑1 (MQ=255)
                                                                                   aGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGcc                                                               >  1:3449637/1‑148 (MQ=255)
                                                                                       gcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGaa                                                           <  2:844492/148‑1 (MQ=255)
                                                                                       gcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGaa                                                           <  2:844157/148‑1 (MQ=255)
                                                                                        cgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGaaa                                                          >  2:32504/1‑148 (MQ=255)
                                                                                           tGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgtg                                                       >  2:688964/1‑148 (MQ=255)
                                                                                            gACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgtgt                                                      <  2:3439585/148‑1 (MQ=255)
                                                                                             aCGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGtgtg                                                     <  2:1009328/148‑1 (MQ=255)
                                                                                             aCGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGtgtg                                                     <  2:251196/148‑1 (MQ=255)
                                                                                                cccAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCtt                                                  >  1:528479/1‑148 (MQ=255)
                                                                                                    gggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTg                                              <  2:766648/148‑1 (MQ=255)
                                                                                                     ggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGc                                             >  2:2566156/1‑148 (MQ=255)
                                                                                                     ggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGc                                             <  1:2888530/148‑1 (MQ=255)
                                                                                                         tggtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCcgcg                                         <  1:707129/148‑1 (MQ=255)
                                                                                                         tggtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCcgcg                                         <  1:704414/148‑1 (MQ=255)
                                                                                                          ggtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGt                                        <  2:1106483/148‑1 (MQ=255)
                                                                                                            tggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTgg                                      >  2:2152906/1‑148 (MQ=255)
                                                                                                             ggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGc                                     >  1:1019137/1‑148 (MQ=255)
                                                                                                             ggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGc                                     >  1:1019130/1‑148 (MQ=255)
                                                                                                                   ggATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCaa                               >  1:2331155/1‑148 (MQ=255)
                                                                                                                    gATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCaaa                              >  2:1958471/1‑148 (MQ=255)
                                                                                                                     aTTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAAt                             >  1:3299242/1‑148 (MQ=255)
                                                                                                                      ttACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATg                            >  1:1151236/1‑148 (MQ=255)
                                                                                                                       tACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGc                           >  1:1438872/1‑148 (MQ=255)
                                                                                                                          ggTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAgg                        >  2:1184266/1‑148 (MQ=255)
                                                                                                                          ggTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAgg                        >  1:2361295/1‑148 (MQ=255)
                                                                                                                           gTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGa                       >  2:2702328/1‑148 (MQ=255)
                                                                                                                            tGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAt                      <  2:822794/148‑1 (MQ=255)
                                                                                                                             gATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAtc                     >  1:727424/1‑148 (MQ=255)
                                                                                                                              aTTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAtct                    <  1:1327817/148‑1 (MQ=255)
                                                                                                                               tttCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGCTCGCGGGTGATTGACAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCCAAGGCAGGAtctc                   >  2:3183410/1‑148 (MQ=255)
                                                                                                                                ttCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCa                  <  1:158254/148‑1 (MQ=255)
                                                                                                                                ttCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCa                  <  1:2050084/148‑1 (MQ=255)
                                                                                                                                ttCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCa                  >  2:1063924/1‑148 (MQ=255)
                                                                                                                                         aCCTTGCCGCACTCAAAAGTCCGTGAGTCTGTCCTGGCGTGCGTGTTTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAATGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgc         <  2:137058/148‑1 (MQ=255)
                                                                                                                                                agcaCTCAAAAGTCATTGATTCTGTCCTGGCGAGCGAGATTAAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCatat  <  1:3532965/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCATAT  >  CP009273/3231205‑3231496

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