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,3490AC97.9% 1179.6 / 3.4 292Q227P (CAG→CCG) ygjRputative NAD(P)‑dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base C (128/157);  total (130/162)
Fisher's exact test for biased strand distribution p-value = 4.67e-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 301 total aligned reads displayed.

GGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCAT  >  CP009273/3231203‑3231494
                                                                                                                                                  |                                                                                                                                                 
ggATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCCGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTgccgc                                                                                                                                                 >  1:2249266/1‑149 (MQ=255)
ggATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTgccgc                                                                                                                                                 >  1:3675581/1‑149 (MQ=255)
ggATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTgccgc                                                                                                                                                 >  1:1427905/1‑149 (MQ=255)
  atatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCAc                                                                                                                                               <  2:3022890/149‑1 (MQ=255)
  atatCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCAc                                                                                                                                               <  2:2622449/149‑1 (MQ=255)
  atatCGATTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCAc                                                                                                                                               <  2:4551017/149‑1 (MQ=255)
     tCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCCGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCa                                                                                                                                            >  1:2510543/1‑149 (MQ=255)
     tCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCa                                                                                                                                            >  2:3617384/1‑149 (MQ=255)
      cGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaa                                                                                                                                           <  1:4148841/149‑1 (MQ=255)
        gTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaaaa                                                                                                                                         >  1:3033452/1‑149 (MQ=255)
         tttttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAg                                                                                                                                        >  1:1025040/1‑149 (MQ=255)
            ttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCa                                                                                                                                     >  1:3163211/1‑149 (MQ=255)
            ttACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCa                                                                                                                                     >  1:2837856/1‑149 (MQ=255)
                tGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGa                                                                                                                                 >  1:2749904/1‑149 (MQ=255)
                tGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGGGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGa                                                                                                                                 >  1:3057563/1‑149 (MQ=255)
                 gTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGAt                                                                                                                                >  2:1363764/1‑149 (MQ=255)
                  tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGAtt                                                                                                                               >  2:2386264/1‑149 (MQ=255)
                  tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGAtt                                                                                                                               >  1:888239/1‑149 (MQ=255)
                   cTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTc                                                                                                                              <  1:2177855/149‑1 (MQ=255)
                   cTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTc                                                                                                                              <  2:4323244/149‑1 (MQ=255)
                   cTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTc                                                                                                                              >  2:2363288/1‑149 (MQ=255)
                    tGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCt                                                                                                                             <  1:1492463/149‑1 (MQ=255)
                     ggCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCCGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTg                                                                                                                            >  1:4187588/1‑149 (MQ=255)
                      gCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGt                                                                                                                           <  2:334570/149‑1 (MQ=255)
                        gTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTcc                                                                                                                         <  2:2789621/149‑1 (MQ=255)
                         tCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCt                                                                                                                        <  1:4592049/149‑1 (MQ=255)
                          cggcggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTg                                                                                                                       <  1:4709264/149‑1 (MQ=255)
                            gcggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgc                                                                                                                     <  2:1844442/149‑1 (MQ=255)
                                 ggCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcgagcg                                                                                                                >  2:2931650/1‑149 (MQ=255)
                                   cGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGcgagcgag                                                                                                              <  2:2900032/149‑1 (MQ=255)
                                       aTTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTc                                                                                                          >  2:1566553/1‑149 (MQ=255)
                                          tGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAgg                                                                                                       <  1:3502379/149‑1 (MQ=255)
                                          tGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAgg                                                                                                       <  1:3764995/149‑1 (MQ=255)
                                               aGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGaa                                                                                                  >  2:1229333/1‑149 (MQ=255)
                                                gCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAg                                                                                                 <  1:3455831/149‑1 (MQ=255)
                                                    aaaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                             <  1:3186778/149‑1 (MQ=255)
                                                    aaaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                             >  2:1476424/1‑149 (MQ=255)
                                                    aaaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                             <  1:3186945/149‑1 (MQ=255)
                                                    aaaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                             >  1:475079/1‑149 (MQ=255)
                                                      aaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGAt                                                                                           <  1:4239713/149‑1 (MQ=255)
                                                      aaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGAt                                                                                           >  2:3512313/1‑149 (MQ=255)
                                                      aaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGAt                                                                                           <  1:704115/149‑1 (MQ=255)
                                                       aGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATc                                                                                          <  1:1149331/149‑1 (MQ=255)
                                                        gCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCg                                                                                         >  2:411645/1‑149 (MQ=255)
                                                        gCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCg                                                                                         >  1:753741/1‑149 (MQ=255)
                                                         cGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGc                                                                                        <  1:1360033/149‑1 (MQ=255)
                                                              aGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCACAAGTCAGTGATTCTGTCCTGGCGAGCGAGAGTCAGGGCGAAGCAGGATCGCTGGTg                                                                                   >  2:2044134/1‑149 (MQ=255)
                                                               ggCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGa                                                                                  <  2:3057939/149‑1 (MQ=255)
                                                               ggCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGa                                                                                  >  2:2358043/1‑149 (MQ=255)
                                                                     cGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaaa                                                                            <  2:1625399/149‑1 (MQ=255)
                                                                      gCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaaaa                                                                           >  2:4406369/1‑149 (MQ=255)
                                                                       ccAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaaaaa                                                                          >  1:4632659/1‑149 (MQ=255)
                                                                         aGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACt                                                                        <  1:3119065/149‑1 (MQ=255)
                                                                          gTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTg                                                                       >  1:2919513/1‑149 (MQ=255)
                                                                          gTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTg                                                                       >  2:1707448/1‑149 (MQ=255)
                                                                            ttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTc                                                                     <  1:3914812/149‑1 (MQ=255)
                                                                              gCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTg                                                                   <  2:185894/149‑1 (MQ=255)
                                                                              gCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTg                                                                   <  1:1074733/149‑1 (MQ=255)
                                                                               cTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGa                                                                  <  2:2420037/149‑1 (MQ=255)
                                                                                tGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGGTGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGaa                                                                 <  2:1284614/149‑1 (MQ=255)
                                                                                tGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGaa                                                                 <  2:3675581/149‑1 (MQ=255)
                                                                                 ggCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGTTGATTGAAAAACTGTCTGAAt                                                                >  2:4614230/1‑149 (MQ=255)
                                                                                 ggCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAAt                                                                >  2:4163348/1‑149 (MQ=255)
                                                                                 ggCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAAt                                                                >  2:4614400/1‑149 (MQ=255)
                                                                                    taGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCATAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGcc                                                             <  2:679056/148‑1 (MQ=255)
                                                                                    aaGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGcc                                                             >  2:3232875/1‑149 (MQ=255)
                                                                                     aGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCa                                                            >  2:531008/1‑149 (MQ=255)
                                                                                     aGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCa                                                            <  1:3233865/149‑1 (MQ=255)
                                                                                      gcggcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAg                                                           >  2:2524723/1‑149 (MQ=255)
                                                                                        ggcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGaa                                                         <  2:1743468/149‑1 (MQ=255)
                                                                                          cgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAg                                                       >  1:1751175/1‑149 (MQ=255)
                                                                                           gTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgt                                                      >  2:3560063/1‑149 (MQ=255)
                                                                                            ttGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgtg                                                     <  2:3968996/149‑1 (MQ=255)
                                                                                              gACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGtgtg                                                   >  1:1283194/1‑149 (MQ=255)
                                                                                                cGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCt                                                 <  1:4563754/149‑1 (MQ=255)
                                                                                                   ccAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCg                                              <  1:1325564/149‑1 (MQ=255)
                                                                                                    cAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGt                                             >  1:2436601/1‑149 (MQ=255)
                                                                                                      gggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGc                                           <  1:2422371/149‑1 (MQ=255)
                                                                                                       ggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGcc                                          <  1:4406369/149‑1 (MQ=255)
                                                                                                           tggtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGt                                      <  2:4262408/149‑1 (MQ=255)
                                                                                                             gtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTgg                                    <  2:1494755/149‑1 (MQ=255)
                                                                                                               ggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCa                                  >  1:1547180/1‑149 (MQ=255)
                                                                                                                gtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAg                                 >  1:531951/1‑149 (MQ=255)
                                                                                                                   aTGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCa                              >  1:3978696/1‑149 (MQ=255)
                                                                                                                     ggATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCaaa                            >  2:1900634/1‑149 (MQ=255)
                                                                                                                        ttACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGc                         <  1:2065038/149‑1 (MQ=255)
                                                                                                                          aCGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAg                       >  2:3634081/1‑149 (MQ=255)
                                                                                                                          aCGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAg                       <  2:644075/149‑1 (MQ=255)
                                                                                                                            ggTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGa                     >  2:2137043/1‑149 (MQ=255)
                                                                                                                            ggTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGa                     <  2:4662829/149‑1 (MQ=255)
                                                                                                                              tGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAtc                   >  2:4217590/1‑149 (MQ=255)
                                                                                                                                aTTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAtctc                 <  1:2992457/149‑1 (MQ=255)
                                                                                                                                aTTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAtctc                 >  2:4747362/1‑149 (MQ=255)
                                                                                                                                       cGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGc          >  1:21875/1‑149 (MQ=255)
                                                                                                                                        gTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGcc         >  1:1970066/1‑149 (MQ=255)
                                                                                                                                        gTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGcc         >  2:3910445/1‑149 (MQ=255)
                                                                                                                                        gTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGcc         <  2:1517869/149‑1 (MQ=255)
                                                                                                                                         tCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCg        >  1:3315204/1‑149 (MQ=255)
                                                                                                                                               tGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCat  <  1:1848635/149‑1 (MQ=255)
                                                                                                                                               tGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCat  <  1:1847772/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCAT  >  CP009273/3231203‑3231494

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