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% 1171.3 / 12.4 309Q227P (CAG→CCG) ygjRputative NAD(P)‑dependent dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base C (132/165);  total (137/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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 312 total aligned reads displayed.

TGGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCAT  >  CP009273/3231202‑3231494
                                                                                                                                                   |                                                                                                                                                 
tGGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTgccgc                                                                                                                                                 <  1:2788510/150‑1 (MQ=255)
  gATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGCTGGATTACGGTGATTTCAGCGTCACCTTGCCGCAc                                                                                                                                               >  2:3144695/1‑150 (MQ=255)
  gATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCAc                                                                                                                                               >  2:1278/1‑150 (MQ=255)
  gATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCAc                                                                                                                                               >  2:2615/1‑150 (MQ=255)
       cGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaaa                                                                                                                                          <  2:4214712/150‑1 (MQ=255)
        ggTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaaaa                                                                                                                                         <  2:541405/150‑1 (MQ=255)
        ggTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCaaaa                                                                                                                                         >  1:3567461/1‑150 (MQ=255)
              tACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGt                                                                                                                                   >  2:2154233/1‑150 (MQ=255)
                  gTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGAtt                                                                                                                               <  2:4010001/150‑1 (MQ=255)
                  gTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACAGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGAtt                                                                                                                               >  1:4174136/1‑150 (MQ=255)
                   tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTc                                                                                                                              >  1:1689898/1‑150 (MQ=255)
                   tCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTc                                                                                                                              <  1:1850437/150‑1 (MQ=255)
                    cTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCt                                                                                                                             <  2:506746/150‑1 (MQ=255)
                        cGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTcc                                                                                                                         >  2:2082113/1‑150 (MQ=255)
                         gTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCt                                                                                                                        <  1:2348127/150‑1 (MQ=255)
                         gTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCt                                                                                                                        <  2:2319341/150‑1 (MQ=255)
                          tCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTg                                                                                                                       <  2:1105817/150‑1 (MQ=255)
                              cggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcga                                                                                                                   <  2:4532367/150‑1 (MQ=255)
                               ggTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcgag                                                                                                                  >  2:2970234/1‑150 (MQ=255)
                                 tGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGgcgagcg                                                                                                                <  1:1103918/150‑1 (MQ=255)
                                        tgttGGAGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                         <  1:1480029/148‑1 (MQ=255)
                                        aTTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCa                                                                                                         <  1:472286/150‑1 (MQ=255)
                                         tttGGTGAGCCGATAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAg                                                                                                        <  2:438569/150‑1 (MQ=255)
                                          ttGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAgg                                                                                                       <  1:4071837/150‑1 (MQ=255)
                                            ggTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGc                                                                                                     >  2:4941373/1‑150 (MQ=255)
                                                aGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAg                                                                                                 <  1:3855078/150‑1 (MQ=255)
                                                  ccGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCa                                                                                               <  2:2525106/150‑1 (MQ=255)
                                                  ccGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCa                                                                                               >  2:467108/1‑150 (MQ=255)
                                                    gAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAaggagg                                                                                             <  2:2482360/150‑1 (MQ=255)
                                                    gAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAaggagg                                                                                             <  2:463240/150‑1 (MQ=255)
                                                    gAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                             <  1:575657/150‑1 (MQ=255)
                                                    gAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAgg                                                                                             <  1:5065858/150‑1 (MQ=255)
                                                     aaaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGa                                                                                            >  2:2738597/1‑150 (MQ=255)
                                                      aaaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGAt                                                                                           <  2:158877/150‑1 (MQ=255)
                                                       aaGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATc                                                                                          <  1:757445/150‑1 (MQ=255)
                                                          cGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCt                                                                                       <  1:3953939/150‑1 (MQ=255)
                                                          cGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCt                                                                                       <  1:1139308/150‑1 (MQ=255)
                                                           gTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTg                                                                                      <  2:4926569/150‑1 (MQ=255)
                                                            tGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTgg                                                                                     >  1:2876734/1‑150 (MQ=255)
                                                             gCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGt                                                                                    <  1:4235381/150‑1 (MQ=255)
                                                              cAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTg                                                                                   <  2:352566/150‑1 (MQ=255)
                                                               aGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGa                                                                                  <  2:3957784/150‑1 (MQ=255)
                                                                ggCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGAt                                                                                 >  2:1876599/1‑150 (MQ=255)
                                                                  cAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTg                                                                               >  2:1774576/1‑150 (MQ=255)
                                                                   aaCCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGa                                                                              >  2:912833/1‑150 (MQ=255)
                                                                   aaCCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGa                                                                              >  2:3407354/1‑150 (MQ=255)
                                                                   aaCCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGa                                                                              >  1:1269665/1‑150 (MQ=255)
                                                                     ccGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaaa                                                                            <  2:2603796/150‑1 (MQ=255)
                                                                       gCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaaaaa                                                                          >  1:631697/1‑150 (MQ=255)
                                                                       gCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGaaaaa                                                                          >  2:5142632/1‑150 (MQ=255)
                                                                           gTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGt                                                                      <  1:3928054/150‑1 (MQ=255)
                                                                            tttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCTCTGGTGATTGAAAAACTGTc                                                                     >  1:914678/1‑150 (MQ=255)
                                                                            tttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTc                                                                     <  1:3334915/150‑1 (MQ=255)
                                                                             ttGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATGACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCt                                                                    <  2:3777574/150‑1 (MQ=255)
                                                                              tGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTg                                                                   <  1:446393/150‑1 (MQ=255)
                                                                               gCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGa                                                                  <  1:2082113/150‑1 (MQ=255)
                                                                                 tGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAAt                                                                >  1:2514621/1‑150 (MQ=255)
                                                                                   gCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGc                                                              <  2:660847/150‑1 (MQ=255)
                                                                                    cAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGcc                                                             <  1:3285355/150‑1 (MQ=255)
                                                                                     aaGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCa                                                            >  1:1915367/1‑150 (MQ=255)
                                                                                     aaGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCa                                                            >  2:3039811/1‑150 (MQ=255)
                                                                                      aGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAg                                                           <  1:268930/150‑1 (MQ=255)
                                                                                       gcggcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGa                                                          <  1:1649136/150‑1 (MQ=255)
                                                                                       gcggcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGa                                                          >  1:5013906/1‑150 (MQ=255)
                                                                                         ggcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGaaa                                                        <  2:4287643/150‑1 (MQ=255)
                                                                                         ggcgTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGaaa                                                        >  1:928262/1‑150 (MQ=255)
                                                                                            gTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgtg                                                     <  1:1386967/150‑1 (MQ=255)
                                                                                             ttGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAgtgt                                                    >  2:1908651/1‑150 (MQ=255)
                                                                                              tGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGtgtg                                                   >  2:3041451/1‑150 (MQ=255)
                                                                                               gACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGc                                                  <  2:3614650/150‑1 (MQ=255)
                                                                                                  gCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTc                                               <  1:2231331/150‑1 (MQ=255)
                                                                                                   cccAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCg                                              <  2:3329130/150‑1 (MQ=255)
                                                                                                    ccAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGt                                             >  2:2786502/1‑150 (MQ=255)
                                                                                                    ccAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGt                                             >  1:443289/1‑150 (MQ=255)
                                                                                                    ccAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGt                                             >  2:4128506/1‑150 (MQ=255)
                                                                                                        ggCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCcg                                         <  2:3063750/150‑1 (MQ=255)
                                                                                                            tggtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTg                                     <  1:1208945/150‑1 (MQ=255)
                                                                                                             ggtggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTgg                                    >  1:2453808/1‑150 (MQ=255)
                                                                                                                ggtgATGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAg                                 <  1:302306/150‑1 (MQ=255)
                                                                                                                     tGGATTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCaaa                            <  2:4722869/150‑1 (MQ=255)
                                                                                                                        aTTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGc                         >  2:3131637/1‑150 (MQ=255)
                                                                                                                        aTTACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGc                         <  2:3661884/150‑1 (MQ=255)
                                                                                                                         ttACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCa                        <  1:496167/150‑1 (MQ=255)
                                                                                                                          tACGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAg                       <  2:1431176/150‑1 (MQ=255)
                                                                                                                            cGGTGATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGa                     <  1:1774576/150‑1 (MQ=255)
                                                                                                                              gtgtTTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAtc                   <  2:5143120/150‑1 (MQ=255)
                                                                                                                                gATTTCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGAtctc                 <  2:2420800/150‑1 (MQ=255)
                                                                                                                                   ttCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGAGTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCAcc              <  2:749681/150‑1 (MQ=255)
                                                                                                                                    tCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCAccc             <  2:1661369/150‑1 (MQ=255)
                                                                                                                                    tCAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCAccc             <  1:5140745/150‑1 (MQ=255)
                                                                                                                                     cAGCGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCa            >  1:3777043/1‑150 (MQ=255)
                                                                                                                                        cGTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGcc         <  1:4609299/150‑1 (MQ=255)
                                                                                                                                         gTCACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCg        >  1:3090408/1‑150 (MQ=255)
                                                                                                                                           cACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgca      >  1:4743264/1‑150 (MQ=255)
                                                                                                                                           cACCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgca      >  1:2947422/1‑150 (MQ=255)
                                                                                                                                            aCCTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgcag     <  2:1347685/150‑1 (MQ=255)
                                                                                                                                             ccTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgcagc    <  1:3633094/150‑1 (MQ=255)
                                                                                                                                              cTTGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCgcagca   <  2:4816360/150‑1 (MQ=255)
                                                                                                                                               ttGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCat  >  1:1259142/1‑150 (MQ=255)
                                                                                                                                               ttGCCGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCat  <  2:959344/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGGATATCGGTTTTTACTGTCTGGCGTCGGCGGTGGCGTTATTTGGTGAGCCGAAAAGCGTGCAGGCAACCGCCAGTTTGCTGGCAAGCGGCGTTGACGCCCAGGGCGTGGTGGTGATGGATTACGGTGATTTCAGCGTCACCTTGCAGCACTCAAAAGTCAGTGATTCTGTCCTGGCGAGCGAGATTCAGGGCGAAGCAGGATCGCTGGTGATTGAAAAACTGTCTGAATGCCAGAAAGTGTGCTTCGTGCCGCGTGGCAGCCAAATGCAGGATCTCACCCAGCCGCAGCAT  >  CP009273/3231202‑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: ‑