Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,174,009 G→T 5.5% R77L (CGC→CTC)  ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,174,0090GT 5.5% 1316.7 / 35.0 354R77L (CGC→CTC) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (189/143);  new base T (14/8);  total (203/151)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

CGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGG  >  CP009273/3173860‑3174158
                                                                                                                                                     |                                                                                                                                                     
 ggCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGc                                                                                                                                                      >  2:1727592/1‑150 (MQ=255)
   cGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAt                                                                                                                                                    <  1:2007552/150‑1 (MQ=255)
    gAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAtt                                                                                                                                                   >  2:105059/1‑150 (MQ=255)
       ccGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCa                                                                                                                                                >  2:5201097/1‑150 (MQ=255)
       ccGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCa                                                                                                                                                <  2:4330406/150‑1 (MQ=255)
         gTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCaaa                                                                                                                                              <  2:1958948/150‑1 (MQ=255)
           aCTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAGGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAac                                                                                                                                            <  2:4237662/150‑1 (MQ=255)
                   gAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTgg                                                                                                                                    >  2:2394574/1‑150 (MQ=255)
                   gAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTgg                                                                                                                                    <  1:3617206/150‑1 (MQ=255)
                      gATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGt                                                                                                                                 >  1:526460/1‑150 (MQ=255)
                           ttACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCTCATTCCAAAACACACCTGGAGTTTtgt                                                                                                                            <  2:3499456/150‑1 (MQ=255)
                           ttACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTtgt                                                                                                                            >  1:3238256/1‑150 (MQ=255)
                            tactacAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTtgtg                                                                                                                           <  1:2828759/150‑1 (MQ=255)
                             actacAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTgc                                                                                                                          >  2:4957927/1‑150 (MQ=255)
                              ctacAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTgcg                                                                                                                         >  2:4254138/1‑150 (MQ=255)
                              ctacAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTgcg                                                                                                                         >  2:1874421/1‑150 (MQ=255)
                                acAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGcc                                                                                                                       >  1:3371294/1‑150 (MQ=255)
                                acAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGcc                                                                                                                       >  1:748494/1‑150 (MQ=255)
                                 catgTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACTCCTGGAGTTTTGTGCGCCa                                                                                                                      <  2:4110648/147‑1 (MQ=255)
                                 cAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCa                                                                                                                      >  1:1164964/1‑150 (MQ=255)
                                  aaGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAg                                                                                                                     >  2:912187/1‑150 (MQ=255)
                                     ttGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCa                                                                                                                  <  2:4019137/150‑1 (MQ=255)
                                      tGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCAt                                                                                                                 >  1:3061133/1‑150 (MQ=255)
                                       gACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATg                                                                                                                <  2:206322/150‑1 (MQ=255)
                                          ccTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACGAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCt                                                                                                             >  2:201819/1‑150 (MQ=255)
                                           cTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTg                                                                                                            >  2:177626/1‑150 (MQ=255)
                                            tCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGa                                                                                                           <  1:1010110/150‑1 (MQ=255)
                                             cGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGAc                                                                                                          <  1:33048/150‑1 (MQ=255)
                                              gCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACg                                                                                                         >  2:2996491/1‑150 (MQ=255)
                                                  aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAACCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCAcc                                                                                                     >  2:533419/1‑150 (MQ=255)
                                                  aGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCAcc                                                                                                     >  2:3409729/1‑150 (MQ=255)
                                                   ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCa                                                                                                    <  1:4074046/150‑1 (MQ=255)
                                                   ggTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGAACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACAAACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCAACa                                                                                                    >  2:2263013/1‑150 (MQ=255)
                                                     ttGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGc                                                                                                  >  1:3310641/1‑150 (MQ=255)
                                                       gAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCTCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCa                                                                                                >  1:4471452/1‑150 (MQ=255)
                                                        aacaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAt                                                                                               <  2:5083518/147‑1 (MQ=255)
                                                        aaaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCAt                                                                                               <  1:4253987/150‑1 (MQ=255)
                                                           aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCCAATTCCGCATTCCAAAACACACCTGGAGTGTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGc                                                                                            >  2:4761806/1‑150 (MQ=255)
                                                            gCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCt                                                                                           >  2:2673887/1‑150 (MQ=255)
                                                            gCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCt                                                                                           >  1:2993558/1‑150 (MQ=255)
                                                                 aagaagTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTAtt                                                                                      >  1:2551537/1‑150 (MQ=255)
                                                                     agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                  >  2:1017029/1‑150 (MQ=255)
                                                                     agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                  >  1:3978128/1‑150 (MQ=255)
                                                                     agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                  >  1:5069068/1‑150 (MQ=255)
                                                                      gTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGt                                                                                 >  2:3896256/1‑150 (MQ=255)
                                                                         aCCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCtt                                                                              >  1:3192648/1‑150 (MQ=255)
                                                                         aCCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCtt                                                                              <  1:2167370/150‑1 (MQ=255)
                                                                          ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                             <  1:1269377/150‑1 (MQ=255)
                                                                            gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGAt                                                                           <  2:2636862/150‑1 (MQ=255)
                                                                             ccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATc                                                                          <  2:5172981/150‑1 (MQ=255)
                                                                               gAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTg                                                                        >  2:2191079/1‑150 (MQ=255)
                                                                                  cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                     <  1:1806368/150‑1 (MQ=255)
                                                                                  cTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCg                                                                     >  2:3888854/1‑150 (MQ=255)
                                                                                    gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                   <  2:761526/150‑1 (MQ=255)
                                                                                    gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                   >  1:922956/1‑150 (MQ=255)
                                                                                      aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCTCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                 >  1:4795478/1‑150 (MQ=255)
                                                                                      aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                 >  1:287193/1‑150 (MQ=255)
                                                                                      aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                 >  1:2925526/1‑150 (MQ=255)
                                                                                       ccAGATGTGCCTGAACGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGa                                                                <  2:3387525/150‑1 (MQ=255)
                                                                                       ccAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGa                                                                >  1:4687386/1‑150 (MQ=255)
                                                                                       ccAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGa                                                                <  1:4830461/150‑1 (MQ=255)
                                                                                        cAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGAt                                                               <  2:316675/150‑1 (MQ=255)
                                                                                        cAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGAt                                                               <  2:1228850/150‑1 (MQ=255)
                                                                                          gATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                             >  1:3466022/1‑150 (MQ=255)
                                                                                          gATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                             >  2:886478/1‑150 (MQ=255)
                                                                                             gtgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCAc                                                          >  1:1470175/1‑150 (MQ=255)
                                                                                              tgCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACt                                                         <  1:1391958/150‑1 (MQ=255)
                                                                                                ccTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTgg                                                       >  1:2063675/1‑150 (MQ=255)
                                                                                                 cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                      <  2:2793685/150‑1 (MQ=255)
                                                                                                 cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                      <  2:3680265/150‑1 (MQ=255)
                                                                                                 cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                      >  2:1351590/1‑150 (MQ=255)
                                                                                                  tGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTg                                                     >  2:3795185/1‑150 (MQ=255)
                                                                                                  agAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTg                                                     <  1:5183874/149‑1 (MQ=255)
                                                                                                       gtggtgGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCt                                                >  2:4522933/1‑150 (MQ=255)
                                                                                                       gtggtgGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCt                                                <  1:1327562/150‑1 (MQ=255)
                                                                                                        tggtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCTCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTa                                               >  1:2686674/1‑150 (MQ=255)
                                                                                                        tggtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTa                                               >  1:1628823/1‑150 (MQ=255)
                                                                                                         ggtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTaa                                              >  2:1003943/1‑150 (MQ=255)
                                                                                                           tggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAAc                                            <  2:4096124/150‑1 (MQ=255)
                                                                                                           aggAAAAAGTGATCGCCAGCGATGAGCTGATGACCCAAGTCCGCAGTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATTGCTGACTCACCAGCCATCGCTTTATTCGAGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAAc                                            <  2:2989379/149‑1 (MQ=255)
                                                                                                                    tGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaa                                   >  2:3744209/1‑150 (MQ=255)
                                                                                                                      aTCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAAt                                 <  1:1057121/150‑1 (MQ=255)
                                                                                                                       tCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATa                                <  1:3332255/150‑1 (MQ=255)
                                                                                                                            aGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                           <  2:48587/150‑1 (MQ=255)
                                                                                                                                aTGAGCTGATGACCAAATTCCGCATTCCAACACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTCAACTTCTGGAAAATAACGCCGATa                       <  2:1125490/150‑1 (MQ=255)
                                                                                                                                 tGAGCTGATGACCAAATTCCTCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATAc                      >  1:496900/1‑150 (MQ=255)
                                                                                                                                 tGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATAc                      >  2:4448522/1‑150 (MQ=255)
                                                                                                                                 tGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATAc                      <  1:382923/150‑1 (MQ=255)
                                                                                                                                 tGAGCTGATGACCAAATTCCGAATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGATTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATAc                      <  1:2368845/150‑1 (MQ=255)
                                                                                                                                       gatgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAAc                <  2:956636/150‑1 (MQ=255)
                                                                                                                                       gatgaCCAAAGTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAAc                <  1:3072544/150‑1 (MQ=255)
                                                                                                                                        atgaCCAAATTCCGCATTCCAATACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGACTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACg               <  2:4872126/150‑1 (MQ=255)
                                                                                                                                        atgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACg               <  2:2275079/150‑1 (MQ=255)
                                                                                                                                        atgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACg               >  1:2273239/1‑150 (MQ=255)
                                                                                                                                            ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc           <  2:3307003/150‑1 (MQ=255)
                                                                                                                                              naaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta         >  2:2501557/2‑150 (MQ=255)
                                                                                                                                              aaaTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACta         >  1:3515332/1‑150 (MQ=255)
                                                                                                                                                aTTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATCCGCCAACGTCACtata       >  2:4795111/1‑150 (MQ=255)
                                                                                                                                                 ttCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATAc      <  2:653202/150‑1 (MQ=255)
                                                                                                                                                    cGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAg   >  2:3611998/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGG  >  CP009273/3173860‑3174158

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: