Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,953,192 C→A 15.9% D866Y (GAC→TAC)  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,1920CA15.9% 831.9 / 109.1 279D866Y (GAC→TAC) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base C (142/91);  new base A (15/30);  total (158/121)
Fisher's exact test for biased strand distribution p-value = 8.91e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

ATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTG  >  CP009273/2953044‑2953341
                                                                                                                                                    |                                                                                                                                                     
aTTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGCAGCGGCGGAACACGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCg                                                                                                                                                      <  2:1197829/150‑1 (MQ=255)
    cACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGa                                                                                                                                                  >  1:1565908/1‑150 (MQ=255)
       cAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGACTAAATGGCTCGGTGTCGGGGAttt                                                                                                                                               >  2:4071130/1‑150 (MQ=255)
           aGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCt                                                                                                                                           <  2:1153870/150‑1 (MQ=255)
            gcgcCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTg                                                                                                                                          <  1:2413844/150‑1 (MQ=255)
               ccATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCt                                                                                                                                       <  1:5072169/150‑1 (MQ=255)
                     ggTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTa                                                                                                                                 >  1:3379499/1‑150 (MQ=255)
                     ggTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTa                                                                                                                                 >  1:3654712/1‑150 (MQ=255)
                     ggTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTa                                                                                                                                 >  2:4206619/1‑150 (MQ=255)
                     ggTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGAGTTCGCTGTCTTCAGTa                                                                                                                                 >  1:2299612/1‑150 (MQ=255)
                      gTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTAc                                                                                                                                >  2:108150/1‑150 (MQ=255)
                       tAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACg                                                                                                                               <  2:3438962/150‑1 (MQ=255)
                        aaaTCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACgg                                                                                                                              >  1:4181535/1‑150 (MQ=255)
                          aTCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGaa                                                                                                                            <  2:3932493/150‑1 (MQ=255)
                           tCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAg                                                                                                                           >  1:5033437/1‑150 (MQ=255)
                                 gCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCAcc                                                                                                                     <  1:95402/150‑1 (MQ=255)
                                      tCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCaa                                                                                                                >  1:1193657/1‑150 (MQ=255)
                                       cGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCaaa                                                                                                               >  1:2902766/1‑150 (MQ=255)
                                                cggAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTgg                                                                                                      <  1:4026967/150‑1 (MQ=255)
                                                   aaCAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGaaa                                                                                                   <  1:2942033/150‑1 (MQ=255)
                                                    aCAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGaaaa                                                                                                  <  2:266945/150‑1 (MQ=255)
                                                      aaGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAAt                                                                                                <  2:4512110/150‑1 (MQ=255)
                                                      aaGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAAt                                                                                                <  1:539143/150‑1 (MQ=255)
                                                              ggCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCAc                                                                                        >  1:438642/1‑150 (MQ=255)
                                                               gCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCAcc                                                                                       >  1:2722905/1‑150 (MQ=255)
                                                               gCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCAcc                                                                                       >  2:4849682/1‑150 (MQ=255)
                                                                catcatCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCg                                                                                      >  2:1953745/1‑150 (MQ=255)
                                                                catcatCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCg                                                                                      >  2:1321440/1‑150 (MQ=255)
                                                                catcatCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCg                                                                                      >  2:57265/1‑150 (MQ=255)
                                                                catcatCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCg                                                                                      <  2:1672611/150‑1 (MQ=255)
                                                                 atcatcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCgg                                                                                     >  1:3493714/1‑150 (MQ=255)
                                                                 atcatcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCgg                                                                                     >  2:1388160/1‑150 (MQ=255)
                                                                 atcatcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCgg                                                                                     >  1:607055/1‑150 (MQ=255)
                                                                 atcatcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCgg                                                                                     <  1:4565271/150‑1 (MQ=255)
                                                                 atcatcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCgg                                                                                     <  1:2657113/150‑1 (MQ=255)
                                                                   catcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGAt                                                                                   >  2:3656195/1‑150 (MQ=255)
                                                                    atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGAtg                                                                                  >  2:1217127/1‑150 (MQ=255)
                                                                    atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGAtg                                                                                  >  2:3813526/1‑150 (MQ=255)
                                                                    atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGAtg                                                                                  <  1:4722307/150‑1 (MQ=255)
                                                                      ccTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGAtgtg                                                                                <  2:826594/150‑1 (MQ=255)
                                                                             aaCCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGaa                                                                         >  2:3988116/1‑150 (MQ=255)
                                                                              aCCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAAt                                                                        <  1:2534417/150‑1 (MQ=255)
                                                                              aCCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAAt                                                                        <  2:4034517/150‑1 (MQ=255)
                                                                               ccAGTGCATTCAATAACGGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATc                                                                       <  2:3207772/150‑1 (MQ=255)
                                                                                  gTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGtt                                                                    >  2:68480/1‑150 (MQ=255)
                                                                                     cATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTa                                                                 >  2:4127816/1‑150 (MQ=255)
                                                                                      aTTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAg                                                                >  1:1524841/1‑150 (MQ=255)
                                                                                      aTTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAg                                                                >  1:1157435/1‑150 (MQ=255)
                                                                                       ttCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGc                                                               <  1:2751181/150‑1 (MQ=255)
                                                                                         cAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGt                                                             <  2:1625735/150‑1 (MQ=255)
                                                                                          aaTAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGtt                                                            <  1:2810052/150‑1 (MQ=255)
                                                                                           aTAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGttt                                                           >  2:2124003/1‑150 (MQ=255)
                                                                                              aCTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCa                                                        >  1:1698996/1‑150 (MQ=255)
                                                                                              aCTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCa                                                        <  1:2772965/150‑1 (MQ=255)
                                                                                               ctgctgATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAg                                                       >  1:3091070/1‑150 (MQ=255)
                                                                                               ctgctgATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAg                                                       >  1:3935860/1‑150 (MQ=255)
                                                                                                 gctgATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCg                                                     <  2:5052636/150‑1 (MQ=255)
                                                                                                  ctgATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCgg                                                    >  2:489063/1‑150 (MQ=255)
                                                                                                     attgattTGATAACGGCTAAGTCCTTCGAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTATTCAGTACGGATGTTCACCTGCAATCGCATCTGGATAAATTCCGGAACCGGATGTGCCCATTATCGTTGTAGAGTTTACTGCGGCAc                                                 <  2:4326542/150‑1 (MQ=255)
                                                                                                     attgattTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCAc                                                 >  2:937275/1‑150 (MQ=255)
                                                                                                       tgatttgatAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACgg                                               >  1:1249852/1‑150 (MQ=255)
                                                                                                       tgatttgatAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACgg                                               >  2:517789/1‑150 (MQ=255)
                                                                                                        gatttgatTACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTAGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCGGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGt                                              <  2:2586098/150‑1 (MQ=255)
                                                                                                        gatgtgatAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGt                                              <  2:4367800/150‑1 (MQ=255)
                                                                                                            tgatAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTcc                                          <  2:252344/150‑1 (MQ=255)
                                                                                                               tAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggt                                       >  2:4017458/1‑150 (MQ=255)
                                                                                                               tAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggt                                       >  1:3032972/1‑150 (MQ=255)
                                                                                                               tAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggt                                       >  2:2993473/1‑150 (MQ=255)
                                                                                                                  cGGCTAAGTCCTTCCAGAATAACTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggtaag                                    <  2:2546079/150‑1 (MQ=255)
                                                                                                                  cGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggtaag                                    >  2:631463/1‑150 (MQ=255)
                                                                                                                   ggCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggtaagg                                   >  1:2755293/1‑150 (MQ=255)
                                                                                                                    gCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggtaaggt                                  <  2:4424285/150‑1 (MQ=255)
                                                                                                                     cTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggtaaggta                                 >  1:2064900/1‑150 (MQ=255)
                                                                                                                      tAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCggtaaggtaa                                >  1:1667446/1‑150 (MQ=255)
                                                                                                                         gTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACg                             >  2:366402/1‑150 (MQ=255)
                                                                                                                            cTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCa                          >  2:2924786/1‑150 (MQ=255)
                                                                                                                                   aaTAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTg                   >  2:4461221/1‑150 (MQ=255)
                                                                                                                                    aTAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGa                  <  2:2453848/150‑1 (MQ=255)
                                                                                                                                    aTAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGa                  <  1:3841127/150‑1 (MQ=255)
                                                                                                                                    aTAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGa                  <  1:3841635/150‑1 (MQ=255)
                                                                                                                                    aTAAATGGCTCGGTGTAGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGa                  <  1:4766753/150‑1 (MQ=255)
                                                                                                                                     tAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGaa                 >  2:6010/1‑150 (MQ=255)
                                                                                                                                     tAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGaa                 >  2:2977149/1‑150 (MQ=255)
                                                                                                                                     tAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGaa                 >  2:4806873/1‑150 (MQ=255)
                                                                                                                                     tAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGaa                 >  1:586028/1‑150 (MQ=255)
                                                                                                                                      aaaTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAAc                <  1:1421964/150‑1 (MQ=255)
                                                                                                                                      aaaTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAAc                <  1:2100166/150‑1 (MQ=255)
                                                                                                                                       aaTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACa               >  2:1836910/1‑150 (MQ=255)
                                                                                                                                       aaTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACa               >  1:2203493/1‑150 (MQ=255)
                                                                                                                                        aTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACaa              >  1:1068908/1‑150 (MQ=255)
                                                                                                                                        aTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACaa              >  1:4279381/1‑150 (MQ=255)
                                                                                                                                        aTGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACaa              >  1:3667443/1‑150 (MQ=255)
                                                                                                                                          ggCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGACCAAAt            >  1:1864483/1‑150 (MQ=255)
                                                                                                                                          ggCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAAt            <  2:4975244/150‑1 (MQ=255)
                                                                                                                                            cTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTc          >  2:1549246/1‑150 (MQ=255)
                                                                                                                                             tCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCa         >  2:2755687/1‑150 (MQ=255)
                                                                                                                                               ggTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGa       >  1:4944746/1‑150 (MQ=255)
                                                                                                                                               ggTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGa       >  2:3759465/1‑150 (MQ=255)
                                                                                                                                               ggTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGa       >  1:2472910/1‑150 (MQ=255)
                                                                                                                                                    cGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAAtgtg  <  1:1113411/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTTCACCAAAAGCGCCATACGGTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGGGATTTCGCTGTCTTCAGTACGGAAGTTCACCTGCAAACGCATCTGGAAAAATGCCCGCACCGGATGTGCCCAGAATCGTTGTAGCGTTTCCAGCGGCACGGTTTCCGGTAAGGTAAACGGCAGCGGCTGAACAAATTCAGAATGTG  >  CP009273/2953044‑2953341

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