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,6650GT96.2% 903.4 / 21.6 265G210G (GGC→GGAhscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base G (2/8);  new base T (136/118);  total (139/126)
Fisher's exact test for biased strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 266 total aligned reads displayed.

GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTC  >  CP009273/2651518‑2651814
                                                                                                                                                   |                                                                                                                                                     
 cTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCa                                                                                                                                                    >  1:1568253/1‑150 (MQ=255)
   actaTCAGTATTGCCCGCCTGCTCGAGATTTTAATCCGCCAGCATTTGTTCTAAATCTTCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGAAATTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCAcc                                                                                                                                                  <  1:4806508/147‑1 (MQ=255)
   aCGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCAcc                                                                                                                                                  <  1:3633085/150‑1 (MQ=255)
   aCGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCAcc                                                                                                                                                  <  2:3245389/150‑1 (MQ=255)
   aCGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATAATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCAcc                                                                                                                                                  >  1:1314891/1‑150 (MQ=255)
      aTCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGACTCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAg                                                                                                                                               >  2:268181/1‑150 (MQ=255)
      aTCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAg                                                                                                                                               >  2:268183/1‑150 (MQ=255)
          ggAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCa                                                                                                                                           <  1:4333896/150‑1 (MQ=255)
               gCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAAc                                                                                                                                      >  1:256241/1‑150 (MQ=255)
                    ccTGCTCGCTAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTAAAACACCCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGa                                                                                                                                 <  1:4057114/150‑1 (MQ=255)
                     gtGCTCGCGAATGTAATCCGCCATCTGTTGGTCTAAATCTTGGCCGCCGAGCGCGGAATCACCGCGGGTTGCCAGCACTTCAAACCCGCCGCGACTTATGCGAATATTGTAAATATAAAACGTCCCTCCACCGTGGTAATAAACGGCGAt                                                                                                                                <  2:4378656/149‑1 (MQ=255)
                      tGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGCGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATc                                                                                                                               <  2:949865/150‑1 (MQ=255)
                      tGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATc                                                                                                                               <  2:3785205/150‑1 (MQ=255)
                      tGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATc                                                                                                                               <  1:4528066/150‑1 (MQ=255)
                        cTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCAc                                                                                                                             <  1:3746480/150‑1 (MQ=255)
                         tCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACg                                                                                                                            >  1:2818991/1‑150 (MQ=255)
                         tCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACg                                                                                                                            <  2:1418761/150‑1 (MQ=255)
                          cgcgAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGc                                                                                                                           >  2:1560350/1‑150 (MQ=255)
                            cgAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCt                                                                                                                         <  2:4037969/150‑1 (MQ=255)
                             gAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCtt                                                                                                                        >  1:4241640/1‑150 (MQ=255)
                              aaTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCATAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                       <  1:1142107/150‑1 (MQ=255)
                              aaTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                       >  1:782072/1‑150 (MQ=255)
                              aaTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                       >  2:66871/1‑150 (MQ=255)
                              aaTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                       >  1:1494372/1‑150 (MQ=255)
                              aaTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTc                                                                                                                       <  1:3267028/150‑1 (MQ=255)
                               aTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTcc                                                                                                                      >  1:3268825/1‑150 (MQ=255)
                               aTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTcc                                                                                                                      >  1:2680011/1‑150 (MQ=255)
                               aTGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTcc                                                                                                                      <  2:1288648/150‑1 (MQ=255)
                                tGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCt                                                                                                                     <  2:182021/150‑1 (MQ=255)
                                  tAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGa                                                                                                                   >  2:1079362/1‑150 (MQ=255)
                                   aaTCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGAc                                                                                                                  >  1:1948755/1‑150 (MQ=255)
                                       cGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGa                                                                                                              <  2:1991513/150‑1 (MQ=255)
                                       cGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGa                                                                                                              <  2:4579535/150‑1 (MQ=255)
                                           agcagAGGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATcc                                                                                                          <  1:3218731/150‑1 (MQ=255)
                                              agATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGc                                                                                                       >  1:3787058/1‑150 (MQ=255)
                                                aTGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGccc                                                                                                     <  2:174497/150‑1 (MQ=255)
                                                      gAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGc                                                                                               >  1:3657684/1‑150 (MQ=255)
                                                      gAAATCAGCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGc                                                                                               <  1:1934166/150‑1 (MQ=255)
                                                       aaaTCATCGCCGCCGCGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCg                                                                                              <  2:4860479/150‑1 (MQ=255)
                                                        aaTCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGa                                                                                             >  1:413679/1‑150 (MQ=255)
                                                           catcGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAg                                                                                          >  2:1066861/1‑150 (MQ=255)
                                                             tcGCCGCCGCGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACCCGCCGCGACTTAAGCGCAGAATGGAAATATCACACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGcc                                                                                        <  1:1354877/150‑1 (MQ=255)
                                                              cgccgcCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCg                                                                                       <  2:3017504/150‑1 (MQ=255)
                                                              cgccgcCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCg                                                                                       <  2:766349/150‑1 (MQ=255)
                                                               gccgccGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGc                                                                                      >  2:3025290/1‑150 (MQ=255)
                                                                ccgccgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCa                                                                                     >  1:3494331/1‑150 (MQ=255)
                                                                  gccgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGc                                                                                   >  2:557325/1‑150 (MQ=255)
                                                                    cgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcgg                                                                                 >  1:312371/1‑150 (MQ=255)
                                                                    cgAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcgg                                                                                 >  1:2925547/1‑150 (MQ=255)
                                                                         gcgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAACCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGcggtcggt                                                                            <  1:4436695/150‑1 (MQ=255)
                                                                          cgGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGtt                                                                           >  1:3631665/1‑150 (MQ=255)
                                                                             aaTCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGt                                                                        <  2:2659085/150‑1 (MQ=255)
                                                                              aTCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAACCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGtt                                                                       <  1:2095871/150‑1 (MQ=255)
                                                                                cACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaa                                                                     >  1:2225096/1‑150 (MQ=255)
                                                                                 aCCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAACCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaag                                                                    <  2:3607828/150‑1 (MQ=255)
                                                                                 aCCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaag                                                                    <  2:3885408/150‑1 (MQ=255)
                                                                                   cgccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagta                                                                  <  1:4416459/150‑1 (MQ=255)
                                                                                   cgccgGTTGCCAGCACTCCAAACCCGCCGCGACTTAAGCTCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagta                                                                  <  1:1000575/150‑1 (MQ=255)
                                                                                    gccgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTtaagtaa                                                                 <  2:2290084/150‑1 (MQ=255)
                                                                                      cgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgc                                                               >  2:3450681/1‑150 (MQ=255)
                                                                                      cgGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgc                                                               >  1:2534566/1‑150 (MQ=255)
                                                                                       ggTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAgcg                                                              >  2:4469624/1‑150 (MQ=255)
                                                                                         ttGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCa                                                            >  2:5129722/1‑150 (MQ=255)
                                                                                           gCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAgg                                                          >  1:1720506/1‑150 (MQ=255)
                                                                                                  cTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGaa                                                   >  2:370271/1‑150 (MQ=255)
                                                                                                  cTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGaa                                                   >  1:5123810/1‑150 (MQ=255)
                                                                                                   ttCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAg                                                  <  2:4744611/150‑1 (MQ=255)
                                                                                                      aaaCACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGcc                                               >  1:1293039/1‑150 (MQ=255)
                                                                                                      aaaCACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGcc                                               >  1:4948546/1‑150 (MQ=255)
                                                                                                       aaCACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGccc                                              <  1:4701981/150‑1 (MQ=255)
                                                                                                          acGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGcc                                           >  2:2855757/1‑150 (MQ=255)
                                                                                                          acGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGcc                                           >  1:636531/1‑150 (MQ=255)
                                                                                                            gccgcGACTTAAGCGCAGAATGGAAATAACAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAg                                         <  1:1824384/150‑1 (MQ=255)
                                                                                                             ccgcGACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGa                                        >  1:1381478/1‑150 (MQ=255)
                                                                                                                 gACTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGAcgcg                                    >  1:4924576/1‑150 (MQ=255)
                                                                                                                   cTTACGCGCAGAATGGAAATATCACACGTCCCTCCACCGAGGTCATAACCGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGAGAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                                  <  1:3969004/150‑1 (MQ=255)
                                                                                                                   cTTAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgcc                                  >  2:5075359/1‑150 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                                >  2:4768667/1‑150 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                                >  2:3840674/1‑150 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                                <  2:4531308/150‑1 (MQ=255)
                                                                                                                     tAAGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGcgccgc                                >  2:1953889/1‑150 (MQ=255)
                                                                                                                      aaGCGCAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCcgcg                               >  2:1725736/1‑150 (MQ=255)
                                                                                                                         cgcAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCt                            <  1:2051372/150‑1 (MQ=255)
                                                                                                                           cAGAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCttt                          <  1:245292/150‑1 (MQ=255)
                                                                                                                             gAATGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTgg                        >  1:4845325/1‑150 (MQ=255)
                                                                                                                              aaTGGAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGt                       <  2:312371/150‑1 (MQ=255)
                                                                                                                                 ggAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGcc                    >  2:2486535/1‑150 (MQ=255)
                                                                                                                                 ggAAATATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGcc                    >  2:3310083/1‑150 (MQ=255)
                                                                                                                                    aaTATCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTg                 >  2:4654523/1‑150 (MQ=255)
                                                                                                                                     atatCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa                <  2:4889473/150‑1 (MQ=255)
                                                                                                                                     atatCAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGa                >  1:3830541/1‑150 (MQ=255)
                                                                                                                                         cAAACGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCt            <  2:4241640/150‑1 (MQ=255)
                                                                                                                                          aaaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCta           <  2:131971/150‑2 (MQ=255)
                                                                                                                                           aaCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTgg          >  1:3337161/1‑150 (MQ=255)
                                                                                                                                           aaCGTCCCTCCACCGAGGTCANAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCAGCAGNGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGANGCTgg          >  2:3354356/1‑150 (MQ=255)
                                                                                                                                            aCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTggg         >  1:4444534/1‑150 (MQ=255)
                                                                                                                                            aCGTCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTggg         <  1:4979325/150‑1 (MQ=255)
                                                                                                                                               tCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTACGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCAt      >  2:2569532/1‑150 (MQ=255)
                                                                                                                                               tCCCTCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCAt      <  2:582702/150‑1 (MQ=255)
                                                                                                                                                   gCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGGTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTc  >  1:1886074/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTACGATCAGGAATGCCCGCCTGCTCGCGAATGTAATCCGCCAGCAGATGGTCGAAATCATCGCCGCCGAGCGCGGAATCACCGCCGGTTGCCAGCACTTCAAACACGCCGCGACTTAAGCGCAGAATGGAAATATCAAACGTCCCGCCACCGAGGTCATAAACGGCGATCACGCCTTCCTGACCGGAATCCAGCCCGTAGGCGATAGCCGCAGCGGTCGGTTCGTTAAGTAAGCGCAGGACGTGAAGGCCCGCCAGACGCGCCGCGTCTTTGGTGCCCTGACGCTGGGCATCGTC  >  CP009273/2651518‑2651814

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: