Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,945,306 C→G 100% G172A (GGC→GCC)  recD ← exonuclease V (RecBCD complex), alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,945,3060CG97.3% 720.4 / 8.1 183G172A (GGC→GCC) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base G (68/110);  total (71/112)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 199 total aligned reads displayed.

CAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGT  >  CP009273/2945161‑2945450
                                                                                                                                                 |                                                                                                                                                
cAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGggcg                                                                                                                                                >  2:1438111/1‑148 (MQ=255)
 aaaGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGggcgc                                                                                                                                               >  1:1270497/1‑148 (MQ=255)
   aGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGgcgccg                                                                                                                                             >  2:3156127/1‑148 (MQ=255)
   aGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGgcgccg                                                                                                                                             >  1:1721142/1‑148 (MQ=255)
   aGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCATTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGgcgccg                                                                                                                                             >  2:2388702/1‑148 (MQ=255)
    gCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCgg                                                                                                                                            >  2:2377537/1‑148 (MQ=255)
     ccTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGa                                                                                                                                           >  1:1285633/1‑148 (MQ=255)
       ttGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGaaa                                                                                                                                         >  2:3551890/1‑148 (MQ=255)
        tGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAAt                                                                                                                                        <  1:369131/148‑1 (MQ=255)
          ccGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCa                                                                                                                                      <  2:753527/148‑1 (MQ=255)
             agagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                   >  2:569601/1‑148 (MQ=255)
             agagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                   >  2:2762466/1‑148 (MQ=255)
              gagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                  >  2:327461/1‑148 (MQ=255)
              gagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                  >  2:3116580/1‑148 (MQ=255)
               agaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCgag                                                                                                                                 <  1:2440951/148‑1 (MQ=255)
                gaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCgaga                                                                                                                                >  1:3222892/1‑148 (MQ=255)
                     ggTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgc                                                                                                                           <  1:2199234/148‑1 (MQ=255)
                     ggTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgc                                                                                                                           <  2:2506593/148‑1 (MQ=255)
                       ttAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccgccg                                                                                                                         <  2:2008274/148‑1 (MQ=255)
                         aaGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCcgcg                                                                                                                       <  1:637127/148‑1 (MQ=255)
                         aaGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCcgcg                                                                                                                       <  1:1026985/148‑1 (MQ=255)
                           gcgcgcGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTc                                                                                                                     <  1:2112434/148‑1 (MQ=255)
                           gcgcgcGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTc                                                                                                                     >  2:67620/1‑148 (MQ=255)
                                  gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCa                                                                                                              >  1:1960582/1‑148 (MQ=255)
                                  gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCa                                                                                                              <  1:483900/148‑1 (MQ=255)
                                     gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTg                                                                                                           <  1:978678/148‑1 (MQ=255)
                                     gCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTg                                                                                                           >  1:2574480/1‑148 (MQ=255)
                                      cTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgc                                                                                                          >  1:1773388/1‑148 (MQ=255)
                                        ttACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgccg                                                                                                        >  2:3047641/1‑148 (MQ=255)
                                              gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCaa                                                                                                  <  1:1195004/148‑1 (MQ=255)
                                               ttGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAAc                                                                                                 >  2:1288985/1‑148 (MQ=255)
                                                tGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACt                                                                                                <  2:2120485/148‑1 (MQ=255)
                                                tGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACt                                                                                                >  1:3504430/1‑148 (MQ=255)
                                                  ttGCAGCCAGACTGATACGGCTGCTTTCTCAGTCGGCCATTTTAATTAACGCTGACAGCAAATTCGCTACTGTTTTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACAGAGTTCCGCCGCGTCAGCGCCACTGCCGCCGCAACttt                                                                                              <  1:1215925/147‑1 (MQ=255)
                                                  gTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACttt                                                                                              <  2:3395170/148‑1 (MQ=255)
                                                     cAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTg                                                                                           >  1:1516748/1‑148 (MQ=255)
                                                          aGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGt                                                                                      >  1:2094973/1‑148 (MQ=255)
                                                          aGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGt                                                                                      <  1:1632107/148‑1 (MQ=255)
                                                            aCGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTa                                                                                    <  1:1402553/148‑1 (MQ=255)
                                                            aCGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTa                                                                                    >  2:1594992/1‑148 (MQ=255)
                                                              ggATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAt                                                                                  <  1:684108/148‑1 (MQ=255)
                                                              ggATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAt                                                                                  <  1:836066/148‑1 (MQ=255)
                                                                aTACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAttt                                                                                >  2:20961/1‑148 (MQ=255)
                                                                  aCGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCa                                                                              >  2:656767/1‑148 (MQ=255)
                                                                  aCGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCa                                                                              >  1:3364329/1‑148 (MQ=255)
                                                                   cGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCAt                                                                             <  2:2799061/148‑1 (MQ=255)
                                                                     gCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCg                                                                           >  1:2306350/1‑148 (MQ=255)
                                                                      cAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGc                                                                          <  2:1324513/148‑1 (MQ=255)
                                                                       aGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCt                                                                         <  1:1403173/148‑1 (MQ=255)
                                                                         cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                       >  2:1198418/1‑148 (MQ=255)
                                                                          gTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTAc                                                                      <  2:622671/148‑1 (MQ=255)
                                                                             cGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTgg                                                                   <  2:2564698/148‑1 (MQ=255)
                                                                                  tCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGaaaaa                                                              <  1:1822864/148‑1 (MQ=255)
                                                                                  tCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGaaaaa                                                              <  2:1446619/148‑1 (MQ=255)
                                                                                    ggCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGt                                                            <  1:972447/148‑1 (MQ=255)
                                                                                     gCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGtt                                                           <  2:2250784/148‑1 (MQ=255)
                                                                                       cATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTg                                                         >  2:93306/1‑148 (MQ=255)
                                                                                        aTTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTGACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGt                                                        <  2:1879810/148‑1 (MQ=255)
                                                                                         tttGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTc                                                       <  2:3500286/148‑1 (MQ=255)
                                                                                         tttGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTc                                                       <  2:3219214/148‑1 (MQ=255)
                                                                                         tttGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTc                                                       <  2:1736992/148‑1 (MQ=255)
                                                                                            gAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAg                                                    <  2:2851462/148‑1 (MQ=255)
                                                                                            gAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAg                                                    <  1:3156127/148‑1 (MQ=255)
                                                                                              aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAggg                                                  <  1:1375309/148‑1 (MQ=255)
                                                                                              aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAggg                                                  <  2:1756247/148‑1 (MQ=255)
                                                                                               ttAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCa                                                     >  1:1027319/1‑144 (MQ=255)
                                                                                               ttAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCa                                                     >  1:2163985/1‑144 (MQ=255)
                                                                                               ttAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCa                                                     <  2:2163985/144‑1 (MQ=255)
                                                                                                     cTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTgcgc                                           <  2:1033355/148‑1 (MQ=255)
                                                                                                     cTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTgcgc                                           >  1:1446601/1‑148 (MQ=255)
                                                                                                        ccAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAg                                        >  1:1088838/1‑148 (MQ=255)
                                                                                                          aGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTa                                      >  2:1011387/1‑148 (MQ=255)
                                                                                                          aGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTa                                      >  2:1848019/1‑148 (MQ=255)
                                                                                                          aGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTa                                      >  2:2654020/1‑148 (MQ=255)
                                                                                                           gCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTag                                     <  1:2388702/148‑1 (MQ=255)
                                                                                                           gCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTag                                     <  1:3012198/148‑1 (MQ=255)
                                                                                                             aaCTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTagag                                   <  1:3361139/148‑1 (MQ=255)
                                                                                                              aCTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGc                                  >  1:3259280/1‑148 (MQ=255)
                                                                                                              aCTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGc                                  <  1:1982783/148‑1 (MQ=255)
                                                                                                              aCTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACATTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGc                                  >  1:3341360/1‑148 (MQ=255)
                                                                                                               cTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAATAGTTTGTCCAGGGTTTGCGCCAGTAGAGCt                                 >  2:521575/1‑148 (MQ=255)
                                                                                                               cTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCt                                 <  1:300556/148‑1 (MQ=255)
                                                                                                               cTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCt                                 >  1:469067/1‑148 (MQ=255)
                                                                                                                   gCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCGGGGTTTGCGCCAGTAGAGCTtcat                             <  2:1721705/148‑1 (MQ=255)
                                                                                                                   gCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcat                             <  2:1721141/148‑1 (MQ=255)
                                                                                                                    cTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcatc                            >  2:2939324/1‑148 (MQ=255)
                                                                                                                    cTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTtcatc                            >  2:1160824/1‑148 (MQ=255)
                                                                                                                      aCGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCaa                          <  1:2868811/148‑1 (MQ=255)
                                                                                                                      aCGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCaa                          <  1:3494968/148‑1 (MQ=255)
                                                                                                                        ggtggtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAAcc                        <  1:1581023/148‑1 (MQ=255)
                                                                                                                           ggtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCa                     <  1:3551890/148‑1 (MQ=255)
                                                                                                                              cGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATg                  <  2:3333606/148‑1 (MQ=255)
                                                                                                                               gTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATgg                 >  1:2008657/1‑148 (MQ=255)
                                                                                                                                    aCCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGa            >  1:2943756/1‑148 (MQ=255)
                                                                                                                                    aCCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGa            <  2:2334215/148‑1 (MQ=255)
                                                                                                                                    aCCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGa            <  1:740859/148‑1 (MQ=255)
                                                                                                                                     ccGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAt           <  2:1087587/148‑1 (MQ=255)
                                                                                                                                      cGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAtt          <  2:1968084/148‑1 (MQ=255)
                                                                                                                                        gTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCa        <  2:1670317/148‑1 (MQ=255)
                                                                                                                                            cAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTc    <  1:2615228/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCGT  >  CP009273/2945161‑2945450

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: