Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,953,398 C→A 8.2% G797V (GGC→GTC)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,953,3980CA 8.2% 1094.1 / 54.0 305G797V (GGC→GTC) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base C (128/151);  new base A (9/17);  total (137/168)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

CCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCA  >  CP009273/2953249‑2953545
                                                                                                                                                     |                                                                                                                                                   
  cGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCct                                                                                                                                                   <  1:1526288/150‑1 (MQ=255)
  cGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCct                                                                                                                                                   <  1:3735738/150‑1 (MQ=255)
   gCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCctg                                                                                                                                                  >  2:4248192/1‑150 (MQ=255)
   gCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCctg                                                                                                                                                  >  2:1646797/1‑150 (MQ=255)
    cACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCctgg                                                                                                                                                 <  1:126678/150‑1 (MQ=255)
        ggATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCctggctgg                                                                                                                                             >  2:4465705/1‑150 (MQ=255)
        ggATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCctggctgg                                                                                                                                             >  2:1162836/1‑150 (MQ=255)
         gATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGt                                                                                                                                            >  2:303519/1‑150 (MQ=255)
          aTGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTa                                                                                                                                           <  1:2393304/150‑1 (MQ=255)
          aTGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTa                                                                                                                                           <  1:1869026/150‑1 (MQ=255)
            gtgCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACACATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGt                                                                                                                                         <  1:131003/150‑1 (MQ=255)
               cccAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGtttt                                                                                                                                      <  2:2547468/150‑1 (MQ=255)
               cccAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGtttt                                                                                                                                      >  1:2510347/1‑150 (MQ=255)
               cccAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGtttt                                                                                                                                      >  1:4001651/1‑150 (MQ=255)
                 cAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTtgt                                                                                                                                    <  1:179831/150‑1 (MQ=255)
                   gAATCGTTGTATCGTTTCAAGCTGCACGGTTTCAGGTAAGTTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACAAGCCTGGCTGGCCGCAGGAAGCCATTCACGAGCATAGCTTTGTCTTTCGACTGGCTGTAAGTTTtgtgt                                                                                                                                  >  2:3128682/1‑149 (MQ=255)
                    aaTCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGATCGCCTGGCTGGTAGTTTTGTGGa                                                                                                                                 >  2:4571063/1‑150 (MQ=255)
                    aaTCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCAATCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGa                                                                                                                                 >  2:4420999/1‑150 (MQ=255)
                     aTCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGAt                                                                                                                                <  1:1223210/150‑1 (MQ=255)
                        gTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCaa                                                                                                                             >  1:3333297/1‑150 (MQ=255)
                          tGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAAc                                                                                                                           >  2:3656795/1‑150 (MQ=255)
                          tGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAAc                                                                                                                           <  1:2659697/150‑1 (MQ=255)
                                    ccAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGAAGTTTTGTGGATCAAACGGCATCCggg                                                                                                                 >  2:4282122/1‑150 (MQ=255)
                                      aGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTa                                                                                                               >  2:1958436/1‑150 (MQ=255)
                                         ggCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATgg                                                                                                            >  2:2080433/1‑150 (MQ=255)
                                           cACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAg                                                                                                          <  1:804832/150‑1 (MQ=255)
                                                tttCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAg                                                                                                     >  2:1386611/1‑150 (MQ=255)
                                                           gtaaACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCt                                                                                          >  1:2910017/1‑150 (MQ=255)
                                                           gtaaACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCt                                                                                          >  1:274369/1‑150 (MQ=255)
                                                            taaACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCAGAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCtt                                                                                         <  2:4428072/150‑1 (MQ=255)
                                                             aaaCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCttt                                                                                        <  2:2653331/150‑1 (MQ=255)
                                                             aaaCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCttt                                                                                        >  2:131597/1‑150 (MQ=255)
                                                             aaaCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCttt                                                                                        <  1:5179128/150‑1 (MQ=255)
                                                             aaaCGGCAGCGGCTGAACAAAGTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCttt                                                                                        <  2:5246679/150‑1 (MQ=255)
                                                              aaCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCtttt                                                                                       <  1:1012120/150‑1 (MQ=255)
                                                               aCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTa                                                                                      <  1:2962154/150‑1 (MQ=255)
                                                               aCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTa                                                                                      <  1:905518/150‑1 (MQ=255)
                                                               aCGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTa                                                                                      <  1:5082300/150‑1 (MQ=255)
                                                                 ggCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTAcc                                                                                    >  2:4130535/1‑150 (MQ=255)
                                                                    aGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCtt                                                                                 <  1:2956564/150‑1 (MQ=255)
                                                                     gCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTg                                                                                >  2:4827282/1‑150 (MQ=255)
                                                                     gCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTg                                                                                <  2:3688495/150‑1 (MQ=255)
                                                                      cGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGc                                                                               >  1:2326596/1‑150 (MQ=255)
                                                                       ggCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGcc                                                                              >  1:1676465/1‑150 (MQ=255)
                                                                            aaCAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCt                                                                         >  1:4624945/1‑150 (MQ=255)
                                                                            aaCAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCt                                                                         >  1:2725822/1‑150 (MQ=255)
                                                                            aaCAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCt                                                                         >  1:439533/1‑150 (MQ=255)
                                                                             aCAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCtt                                                                        <  2:3537585/150‑1 (MQ=255)
                                                                               aaaTTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtc                                                                      >  2:1615171/1‑150 (MQ=255)
                                                                               aaaTTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtc                                                                      <  2:3997299/150‑1 (MQ=255)
                                                                               aaaTTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtc                                                                      <  1:1817599/150‑1 (MQ=255)
                                                                               aaaTTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGATTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTAtc                                                                      >  2:2376256/1‑150 (MQ=255)
                                                                                 aTTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcat                                                                    >  1:113548/1‑150 (MQ=255)
                                                                                 aTTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcat                                                                    <  2:1792423/150‑1 (MQ=255)
                                                                                   tCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcatca                                                                  <  2:1005030/150‑1 (MQ=255)
                                                                                   tCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcatca                                                                  <  2:2102676/150‑1 (MQ=255)
                                                                                   tCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcatca                                                                  <  1:4691884/150‑1 (MQ=255)
                                                                                   tCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTtcatca                                                                  <  2:4702800/150‑1 (MQ=255)
                                                                                    cAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTACGATGCGCTTTTACCCTTGCCTCGCTTTCATcac                                                                 >  2:626866/1‑150 (MQ=255)
                                                                                      gAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAg                                                               <  1:1596184/150‑1 (MQ=255)
                                                                                      gAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAg                                                               >  2:141938/1‑150 (MQ=255)
                                                                                            tcTTTACCAGCCTGGCTGGCAGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAgc                                                         <  2:2771512/149‑1 (MQ=255)
                                                                                                   cAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCAt                                                  >  1:627059/1‑150 (MQ=255)
                                                                                                     gCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCg                                                >  2:2623401/1‑150 (MQ=255)
                                                                                                     gCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCg                                                <  1:1157799/150‑1 (MQ=255)
                                                                                                     gCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCg                                                <  1:3189795/150‑1 (MQ=255)
                                                                                                       ctggctggCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGcc                                              >  1:4564810/1‑150 (MQ=255)
                                                                                                        tggctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGccc                                             <  1:3464642/150‑1 (MQ=255)
                                                                                                        tggctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGccc                                             >  1:4041712/1‑150 (MQ=255)
                                                                                                        tggctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGccc                                             <  1:3464550/150‑1 (MQ=255)
                                                                                                        tggctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTATTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGccc                                             <  2:3907348/150‑1 (MQ=255)
                                                                                                         ggctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCg                                            >  2:1126623/1‑150 (MQ=255)
                                                                                                          gctggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCgg                                           <  1:4433461/150‑1 (MQ=255)
                                                                                                            tggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTa                                         <  1:2851972/150‑1 (MQ=255)
                                                                                                            tggcCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTa                                         >  1:1630830/1‑150 (MQ=255)
                                                                                                               ccGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGAt                                      <  1:4978410/150‑1 (MQ=255)
                                                                                                                cGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCTCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATa                                     <  1:1540773/150‑1 (MQ=255)
                                                                                                                cGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATa                                     <  1:2181549/150‑1 (MQ=255)
                                                                                                                  cAGGTAGCCATTCACTAGCATAGCATTGTCTTTCGCCTGTCTGGTAGTTTTGTGGATCAAACGTCATCCGTTTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAAt                                   <  2:4907323/150‑1 (MQ=255)
                                                                                                                  cAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAAt                                   <  1:3065685/150‑1 (MQ=255)
                                                                                                                    ggTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGa                                 <  1:4938081/150‑1 (MQ=255)
                                                                                                                        gCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACttt                             >  1:4273693/1‑150 (MQ=255)
                                                                                                                        gCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACttt                             <  1:949336/150‑1 (MQ=255)
                                                                                                                         ccATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTg                            <  1:2535329/150‑1 (MQ=255)
                                                                                                                            ttCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGccc                         >  1:134835/1‑150 (MQ=255)
                                                                                                                             tCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCg                        <  1:4120194/150‑1 (MQ=255)
                                                                                                                               aCGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCTTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGAt                      <  2:2256170/150‑1 (MQ=255)
                                                                                                                                cGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATg                     <  2:3056027/150‑1 (MQ=255)
                                                                                                                                 gAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGt                    >  2:728224/1‑150 (MQ=255)
                                                                                                                                    cagaGCTTTGTCGTTCGCCTGGGTGGTAGTTTTGTGGATCAAACGGCATCCGGGTTTGGAGGCAAGTAAGTTGCGCTTTTACCCGTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCGCGGTAGTTAATGACTTTGCCCGATGTAGt                 <  1:4842803/147‑1 (MQ=255)
                                                                                                                                    cATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGt                 <  1:538769/150‑1 (MQ=255)
                                                                                                                                      tAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCg               >  1:2413711/1‑150 (MQ=255)
                                                                                                                                       aGCTTTGTCGTTCGACTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGa              <  1:3760373/150‑1 (MQ=255)
                                                                                                                                        gCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGAt             >  1:504541/1‑150 (MQ=255)
                                                                                                                                          tttGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCa           <  1:4809620/150‑1 (MQ=255)
                                                                                                                                            tGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGt         <  2:3387426/150‑1 (MQ=255)
                                                                                                                                            tGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGACGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTATCATCACAGTATAGCGATTCATCGCCCGTTAGATAATGAATATGCCCGATGTAGTCGATCAGt         >  1:3781146/1‑150 (MQ=255)
                                                                                                                                             gTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGtt        <  1:3509040/150‑1 (MQ=255)
                                                                                                                                                  tCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGc   >  2:673707/1‑150 (MQ=255)
                                                                                                                                                   cGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCa  >  2:2470446/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTGCTTTACCAGCCTGGCTGGCCGCAGGTAGCCATTCACGAGCATAGCTTTGTCGTTCGCCTGGCTGGTAGTTTTGTGGATCAAACGGCATCCGGGTATGGAGGCAAGTAAGATGCGCTTTTACCCTTGCCTCGCTTTCATCACAGTTGAGCGCTTCATCGCCCGGTAGATAATGACTTTGCCCGATGTAGTCGATCAGTTCCTGCA  >  CP009273/2953249‑2953545

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: