Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,651,665 G→T 100% 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,6650GT97.0% 1053.1 / 18.9 316G210G (GGC→GGAhscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base T (185/121);  total (191/125)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 316 total aligned reads displayed.

GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATC  >  CP009273/2651518‑2651811
                                                                                                                                                   |                                                                                                                                                  
gCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTcc                                                                                                                                                  <  2:1450225/150‑1 (MQ=255)
 cTACGATCAGGAATGCCCGCCTGCTCGCGATTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATAGCAAACGTCCCTCCa                                                                                                                                                 <  2:2024249/150‑1 (MQ=255)
 cTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCGTCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCa                                                                                                                                                 <  1:437740/150‑1 (MQ=255)
    cGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCg                                                                                                                                              <  2:4608992/150‑1 (MQ=255)
      aTCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAg                                                                                                                                            >  2:4080466/1‑150 (MQ=255)
        cAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGt                                                                                                                                          >  2:4234869/1‑150 (MQ=255)
         aGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTc                                                                                                                                         <  1:893393/150‑1 (MQ=255)
           gAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCAt                                                                                                                                       >  1:2720450/1‑150 (MQ=255)
           gAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCAt                                                                                                                                       >  1:1942665/1‑150 (MQ=255)
            aaTTCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATa                                                                                                                                      >  2:1442144/1‑150 (MQ=255)
            aaTGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCACACGTCCCTCCACCGAGGTCATa                                                                                                                                      <  1:574991/150‑1 (MQ=255)
              tGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATaaa                                                                                                                                    >  1:1110603/1‑150 (MQ=255)
               gCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAATGTCCCTCCACCGAGGTCATAAAc                                                                                                                                   >  1:250123/1‑150 (MQ=255)
               gCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAAc                                                                                                                                   <  2:3666535/150‑1 (MQ=255)
                 ccGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACgg                                                                                                                                 >  2:2311394/1‑150 (MQ=255)
                  cGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGc                                                                                                                                <  1:2204684/150‑1 (MQ=255)
                  cGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGc                                                                                                                                >  2:2366462/1‑150 (MQ=255)
                     cTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGAt                                                                                                                             >  2:4128833/1‑150 (MQ=255)
                      tGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATc                                                                                                                            >  2:1015261/1‑150 (MQ=255)
                         tCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGACTCACCGCCGGTTGCCAGCACTTCAAACACGCCTCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACg                                                                                                                         <  1:215235/150‑1 (MQ=255)
                               aTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGTACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTcc                                                                                                                   >  2:3077019/1‑150 (MQ=255)
                                  tAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGa                                                                                                                >  2:3456200/1‑150 (MQ=255)
                                   aaTCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGATGTCATAAACGGCGATCACGCCTTCCTGAc                                                                                                               <  2:1292061/150‑1 (MQ=255)
                                    aTCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACTGCGATCACGCCTTCCTGAcc                                                                                                              >  1:3767722/1‑150 (MQ=255)
                                    aTCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGAcc                                                                                                              >  2:2757898/1‑150 (MQ=255)
                                     tCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCg                                                                                                             >  1:3214500/1‑150 (MQ=255)
                                          cagcagATGTTCGAAATCATCGCCGCCGAGCGCTGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCAACGAGGTCATAAACGGCGATCACGCCTTCATGACCGGAATc                                                                                                        >  2:3966031/1‑150 (MQ=255)
                                          cagcagATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATc                                                                                                        >  1:1011229/1‑150 (MQ=255)
                                          cagcagATGGTCGAAATCATCGACGCCGAGCGCGGAATCACCGCCGGTTGCCAGCAATTCAAACACGCCGCGAATTAAGCGGTGAATGGAAATATCAAACGTCGCTCCACAGAGGTCATAAACGGAGATCACGCCTTGCTGACCGTAAtg                                                                                                        >  2:1403299/1‑149 (MQ=255)
                                             cagATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAg                                                                                                     <  1:4607770/150‑1 (MQ=255)
                                                   gTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTa                                                                                               >  2:2849600/1‑150 (MQ=255)
                                                   gTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTa                                                                                               >  2:1527507/1‑150 (MQ=255)
                                                     cGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAgg                                                                                             >  2:884230/1‑150 (MQ=255)
                                                      gAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGc                                                                                            >  2:1792559/1‑150 (MQ=255)
                                                      gAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGc                                                                                            >  2:2604005/1‑150 (MQ=255)
                                                       aaaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGTCTTCCTGACCGGAATCCAGCCCGTAGGCg                                                                                           >  2:1379252/1‑150 (MQ=255)
                                                       aaaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCg                                                                                           >  2:3106981/1‑150 (MQ=255)
                                                        aaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGa                                                                                          >  1:2122564/1‑150 (MQ=255)
                                                          tcatcGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATa                                                                                        <  1:3272448/150‑1 (MQ=255)
                                                          tcatcGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATa                                                                                        >  2:4347413/1‑150 (MQ=255)
                                                             tcGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGTCCGGAATCCAGCCCGTAGGCGATAGcc                                                                                     >  1:3560875/1‑150 (MQ=255)
                                                              cgccgcCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCg                                                                                    <  2:2407952/150‑1 (MQ=255)
                                                                 cgccgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAg                                                                                 <  2:1193905/150‑1 (MQ=255)
                                                                  gccgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGc                                                                                <  2:2758804/150‑1 (MQ=255)
                                                                    cgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcgg                                                                              >  1:1996415/1‑150 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                             >  2:2903985/1‑150 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                             >  1:3128413/1‑150 (MQ=255)
                                                                      aGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtc                                                                            >  2:2223310/1‑150 (MQ=255)
                                                                      aGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtc                                                                            <  1:4464030/150‑1 (MQ=255)
                                                                       gcgcGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtcg                                                                           <  1:3944475/150‑1 (MQ=255)
                                                                       gcgcGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtcg                                                                           <  1:3645389/150‑1 (MQ=255)
                                                                          cgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGtt                                                                        >  1:3803681/1‑150 (MQ=255)
                                                                          cgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGtt                                                                        >  1:1218239/1‑150 (MQ=255)
                                                                          cgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGtt                                                                        >  1:4033553/1‑150 (MQ=255)
                                                                            gAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCg                                                                      >  2:3490335/1‑150 (MQ=255)
                                                                              aTCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGtt                                                                    <  2:504383/150‑1 (MQ=255)
                                                                                     ccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagtaag                                                             <  1:408454/150‑1 (MQ=255)
                                                                                          tGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAg                                                        >  1:1549934/1‑150 (MQ=255)
                                                                                          tGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAg                                                        >  1:1809473/1‑150 (MQ=255)
                                                                                            ccAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGa                                                      <  1:3388976/150‑1 (MQ=255)
                                                                                              aGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACg                                                    <  1:3841548/150‑1 (MQ=255)
                                                                                              aGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACg                                                    >  2:3854591/1‑150 (MQ=255)
                                                                                                 aCTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGa                                                 <  2:3389803/150‑1 (MQ=255)
                                                                                                 aCTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGa                                                 >  2:3055115/1‑150 (MQ=255)
                                                                                                    tCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAgg                                              >  2:4170504/1‑150 (MQ=255)
                                                                                                      aaaCACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGcc                                            <  1:431103/150‑1 (MQ=255)
                                                                                                         cacGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGc                                         >  2:3609354/1‑150 (MQ=255)
                                                                                                             ccgcGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGa                                     <  2:4484428/150‑1 (MQ=255)
                                                                                                              cgcgACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAc                                    >  2:3225040/1‑150 (MQ=255)
                                                                                                              cgcgACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAc                                    >  1:1937051/1‑150 (MQ=255)
                                                                                                              cgcgACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAc                                    <  2:2518405/150‑1 (MQ=255)
                                                                                                                 gACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAcgcg                                 <  2:2434077/150‑1 (MQ=255)
                                                                                                                  aCTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACgcgc                                >  1:2309825/1‑150 (MQ=255)
                                                                                                                        gcgcAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTc                          <  2:3923335/150‑1 (MQ=255)
                                                                                                                           cAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCttt                       >  1:341340/1‑150 (MQ=255)
                                                                                                                            aGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTg                      >  2:1567147/1‑150 (MQ=255)
                                                                                                                             gAATGGAAATATCAAACGTCCCTCCACCGAGGTAATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTAGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTATTTgg                     <  2:1411177/150‑1 (MQ=255)
                                                                                                                              aaTGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                    >  1:642551/1‑150 (MQ=255)
                                                                                                                              aaTGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                    <  2:1429560/150‑1 (MQ=255)
                                                                                                                               aTGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTg                   >  1:26044/1‑150 (MQ=255)
                                                                                                                                 ggAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGcc                 <  2:827884/150‑1 (MQ=255)
                                                                                                                                  gAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGccc                >  1:1141301/1‑150 (MQ=255)
                                                                                                                                    aaTATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGTAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTg              <  1:956988/150‑1 (MQ=255)
                                                                                                                                     atatCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa             >  2:1815013/1‑150 (MQ=255)
                                                                                                                                      tatCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGAc            >  1:3023426/1‑150 (MQ=255)
                                                                                                                                        tCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGc          >  2:1756379/1‑150 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg        >  2:1032896/1‑150 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg        <  1:723987/150‑1 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg        <  1:2774151/150‑1 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg        >  2:3170695/1‑150 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGACGGTTCGTTAAGTAAGCGCAGGACGTCAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg        >  1:516139/1‑150 (MQ=255)
                                                                                                                                           aaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg       >  1:1411044/1‑150 (MQ=255)
                                                                                                                                           aaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg       >  1:4577156/1‑150 (MQ=255)
                                                                                                                                           aaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg       <  2:4567650/150‑1 (MQ=255)
                                                                                                                                           aaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg       <  2:843743/150‑1 (MQ=255)
                                                                                                                                             cGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc     <  2:1122772/150‑1 (MQ=255)
                                                                                                                                             cGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc     <  2:2278908/150‑1 (MQ=255)
                                                                                                                                             cGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc     <  1:3795413/150‑1 (MQ=255)
                                                                                                                                               ccccTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCAt   >  1:159345/2‑150 (MQ=255)
                                                                                                                                                cccTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATc  >  1:3816290/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATC  >  CP009273/2651518‑2651811

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: