Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,651,665 G→T 75.2% G210G (GGC→GGA hscA ← DnaK‑like molecular chaperone specific for IscU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,651,6650GT75.2% 675.3 / inf 387G210G (GGC→GGAhscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base G (45/51);  new base T (185/105);  total (231/156)
Fisher's exact test for biased strand distribution p-value = 3.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 387 total aligned reads displayed.

GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGT  >  CP009273/2651518‑2651813
                                                                                                                                                   |                                                                                                                                                    
gCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTcc                                                                                                                                                    <  1:3572646/150‑1 (MQ=255)
 cTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCa                                                                                                                                                   >  1:3775497/1‑150 (MQ=255)
     gATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGa                                                                                                                                               >  1:677610/1‑150 (MQ=255)
        cAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGt                                                                                                                                            >  1:2669018/1‑150 (MQ=255)
         aGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTc                                                                                                                                           >  1:3966263/1‑150 (MQ=255)
          ggAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCa                                                                                                                                          >  1:1601134/1‑150 (MQ=255)
               gCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAAc                                                                                                                                     <  2:1801394/150‑1 (MQ=255)
               gCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAAc                                                                                                                                     >  1:3179641/1‑150 (MQ=255)
               gCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAAc                                                                                                                                     >  1:1503599/1‑150 (MQ=255)
                 ccGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACgg                                                                                                                                   >  1:454106/1‑150 (MQ=255)
                  cGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGc                                                                                                                                  >  1:2301483/1‑150 (MQ=255)
                  cGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGc                                                                                                                                  <  2:3534106/150‑1 (MQ=255)
                   gCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGACATATCAATCGTCCCGCCACCGAGGTCATAAACGGCg                                                                                                                                 >  2:5033046/1‑150 (MQ=255)
                     cTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGAt                                                                                                                               >  2:3076812/1‑150 (MQ=255)
                     cTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGAt                                                                                                                               >  1:1642452/1‑150 (MQ=255)
                      tGCTCGCGAATGTAATCCGCCAGCAGAAGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATc                                                                                                                              >  2:1672442/1‑150 (MQ=255)
                       gCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCa                                                                                                                             <  2:5526804/150‑1 (MQ=255)
                        cTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCAc                                                                                                                            >  2:2389323/1‑150 (MQ=255)
                        cTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCAc                                                                                                                            >  2:5163002/1‑150 (MQ=255)
                        cTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCAc                                                                                                                            <  1:4681454/150‑1 (MQ=255)
                          cgcgAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGc                                                                                                                          >  2:5476807/1‑150 (MQ=255)
                              aaTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                      >  2:1447747/1‑150 (MQ=255)
                                       cGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGa                                                                                                             >  1:4166402/1‑150 (MQ=255)
                                         ccAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAAt                                                                                                           >  1:2218006/1‑150 (MQ=255)
                                            gcagATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCa                                                                                                        <  2:1166372/150‑1 (MQ=255)
                                            gcagATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCa                                                                                                        <  1:2067049/150‑1 (MQ=255)
                                                aTGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGccc                                                                                                    >  1:1701829/1‑150 (MQ=255)
                                                aTGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGccc                                                                                                    >  1:3486200/1‑150 (MQ=255)
                                                aTGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGccc                                                                                                    >  1:1133391/1‑150 (MQ=255)
                                                  ggTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGt                                                                                                  >  2:2488817/1‑150 (MQ=255)
                                                   gTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTa                                                                                                 >  1:3621269/1‑150 (MQ=255)
                                                      gAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGc                                                                                              <  1:2958041/150‑1 (MQ=255)
                                                      gAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGc                                                                                              >  1:1028024/1‑150 (MQ=255)
                                                       aaaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCg                                                                                             >  1:1996904/1‑150 (MQ=255)
                                                            atcGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGc                                                                                        <  1:2737459/150‑1 (MQ=255)
                                                                 cgccgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAg                                                                                   >  2:1237083/1‑150 (MQ=255)
                                                                  gccgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGc                                                                                  <  1:1397653/150‑1 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                               >  1:259269/1‑150 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                               >  1:931808/1‑150 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                               >  1:5417707/1‑150 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                               >  1:4916261/1‑150 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                               <  1:4023106/150‑1 (MQ=255)
                                                                      aGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtc                                                                              <  1:2803932/150‑1 (MQ=255)
                                                                      aGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtc                                                                              <  1:4428853/150‑1 (MQ=255)
                                                                      aGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtc                                                                              <  1:5076486/150‑1 (MQ=255)
                                                                       gcgcGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtcg                                                                             <  2:597883/150‑1 (MQ=255)
                                                                           ggAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTc                                                                         >  1:5047079/1‑150 (MQ=255)
                                                                              aTCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGAAAGCCGCAGCGGTCGCTTCGtt                                                                      >  2:3766710/1‑150 (MQ=255)
                                                                                 aCCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaag                                                                   <  2:92248/150‑1 (MQ=255)
                                                                                  ccgccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagt                                                                  <  2:3662288/150‑1 (MQ=255)
                                                                                  ccgccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagt                                                                  <  2:4688836/150‑1 (MQ=255)
                                                                                     ccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagtaag                                                               >  1:4982054/1‑150 (MQ=255)
                                                                                     ccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagtaag                                                               <  2:3336094/150‑1 (MQ=255)
                                                                                      cgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgc                                                              >  2:218287/1‑150 (MQ=255)
                                                                                       ggTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgcg                                                             >  1:4457007/1‑150 (MQ=255)
                                                                                         ttGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCa                                                           <  1:1486027/150‑1 (MQ=255)
                                                                                          tGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAg                                                          >  2:3389191/1‑150 (MQ=255)
                                                                                           gCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAgg                                                         <  1:5214925/150‑1 (MQ=255)
                                                                                            ccAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAGCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGa                                                        >  1:5139136/1‑150 (MQ=255)
                                                                                            ccAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGAAGGCGAAAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGa                                                        >  1:3987531/1‑150 (MQ=255)
                                                                                             cAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGAc                                                       <  2:4371109/150‑1 (MQ=255)
                                                                                              aGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACg                                                      <  2:4600532/150‑1 (MQ=255)
                                                                                                cACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTg                                                    <  1:1494690/150‑1 (MQ=255)
                                                                                                cACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTg                                                    <  2:2301473/150‑1 (MQ=255)
                                                                                                cACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTg                                                    >  2:2646513/1‑150 (MQ=255)
                                                                                                  cTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGaa                                                  >  2:40219/1‑150 (MQ=255)
                                                                                                   ttCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAg                                                 <  1:4632837/150‑1 (MQ=255)
                                                                                                    tCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAgg                                                >  1:1328829/1‑150 (MQ=255)
                                                                                                     cAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGc                                               >  1:3307960/1‑150 (MQ=255)
                                                                                                      aaaCACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGcc                                              <  2:3179641/150‑1 (MQ=255)
                                                                                                       aaCACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGACTCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGccc                                             >  2:4216473/1‑150 (MQ=255)
                                                                                                       aaCACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGccc                                             >  2:1768151/1‑150 (MQ=255)
                                                                                                        acacGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCg                                            <  2:631032/150‑1 (MQ=255)
                                                                                                         cacGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGc                                           >  1:1015907/1‑150 (MQ=255)
                                                                                                         cacGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGc                                           >  1:2999602/1‑150 (MQ=255)
                                                                                                           cgactctaCTTATGCGCAGAATGGATATTTCAAACGTCCCGCCTACGTGGTCTTAAACGGCGATCACGCCTTCCTGACCGTAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTATGTAAGCGAAGGACGTGAAGGCCCGCCa                                         <  1:2415777/143‑1 (MQ=255)
                                                                                                            gccgcGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAg                                        >  1:5597351/1‑150 (MQ=255)
                                                                                                              cgcgACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAc                                      >  1:1898711/1‑150 (MQ=255)
                                                                                                              cgcgACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAc                                      >  1:306059/1‑150 (MQ=255)
                                                                                                                 gACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAcgcg                                   >  2:488199/1‑150 (MQ=255)
                                                                                                                    ttAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGATAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccg                                <  2:694655/150‑1 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                               >  2:5383126/1‑150 (MQ=255)
                                                                                                                      aaGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCcgcg                              >  2:50058/1‑150 (MQ=255)
                                                                                                                          gcAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCtt                          >  2:1748016/1‑150 (MQ=255)
                                                                                                                             gAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTgg                       >  1:2120012/1‑150 (MQ=255)
                                                                                                                                tGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCGCGCCAGACGCGCCGCGTCTTTGGTGc                    >  2:3415454/1‑150 (MQ=255)
                                                                                                                                 ggAAATATCAAACGTCCCGCCACCGAGGTCATACACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGcc                   <  2:878316/150‑1 (MQ=255)
                                                                                                                                  gAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                  <  2:2835928/150‑1 (MQ=255)
                                                                                                                                    aaTATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCATCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACTCGCCGCCTCTATGGTGCCCTg                >  2:1861434/1‑150 (MQ=255)
                                                                                                                                    aaTATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTg                <  1:1534543/150‑1 (MQ=255)
                                                                                                                                      tatCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGAc              <  2:2154013/150‑1 (MQ=255)
                                                                                                                                       atCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACg             >  1:4250859/1‑150 (MQ=255)
                                                                                                                                         cAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCt           <  2:4516666/150‑1 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg          >  1:5470380/1‑150 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg          <  2:3046802/150‑1 (MQ=255)
                                                                                                                                           aaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg         >  1:2143732/1‑150 (MQ=255)
                                                                                                                                             cGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc       <  2:2052676/150‑1 (MQ=255)
                                                                                                                                              gTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCa      >  1:4151257/1‑150 (MQ=255)
                                                                                                                                              gTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCa      >  1:1243017/1‑150 (MQ=255)
                                                                                                                                              gTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCa      <  1:2206390/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGT  >  CP009273/2651518‑2651813

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: