Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,651,665 G→T 90.4% 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,6650GT90.4% 968.8 / 92.7 345G210G (GGC→GGAhscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base G (19/14);  new base T (166/145);  total (185/160)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 345 total aligned reads displayed.

GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATC  >  CP009273/2651518‑2651811
                                                                                                                                                   |                                                                                                                                                  
gCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTcc                                                                                                                                                  >  2:3033938/1‑150 (MQ=255)
 cTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCa                                                                                                                                                 >  1:1515866/1‑150 (MQ=255)
 cTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCa                                                                                                                                                 <  1:3566090/150‑1 (MQ=255)
 cTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCa                                                                                                                                                 <  2:3833227/150‑1 (MQ=255)
    cGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCg                                                                                                                                              >  1:3824354/1‑150 (MQ=255)
     gATCAGGATTGCCCGCCTGCTCGCGATTGTAATCCGCCAGCAGTTGGTCGAAATCATCGCCGCGGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGa                                                                                                                                             <  1:1559487/150‑1 (MQ=255)
     gATCAGGAATGCCCGCCTGCTCGCGATTGTAATCCGCCAGCAGTTGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCGGGTTGCCAGCAATTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGa                                                                                                                                             <  1:3746412/150‑1 (MQ=255)
         aGGAATGCCCGCCTGCTCGCTAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCAATTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTc                                                                                                                                         <  1:280188/150‑1 (MQ=255)
           gattGCCCGCCTGCTCGCTAATGTAATCCGCCAGCAGATGGTCGAAATCTTCGCCGCCGCGCGCGGATTCACCGCCGGTTGCCAGCACTTAAAACACGCCGAGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCAt                                                                                                                                       <  1:132181/147‑1 (MQ=255)
            aaTGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATa                                                                                                                                      <  2:3554464/150‑1 (MQ=255)
             aTGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATaa                                                                                                                                     <  2:4350106/150‑1 (MQ=255)
                cccGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACg                                                                                                                                  >  1:2310578/1‑150 (MQ=255)
                        cTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCAc                                                                                                                          >  1:4313181/1‑150 (MQ=255)
                        cTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCAc                                                                                                                          >  1:2124733/1‑150 (MQ=255)
                         tCGCGATTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCATGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACg                                                                                                                         <  1:4356133/150‑1 (MQ=255)
                         tCGCGAATGTAATCCGCCAGCAGATGGTCTAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACg                                                                                                                         <  2:268436/150‑1 (MQ=255)
                           gcgAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGcc                                                                                                                       >  2:92092/1‑150 (MQ=255)
                            cgAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCt                                                                                                                      >  2:2797972/1‑150 (MQ=255)
                              attttAATCCGCCAGATGTTGGTCGAAATCATCGCCGCCGAGCGCGTAATCACCGCCGGTTGCCAGCACTTAAAACCCGCCGCGACTTAAGCGCAGAATGTAAATATCAACCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                    <  1:4211181/146‑1 (MQ=255)
                              aaTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                    >  1:1789984/1‑150 (MQ=255)
                                tttAATCCGCCAGCTGTTGGTCTAAAACATCGCCGCCGAGCGCGTAATCACCGCGGGTTGCCAGAACTTCAAACACGCCGCGACTTAAGCGCAGATTGGAAATATCAACCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCt                                                                                                                  <  1:3005279/148‑1 (MQ=255)
                                tGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCt                                                                                                                  <  2:3255881/150‑1 (MQ=255)
                                tGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCt                                                                                                                  >  2:3197433/1‑150 (MQ=255)
                                tGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCt                                                                                                                  >  1:1650894/1‑150 (MQ=255)
                                 gTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTg                                                                                                                 >  1:4014978/1‑150 (MQ=255)
                                 gTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTg                                                                                                                 <  2:690391/150‑1 (MQ=255)
                                  tAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGa                                                                                                                <  2:1201894/150‑1 (MQ=255)
                                  tAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGa                                                                                                                >  2:4060438/1‑150 (MQ=255)
                                  tAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGa                                                                                                                <  2:2514093/150‑1 (MQ=255)
                                  tAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGa                                                                                                                >  1:1232376/1‑150 (MQ=255)
                                      ccGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCgg                                                                                                            >  1:583856/1‑150 (MQ=255)
                                      ccGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCgg                                                                                                            >  1:3872068/1‑150 (MQ=255)
                                       cGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAACTATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGa                                                                                                           >  2:4014726/1‑150 (MQ=255)
                                       cGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGa                                                                                                           <  2:3911699/150‑1 (MQ=255)
                                           agcagATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATcc                                                                                                       >  1:4033607/1‑150 (MQ=255)
                                                 tGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCg                                                                                                 >  1:786201/1‑150 (MQ=255)
                                                   gTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTa                                                                                               <  2:4181960/150‑1 (MQ=255)
                                                       aaaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACTTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTACCTGACCGGAATCCAGCCCGTAGGCg                                                                                           >  2:3706478/1‑150 (MQ=255)
                                                       aaaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCg                                                                                           >  2:3419551/1‑150 (MQ=255)
                                                        aaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGa                                                                                          >  1:2419122/1‑150 (MQ=255)
                                                            atcGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGc                                                                                      >  1:2347871/1‑150 (MQ=255)
                                                                     gAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggt                                                                             >  1:812729/1‑150 (MQ=255)
                                                                      aGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtc                                                                            <  2:651850/150‑1 (MQ=255)
                                                                         gcgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtcggt                                                                         <  2:3802678/150‑1 (MQ=255)
                                                                         gcgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtcggt                                                                         <  2:3355927/150‑1 (MQ=255)
                                                                          cgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGtt                                                                        >  1:4165503/1‑150 (MQ=255)
                                                                            gAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCg                                                                      >  2:2837622/1‑150 (MQ=255)
                                                                            gAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCg                                                                      >  2:1851973/1‑150 (MQ=255)
                                                                              aacACCGCCGGTTGCCAGAAATAAAAACACGCCGCTACTTAAGCGCATAATGTAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGtt                                                                    <  1:1481153/148‑1 (MQ=255)
                                                                              aTCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGtt                                                                    <  2:2892357/150‑1 (MQ=255)
                                                                               tCACCGCCGGTTGCCAGCACTTCAAACCCGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTta                                                                   <  1:3205839/150‑1 (MQ=255)
                                                                                 aCCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaag                                                                 <  1:1026230/150‑1 (MQ=255)
                                                                                 aCCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGGtaag                                                                 >  1:397065/1‑150 (MQ=255)
                                                                                    gccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagtaa                                                              <  2:1128151/150‑1 (MQ=255)
                                                                                     ccgGTTGCCTGCACTTCAAACACTCCGCGACTTAAGCGCAGAATGGATATATCAAACGTCCCTCGACCGTGTTCTTAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTttagtaat                                                             <  2:1975951/150‑2 (MQ=255)
                                                                                      cgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgc                                                            >  1:3300927/1‑150 (MQ=255)
                                                                                      cgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgc                                                            >  2:1607075/1‑150 (MQ=255)
                                                                                        gTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgcgc                                                          >  2:2922384/1‑150 (MQ=255)
                                                                                            ccAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGa                                                      <  1:3057389/150‑1 (MQ=255)
                                                                                             cAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGAc                                                     >  2:2865578/1‑150 (MQ=255)
                                                                                              aGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACg                                                    >  2:4174470/1‑150 (MQ=255)
                                                                                              aGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACg                                                    >  2:2469831/1‑150 (MQ=255)
                                                                                              aGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACg                                                    <  1:1955118/150‑1 (MQ=255)
                                                                                                cACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTg                                                  <  2:1810961/150‑1 (MQ=255)
                                                                                                  cTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGaa                                                >  1:2496753/1‑150 (MQ=255)
                                                                                                  cTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGaa                                                >  2:2022144/1‑150 (MQ=255)
                                                                                                   ttCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAg                                               >  1:4105729/1‑150 (MQ=255)
                                                                                                   ttCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAg                                               >  2:3133577/1‑150 (MQ=255)
                                                                                                   ttCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAg                                               <  2:551204/150‑1 (MQ=255)
                                                                                                     cAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGc                                             <  1:891726/150‑1 (MQ=255)
                                                                                                           cgccgcGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCa                                       <  2:2284169/150‑1 (MQ=255)
                                                                                                              cgcgACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAc                                    >  1:699596/1‑150 (MQ=255)
                                                                                                                  aCTTAAGCGCAGAATGGAAATATCAAACGTTCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACgcgc                                >  1:1449590/1‑150 (MQ=255)
                                                                                                                   cTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                               <  1:2607053/150‑1 (MQ=255)
                                                                                                                   cTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                               >  2:531012/1‑150 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                             >  1:2588345/1‑150 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                             >  1:2739672/1‑150 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                             >  1:3114591/1‑150 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                             >  2:1394297/1‑150 (MQ=255)
                                                                                                                      aaGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCcgcg                            >  1:2062276/1‑150 (MQ=255)
                                                                                                                       aGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGt                           >  2:2249886/1‑150 (MQ=255)
                                                                                                                        gcgcAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTc                          <  2:348996/150‑1 (MQ=255)
                                                                                                                        gcgcAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAACCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTc                          <  2:337771/150‑1 (MQ=255)
                                                                                                                            ataaTGTAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTg                      <  1:2927844/148‑1 (MQ=255)
                                                                                                                            ataaTGTAAATAACAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTg                      <  1:755263/148‑1 (MQ=255)
                                                                                                                            aGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTg                      <  1:1888318/150‑1 (MQ=255)
                                                                                                                             gAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTgg                     >  1:2696265/1‑150 (MQ=255)
                                                                                                                              aaTGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                    >  1:986868/1‑150 (MQ=255)
                                                                                                                              aaTGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                    <  2:1069449/150‑1 (MQ=255)
                                                                                                                                     atatCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa             >  1:2908948/1‑150 (MQ=255)
                                                                                                                                      tatCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTGAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGAc            >  2:625324/1‑150 (MQ=255)
                                                                                                                                       atCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACg           >  1:3520946/1‑150 (MQ=255)
                                                                                                                                       atCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACg           <  1:1804256/150‑1 (MQ=255)
                                                                                                                                        tCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCTATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGc          >  1:102709/1‑150 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg        >  1:1984361/1‑150 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTg        <  2:1525903/150‑1 (MQ=255)
                                                                                                                                           aaCGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg       >  2:742605/1‑150 (MQ=255)
                                                                                                                                            aCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgtg      <  2:2084253/150‑3 (MQ=255)
                                                                                                                                            aCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTggg      <  2:812729/150‑1 (MQ=255)
                                                                                                                                             cGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGc     >  2:2574292/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATC  >  CP009273/2651518‑2651811

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: