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,3060CG96.9% 755.2 / 9.9 191G172A (GGC→GCC) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base G (89/96);  total (93/98)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 196 total aligned reads displayed.

CAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCG  >  CP009273/2945161‑2945449
                                                                                                                                                 |                                                                                                                                               
  aaGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGggcgcc                                                                                                                                             >  2:1117996/1‑148 (MQ=255)
      cTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGaa                                                                                                                                         <  1:3564572/148‑1 (MQ=255)
        tGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAAt                                                                                                                                       >  2:2659075/1‑148 (MQ=255)
        tGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAAt                                                                                                                                       >  2:883168/1‑148 (MQ=255)
          ccGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAACTCa                                                                                                                                     >  2:2764538/1‑148 (MQ=255)
          ccGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCa                                                                                                                                     <  1:1081727/148‑1 (MQ=255)
          ccGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCa                                                                                                                                     <  1:1081675/148‑1 (MQ=255)
             agagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                  >  2:2968774/1‑148 (MQ=255)
             agagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                  >  2:2249467/1‑148 (MQ=255)
             agagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCg                                                                                                                                  <  1:522975/148‑1 (MQ=255)
              gagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                 >  1:3699511/1‑148 (MQ=255)
              gagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                 <  2:203830/148‑1 (MQ=255)
              gagaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                 >  1:2007992/1‑148 (MQ=255)
              gagaTTAGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATAAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCga                                                                                                                                 >  1:1722447/1‑148 (MQ=255)
                gaTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCgaga                                                                                                                               <  2:3574343/148‑1 (MQ=255)
                 aTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGAt                                                                                                                              >  2:893205/1‑148 (MQ=255)
                 aTTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGAt                                                                                                                              >  2:2298644/1‑148 (MQ=255)
                   tCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATcc                                                                                                                            >  1:1777373/1‑148 (MQ=255)
                    cGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccg                                                                                                                           <  2:2610337/148‑1 (MQ=255)
                    cGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATccg                                                                                                                           <  1:2820463/148‑1 (MQ=255)
                        tAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCcgccgc                                                                                                                       >  1:3720586/1‑148 (MQ=255)
                        tAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCcgccgc                                                                                                                       >  1:610009/1‑148 (MQ=255)
                         aaGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCcgcg                                                                                                                      >  2:1796081/1‑148 (MQ=255)
                                cgGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAgcgc                                                                                                               <  1:787077/148‑1 (MQ=255)
                                  gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCa                                                                                                             <  2:47166/148‑1 (MQ=255)
                                  gCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCa                                                                                                             <  1:2659075/148‑1 (MQ=255)
                                       tttACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                        >  2:3545061/1‑148 (MQ=255)
                                       tttACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgcc                                                                                                        >  2:1649512/1‑148 (MQ=255)
                                        ttACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgccg                                                                                                       >  1:2518384/1‑148 (MQ=255)
                                        ttACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgccg                                                                                                       <  2:2355771/148‑1 (MQ=255)
                                        ttACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgccg                                                                                                       <  2:2356031/148‑1 (MQ=255)
                                        ttACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgccg                                                                                                       <  2:2356032/148‑1 (MQ=255)
                                          aCCCGTTGGTGCAGACAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTGTTCCCGGTGCCCGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTgccgcc                                                                                                     >  2:1973006/1‑148 (MQ=255)
                                           cccGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGccgccg                                                                                                    >  1:2312647/1‑148 (MQ=255)
                                              gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCGGTCGGCCATTTGAAGTACCGCTGCCAGCAACTTCGCTACGGGGGTCGTTTTACCGGTGCCAGGGGCGCCGGTATTCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCaa                                                                                                 <  2:2336102/148‑1 (MQ=255)
                                              gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCaa                                                                                                 <  2:330769/148‑1 (MQ=255)
                                              gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCaa                                                                                                 <  2:94256/148‑1 (MQ=255)
                                              gTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCaa                                                                                                 <  1:2102803/148‑1 (MQ=255)
                                                tGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACt                                                                                               <  1:1464711/148‑1 (MQ=255)
                                                tGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACt                                                                                               <  1:572250/148‑1 (MQ=255)
                                                      aGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGc                                                                                         <  1:2579508/148‑1 (MQ=255)
                                                        ccAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCa                                                                                       >  1:1196268/1‑148 (MQ=255)
                                                          aGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGt                                                                                     <  1:3359954/148‑1 (MQ=255)
                                                           gACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGtt                                                                                    >  2:248419/1‑148 (MQ=255)
                                                                aTACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAttt                                                                               >  1:2328055/1‑148 (MQ=255)
                                                                aTACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAAttt                                                                               >  1:2318958/1‑148 (MQ=255)
                                                                  aCGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCa                                                                             >  2:3201225/1‑148 (MQ=255)
                                                                   cGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCAt                                                                            <  1:673923/148‑1 (MQ=255)
                                                                   cGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCAt                                                                            <  1:748762/148‑1 (MQ=255)
                                                                     gCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCg                                                                          <  2:2514902/148‑1 (MQ=255)
                                                                      cAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGc                                                                         <  2:440280/148‑1 (MQ=255)
                                                                       aGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCt                                                                        <  1:688772/148‑1 (MQ=255)
                                                                        gCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCtt                                                                       <  1:418979/148‑1 (MQ=255)
                                                                         cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                      <  1:1871723/148‑1 (MQ=255)
                                                                         cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                      <  1:2124279/148‑1 (MQ=255)
                                                                         cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                      <  1:2225060/148‑1 (MQ=255)
                                                                         cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                      <  2:833037/148‑1 (MQ=255)
                                                                         cGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTa                                                                      <  2:1836837/148‑1 (MQ=255)
                                                                           ttCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACt                                                                    >  1:2836531/1‑148 (MQ=255)
                                                                            tCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTg                                                                   >  1:1317975/1‑148 (MQ=255)
                                                                                  tCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGaaaaa                                                             >  2:721640/1‑148 (MQ=255)
                                                                                  tCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGaaaaa                                                             >  2:1292624/1‑148 (MQ=255)
                                                                                   cGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGTCAGTTAATTTCATCGCTTACTGGAAAAAg                                                            >  1:2275848/1‑148 (MQ=255)
                                                                                    ggCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGt                                                           <  2:2007992/148‑1 (MQ=255)
                                                                                    ggCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGt                                                           <  2:1722447/148‑1 (MQ=255)
                                                                                       cATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTg                                                        >  1:235873/1‑148 (MQ=255)
                                                                                         tttGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTc                                                      >  1:633776/1‑148 (MQ=255)
                                                                                         tttGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTc                                                      >  1:389800/1‑148 (MQ=255)
                                                                                          ttGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTcc                                                     >  1:1546330/1‑148 (MQ=255)
                                                                                           tGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCa                                                    <  1:1805346/148‑1 (MQ=255)
                                                                                             aaTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAgg                                                  <  1:857463/148‑1 (MQ=255)
                                                                                              aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAggg                                                 <  1:2461929/148‑1 (MQ=255)
                                                                                              aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAggg                                                 >  1:3486054/1‑148 (MQ=255)
                                                                                              aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAggg                                                 <  1:342404/148‑1 (MQ=255)
                                                                                              aTTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTATTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAggg                                                 >  2:1640240/1‑148 (MQ=255)
                                                                                               ttAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTTCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGt                                                >  1:3563724/1‑148 (MQ=255)
                                                                                               ttAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGt                                                >  1:1623493/1‑148 (MQ=255)
                                                                                               ttAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGt                                                >  2:3590882/1‑148 (MQ=255)
                                                                                               ttAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGt                                                >  1:2174096/1‑148 (MQ=255)
                                                                                                tAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGtt                                               >  2:3671398/1‑148 (MQ=255)
                                                                                                tAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGtt                                               >  2:36135/1‑148 (MQ=255)
                                                                                                tAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGtt                                               >  1:291573/1‑148 (MQ=255)
                                                                                                tAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGtt                                               >  1:1349337/1‑148 (MQ=255)
                                                                                                  aCGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTg                                             <  1:379080/148‑1 (MQ=255)
                                                                                                    gCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTgcg                                           <  1:2726498/148‑1 (MQ=255)
                                                                                                     cTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTgcgc                                          >  1:314389/1‑148 (MQ=255)
                                                                                                       gCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCa                                        >  1:2453249/1‑148 (MQ=255)
                                                                                                       gCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCa                                        >  1:2452037/1‑148 (MQ=255)
                                                                                                       gCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCAACGCTTACTGGAAAAAGTTTGTCCAGGGTTTTCGCCa                                        >  1:886584/1‑148 (MQ=255)
                                                                                                         cAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGt                                      <  1:1643447/148‑1 (MQ=255)
                                                                                                            cAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTaga                                   >  2:2311886/1‑148 (MQ=255)
                                                                                                              aCTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGc                                 <  2:1054151/148‑1 (MQ=255)
                                                                                                              aCTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCCCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGc                                 <  1:3118158/148‑1 (MQ=255)
                                                                                                                           ggtCGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCa                    <  2:761370/148‑1 (MQ=255)
                                                                                                                              cGTTTTACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATg                 >  2:1607391/1‑148 (MQ=255)
                                                                                                                                 tttACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCa              <  2:2950376/148‑1 (MQ=255)
                                                                                                                                 tttACCGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCa              >  2:1469701/1‑148 (MQ=255)
                                                                                                                                      cGGTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGAtt         <  2:1546330/148‑1 (MQ=255)
                                                                                                                                       ggTGCCAGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTc        <  1:1818453/148‑1 (MQ=255)
                                                                                                                                             aGGGGCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCg  >  2:505993/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCAGCCAGACGGATACGGCAGCGTTCGCCGTCGGCCATTTGAATTAACGCTGCCAGCAACTTCGCTACGGTGGTCGTTTTACCGGTGCCAGGGCCGCCGGAAATCACCGAGATCCGCCGCGTCAGCGCCACTGCCGCCGCAACTTTTTGCCAGTTAATTTCATCGCTTACTGGAAAAAGTTTGTCCAGGGTTTGCGCCAGTAGAGCTTCATCAACCTCAATGGCATGATTCACTTCG  >  CP009273/2945161‑2945449

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