Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,434,587 G→T 30.8% S152* (TCG→TAG)  pfo ← pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,434,5870GT30.8% 622.1 / 247.8 308S152* (TCG→TAG) pfopyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (111/102);  new base T (62/33);  total (173/135)
Fisher's exact test for biased strand distribution p-value = 3.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

GGCGCGATGAGCATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAAC  >  CP009273/1434439‑1434736
                                                                                                                                                    |                                                                                                                                                     
   gcgATGAGCATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAAt                                                                                                                                                   >  1:601424/1‑150 (MQ=255)
    cgATGAGCATCAATTTCGACCTGCGGCATGAGATCAAGAAACGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATg                                                                                                                                                  >  2:3030347/1‑150 (MQ=255)
      aTGAGCATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATgag                                                                                                                                                <  2:805461/150‑1 (MQ=255)
      aTGAGCATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATgag                                                                                                                                                >  1:1229471/1‑150 (MQ=255)
       tGAGCATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATgaga                                                                                                                                               <  1:329900/150‑1 (MQ=255)
         aGCATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGc                                                                                                                                             >  2:533295/1‑150 (MQ=255)
           cATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCaa                                                                                                                                           <  1:224214/150‑1 (MQ=255)
           cATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCaa                                                                                                                                           >  1:1957848/1‑150 (MQ=255)
            aTCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCaaa                                                                                                                                          >  1:3494182/1‑150 (MQ=255)
             tCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAg                                                                                                                                         <  1:812463/150‑1 (MQ=255)
              cAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGt                                                                                                                                        <  1:3890239/150‑1 (MQ=255)
               aaTTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTc                                                                                                                                       >  1:4538967/1‑150 (MQ=255)
                 tttCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCtt                                                                                                                                     <  2:3485746/150‑1 (MQ=255)
                    cGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTgcg                                                                                                                                  <  1:4006682/150‑1 (MQ=255)
                    cGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTgcg                                                                                                                                  <  2:2047509/150‑1 (MQ=255)
                     gACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTgcgc                                                                                                                                 <  2:1717645/150‑1 (MQ=255)
                     gACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTgcgc                                                                                                                                 >  2:1015893/1‑150 (MQ=255)
                     gACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTgcgc                                                                                                                                 >  1:667936/1‑150 (MQ=255)
                     gACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAAAAAATTGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTgcgc                                                                                                                                 >  2:3381724/1‑150 (MQ=255)
                           cGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTg                                                                                                                           <  2:1525000/150‑1 (MQ=255)
                             gCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGa                                                                                                                         >  2:3745731/1‑150 (MQ=255)
                             gCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGAACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGa                                                                                                                         >  2:3958855/1‑150 (MQ=255)
                              cATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGAc                                                                                                                        >  2:1160957/1‑150 (MQ=255)
                              cATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGAc                                                                                                                        <  2:1959194/150‑1 (MQ=255)
                                  agaTCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGttt                                                                                                                    <  2:3578971/150‑1 (MQ=255)
                                   gaTCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTtg                                                                                                                   <  2:4247854/150‑1 (MQ=255)
                                    aTCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTtgc                                                                                                                  <  2:2872837/150‑1 (MQ=255)
                                    aTCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTtgc                                                                                                                  >  2:4020016/1‑150 (MQ=255)
                                      cAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTtgctg                                                                                                                >  2:4189727/1‑150 (MQ=255)
                                      cAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTtgctg                                                                                                                >  2:1135269/1‑150 (MQ=255)
                                        aGAATCGTGTCATCGGCCAGCGGGACAATTTTTTTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGca                                                                                                              <  1:3013150/150‑1 (MQ=255)
                                        aGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGca                                                                                                              >  1:1137454/1‑150 (MQ=255)
                                          aaTCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGcaca                                                                                                            >  1:3727655/1‑150 (MQ=255)
                                              gtgtCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACacaac                                                                                                        >  2:469122/1‑150 (MQ=255)
                                               tgtCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACacaaca                                                                                                       >  1:751699/1‑150 (MQ=255)
                                                 tCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATc                                                                                                     <  2:3929317/150‑1 (MQ=255)
                                                 tCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATc                                                                                                     <  2:1355597/150‑1 (MQ=255)
                                                   aTCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATcgc                                                                                                   >  1:1758817/1‑150 (MQ=255)
                                                   aTCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATcgc                                                                                                   >  1:4382470/1‑150 (MQ=255)
                                                     cGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCgcgc                                                                                                 >  2:4327644/1‑150 (MQ=255)
                                                      ggCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCa                                                                                                >  2:2599454/1‑150 (MQ=255)
                                                      ggCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATTGCGCa                                                                                                >  2:186875/1‑150 (MQ=255)
                                                       gCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCaa                                                                                               <  2:753638/150‑1 (MQ=255)
                                                              ggACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTATCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTAGCTGCACACAACAACGCGCAACCCGTCt                                                                                        >  1:3243048/1‑150 (MQ=255)
                                                                     tttATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCAc                                                                                 >  1:2200009/1‑150 (MQ=255)
                                                                       tATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCg                                                                               >  2:176899/1‑150 (MQ=255)
                                                                       tATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCg                                                                               >  2:3394876/1‑150 (MQ=255)
                                                                          tGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCa                                                                            >  1:1024041/1‑150 (MQ=255)
                                                                              ttCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAAc                                                                        <  1:1448286/150‑1 (MQ=255)
                                                                                cGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGt                                                                      >  1:2798718/1‑150 (MQ=255)
                                                                                 gTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTc                                                                     <  2:1654836/150‑1 (MQ=255)
                                                                                   gggACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCgg                                                                   >  1:4093858/1‑150 (MQ=255)
                                                                                   gggACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCgg                                                                   >  1:1763971/1‑150 (MQ=255)
                                                                                     gACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGaa                                                                 >  2:3143194/1‑150 (MQ=255)
                                                                                     gACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGaa                                                                 <  2:4625458/150‑1 (MQ=255)
                                                                                     gACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGaa                                                                 >  2:1421087/1‑150 (MQ=255)
                                                                                      aCGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAAt                                                                <  1:1412604/150‑1 (MQ=255)
                                                                                      aCGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAAt                                                                <  1:1237296/150‑1 (MQ=255)
                                                                                       cGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATg                                                               <  2:601424/150‑1 (MQ=255)
                                                                                            ggAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGc                                                          >  1:1091977/1‑150 (MQ=255)
                                                                                            ggAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGc                                                          >  1:2417918/1‑150 (MQ=255)
                                                                                             gAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGcc                                                         <  1:723445/150‑1 (MQ=255)
                                                                                             gAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGcc                                                         >  1:2892863/1‑150 (MQ=255)
                                                                                               aaCCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCaa                                                       >  2:2022055/1‑150 (MQ=255)
                                                                                                 ccATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCaaaa                                                     >  1:4296122/1‑150 (MQ=255)
                                                                                                  cATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCaaaaa                                                    >  2:676354/1‑150 (MQ=255)
                                                                                                   aTCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAAt                                                   >  1:360185/1‑150 (MQ=255)
                                                                                                   aTCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAAt                                                   <  2:3860476/150‑1 (MQ=255)
                                                                                                     cAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATag                                                 >  1:244926/1‑150 (MQ=255)
                                                                                                      aaaGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATaga                                                >  1:1500080/1‑150 (MQ=255)
                                                                                                      aaaGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATaga                                                <  1:3517049/150‑1 (MQ=255)
                                                                                                       aaGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATagag                                               >  2:3583459/1‑150 (MQ=255)
                                                                                                       aaGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATagag                                               >  1:527927/1‑150 (MQ=255)
                                                                                                         gAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATagagag                                             >  1:1191759/1‑150 (MQ=255)
                                                                                                          aaaTGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGt                                            <  2:2179825/150‑1 (MQ=255)
                                                                                                           aaTGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTg                                           >  2:4149251/1‑150 (MQ=255)
                                                                                                            aTGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGc                                          >  1:2970079/1‑150 (MQ=255)
                                                                                                            aTGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGc                                          >  1:4633253/1‑150 (MQ=255)
                                                                                                             tGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCa                                         <  2:3269996/150‑1 (MQ=255)
                                                                                                               aaTAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAtg                                       <  1:4364404/150‑1 (MQ=255)
                                                                                                                  aaaTGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCAtgtgt                                    >  1:2266355/1‑150 (MQ=255)
                                                                                                                    aTGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTgg                                  <  2:2744659/150‑1 (MQ=255)
                                                                                                                     tGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGc                                 >  1:2071124/1‑150 (MQ=255)
                                                                                                                     tGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGc                                 <  2:4011460/150‑1 (MQ=255)
                                                                                                                         aCGCGGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACg                             <  2:977211/150‑1 (MQ=255)
                                                                                                                           gcgGCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGt                           >  1:1722082/1‑150 (MQ=255)
                                                                                                                             ggCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTAc                         >  2:2549095/1‑150 (MQ=255)
                                                                                                                              gCTTTTCAGCGTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACg                        >  1:2645167/1‑150 (MQ=255)
                                                                                                                              gCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACg                        <  1:2378329/150‑1 (MQ=255)
                                                                                                                               cTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGt                       >  1:3504604/1‑150 (MQ=255)
                                                                                                                                     aGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCt                 >  1:3471372/1‑150 (MQ=255)
                                                                                                                                       cGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTAc               >  2:3298037/1‑150 (MQ=255)
                                                                                                                                        gTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACa              >  2:4265903/1‑150 (MQ=255)
                                                                                                                                        gTCGCGATTTGCTAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACa              >  2:3461778/1‑150 (MQ=255)
                                                                                                                                         tCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACAt             >  1:3342137/1‑150 (MQ=255)
                                                                                                                                           gcgATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGc           <  2:2058614/150‑1 (MQ=255)
                                                                                                                                             gATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAg         >  1:3424330/1‑150 (MQ=255)
                                                                                                                                               tttGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGa       <  2:3543074/150‑1 (MQ=255)
                                                                                                                                                ttGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGAc      <  1:3167022/150‑1 (MQ=255)
                                                                                                                                                    gAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAAc  <  2:1229672/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGCGATGAGCATCAATTTCGACCTGCGGCATGAGATCAAGAATCGTGTCATCGGCCAGCGGGACAATTTTATTGATTTCGTGGGACGTGCGGAAACCATCAAAGAAATGAATAAATGGCACGCGGCTTTTCAGCGTCGCGATTTGCGAAATGAGAGCAAAGTCTTGCGCTTCCTGGACGTTTGCTGCACACAACATCGCGCAACCCGTCTGGCGCACCGCCATAACGTCGGAATGATCGCCAAAAATAGAGAGTGCATGTGTGGCAACGGTACGTGCCGCTACATGCAGGACAAAC  >  CP009273/1434439‑1434736

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: