Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,930,891 G→T 6.2% D80Y (GAT→TAT)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,930,8910GT 6.2% 1042.6 / 31.9 281D80Y (GAT→TAT) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base G (117/146);  new base T (13/5);  total (130/151)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

CCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGT  >  CP009273/3930742‑3931040
                                                                                                                                                     |                                                                                                                                                     
 cGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGa                                                                                                                                                      <  2:5136817/150‑1 (MQ=255)
   tAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTtata                                                                                                                                                    >  2:4266489/1‑150 (MQ=255)
       cGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                                                                                >  2:4501256/1‑150 (MQ=255)
       cGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAt                                                                                                                                                >  2:3763742/1‑150 (MQ=255)
         gTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACGTATAACATTg                                                                                                                                              >  2:4940084/1‑150 (MQ=255)
          tAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGa                                                                                                                                             >  1:4251222/1‑150 (MQ=255)
          tAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGa                                                                                                                                             >  1:2356043/1‑150 (MQ=255)
           aaCCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGat                                                                                                                                            <  2:1451849/150‑1 (MQ=255)
           aaCCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGat                                                                                                                                            <  1:3433798/150‑1 (MQ=255)
            aCCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGACTATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGata                                                                                                                                           >  1:613470/1‑150 (MQ=255)
            aCCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGata                                                                                                                                           >  1:669773/1‑150 (MQ=255)
             ccACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGatat                                                                                                                                          <  2:928868/150‑1 (MQ=255)
               aCTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTa                                                                                                                                        <  2:3383769/150‑1 (MQ=255)
                 tATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACt                                                                                                                                      >  1:1406094/1‑150 (MQ=255)
                 tATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACt                                                                                                                                      <  1:791630/150‑1 (MQ=255)
                  aTCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTc                                                                                                                                     >  2:3942044/1‑150 (MQ=255)
                   tCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTcc                                                                                                                                    <  1:4061700/150‑1 (MQ=255)
                    cAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCg                                                                                                                                   <  1:971458/150‑1 (MQ=255)
                      ggTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGt                                                                                                                                 <  1:319178/150‑1 (MQ=255)
                       gTTGCATTTGGCGGCAAAGGCGCTAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTc                                                                                                                                <  2:3128690/150‑1 (MQ=255)
                        ttGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATCTTACTCCGGTCa                                                                                                                               >  1:3143041/1‑150 (MQ=255)
                        ttGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCa                                                                                                                               <  1:894430/150‑1 (MQ=255)
                        ttGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCa                                                                                                                               <  1:1174243/150‑1 (MQ=255)
                        ttGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCACTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCa                                                                                                                               <  1:1174299/150‑1 (MQ=255)
                          gCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGc                                                                                                                             <  2:1178876/150‑1 (MQ=255)
                             tttGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTg                                                                                                                          <  1:950147/150‑1 (MQ=255)
                                ggcggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATc                                                                                                                       <  2:3615992/150‑1 (MQ=255)
                                   ggcAAAGGCGCGAATCAGGGTGTGGCTGCTTGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTTATAACATTGATATTACTCCGGTCAGCGTGATCaaa                                                                                                                    <  2:1540830/150‑1 (MQ=255)
                                   ggcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCaaa                                                                                                                    >  1:1500327/1‑150 (MQ=255)
                                    gcAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAg                                                                                                                   <  1:1793661/150‑1 (MQ=255)
                                      aaaGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACAGTGATATTACTCCGGTCAGCGTGATCAATGGc                                                                                                                 >  1:4839065/1‑150 (MQ=255)
                                        aGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGa                                                                                                               >  2:3295348/1‑150 (MQ=255)
                                        aGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGa                                                                                                               >  1:493238/1‑150 (MQ=255)
                                        aGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATAGATATTACTCCGGTCAGCGTGATCAAAGGCGa                                                                                                               >  1:4287662/1‑150 (MQ=255)
                                         gnCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGaa                                                                                                              >  1:3566197/1‑150 (MQ=255)
                                          gcgcGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAAt                                                                                                             <  2:1718932/150‑1 (MQ=255)
                                           cgcgAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATc                                                                                                            <  1:3261931/150‑1 (MQ=255)
                                                aTCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAg                                                                                                       <  1:743558/150‑1 (MQ=255)
                                                 tCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAgg                                                                                                      >  2:1975393/1‑150 (MQ=255)
                                                  cAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgt                                                                                                     >  1:2933306/1‑150 (MQ=255)
                                                   aGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTTATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgtg                                                                                                    <  1:3717964/150‑1 (MQ=255)
                                                    ggCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCTTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgtgt                                                                                                   <  1:1734539/150‑1 (MQ=255)
                                                    ggCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTTATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGgtgt                                                                                                   <  2:3303172/150‑1 (MQ=255)
                                                         tgGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCTATGGCGAATCAACAGGTGTGgcgc                                                                                              >  2:2675169/1‑150 (MQ=255)
                                                          ggCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCt                                                                                             <  2:2546211/150‑1 (MQ=255)
                                                          ggCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGAAGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCt                                                                                             >  1:1254458/1‑150 (MQ=255)
                                                           gctgctGGCCGTAGCGGTTCGAATATAGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCAACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTg                                                                                            <  2:3484132/150‑1 (MQ=255)
                                                             tgctgGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTTATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGAt                                                                                          >  1:1704985/1‑150 (MQ=255)
                                                             tgctgGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGAt                                                                                          >  1:2647041/1‑150 (MQ=255)
                                                                cgggCGTAGCGGTGCGACGATCGCGTGTATAGCCAGGACGGGTGATGCCAGCAGGGGTGAGAGCGGTCGCCAGCAGCTCGCCCCTGATAACATCGATTGTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGAtttt                                                                                       <  2:3553282/149‑1 (MQ=255)
                                                                 gggCGTAGCGGTGCGAAAATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGAttttt                                                                                      <  1:4966413/150‑1 (MQ=255)
                                                                  ggCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTg                                                                                     <  1:3929830/150‑1 (MQ=255)
                                                                    cGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGtt                                                                                   <  1:3016589/150‑1 (MQ=255)
                                                                     gTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTa                                                                                  >  2:814583/1‑150 (MQ=255)
                                                                       aGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAAt                                                                                >  2:1790379/1‑150 (MQ=255)
                                                                       aGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAAt                                                                                >  2:2155966/1‑150 (MQ=255)
                                                                        gCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATg                                                                               >  2:4100524/1‑150 (MQ=255)
                                                                        gCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATg                                                                               <  2:4688030/150‑1 (MQ=255)
                                                                        gCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATg                                                                               <  2:2081537/150‑1 (MQ=255)
                                                                               gAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAgg                                                                        <  1:3388957/150‑1 (MQ=255)
                                                                                 atatCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTg                                                                      <  1:2078101/150‑1 (MQ=255)
                                                                                 atatCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTg                                                                      <  2:2254192/150‑1 (MQ=255)
                                                                                   atCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTgag                                                                    >  2:1959317/1‑150 (MQ=255)
                                                                                       cgTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATg                                                                <  2:234046/150‑1 (MQ=255)
                                                                                        gTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACCTTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGt                                                               <  2:4185293/150‑1 (MQ=255)
                                                                                        gTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGt                                                               >  2:2457658/1‑150 (MQ=255)
                                                                                          ttATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTTATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCa                                                             >  2:2929612/1‑150 (MQ=255)
                                                                                            aTTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATc                                                           >  1:2456161/1‑150 (MQ=255)
                                                                                              tGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGAGAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCgg                                                         <  1:5071736/150‑1 (MQ=255)
                                                                                               gCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAAAGTCATCGGt                                                        >  1:2813319/1‑150 (MQ=255)
                                                                                                ccTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGACTCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCAGCGGTa                                                       >  2:4590406/1‑150 (MQ=255)
                                                                                                ccTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTa                                                       >  1:4560011/1‑150 (MQ=255)
                                                                                                   gTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTc                                                    <  1:5208675/150‑1 (MQ=255)
                                                                                                    tACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCa                                                   >  2:5071358/1‑150 (MQ=255)
                                                                                                      cGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATg                                                 <  2:3040653/150‑1 (MQ=255)
                                                                                                        ggTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGcc                                               <  2:626655/150‑1 (MQ=255)
                                                                                                        ggTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATAAAAGGCGAATCAAAATGTTTGGCGCTGATTTTTGTTAATGGCGTAGGTGAGAATGTCATCTGTTTTCATGcc                                               <  2:869828/150‑1 (MQ=255)
                                                                                                           gatgaCAGCATTGGTGAGATCGTTCGCCAGCAGCACGACACAGATAACATTGATATTACTCCGGTCATCGTGATCAAATCCGAATCAACAGTATAGGCGCTGATATTTGTTAATCGCGAATGAGACAAAGTCATCGGTATTAATGCCTgc                                            >  2:2554337/1‑150 (MQ=255)
                                                                                                            atgaCAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGgcg                                           >  2:4766360/1‑150 (MQ=255)
                                                                                                              gaCAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCt                                         <  1:1558216/150‑1 (MQ=255)
                                                                                                               aCAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTa                                        <  2:4387627/150‑1 (MQ=255)
                                                                                                                cAGCATTGGAGAGAGCGTTCGCCATCAGATCGCCACTGATAACATAGATAATACTACGGTCAGAGAGATAAAAGGCGAATCAACAGGTGTGGCGCAGATTTTTGTTAAATGCTCAGATGATAATGTCAGCGGTATACATGCAGGCGCTaa                                       >  2:4873661/1‑150 (MQ=255)
                                                                                                                   cATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAAtgc                                    >  1:3582641/1‑150 (MQ=255)
                                                                                                                       ggTGAGAGCGTTCGCCAGCAGCTCGCCACTGATACCATTGATATTACTCCGGTCAGCGTGATCAAAGGCTAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGcc                                <  2:3989489/150‑1 (MQ=255)
                                                                                                                          gagaGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCtt                             <  1:1241865/150‑1 (MQ=255)
                                                                                                                           agagCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCttt                            <  2:3985670/150‑1 (MQ=255)
                                                                                                                           agagCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCttt                            <  1:4786949/150‑1 (MQ=255)
                                                                                                                               cGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTcccc                        <  1:2133656/150‑1 (MQ=255)
                                                                                                                                gTTCGCCAGCAGCTCGCCACTTATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCg                       >  1:4409096/1‑150 (MQ=255)
                                                                                                                                   cGCAAGCAGCTCGCAACTTATAACATTTATATTACACCGGTCAGCGTGATCAAAGGCGAATCAAAAGGTGTGCCGCTGATTTATGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCggag                    >  2:56587/1‑148 (MQ=255)
                                                                                                                                      cagcagCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCtg                 >  1:4667521/1‑150 (MQ=255)
                                                                                                                                      cagcagCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCtg                 >  1:2211465/1‑150 (MQ=255)
                                                                                                                                          agcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGGGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCtggtgg             >  1:2768925/1‑150 (MQ=255)
                                                                                                                                              cGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAgc         >  1:4882927/1‑150 (MQ=255)
                                                                                                                                              cGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAgc         >  1:4168541/1‑150 (MQ=255)
                                                                                                                                              cGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAgc         <  1:3989016/150‑1 (MQ=255)
                                                                                                                                                  aCTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTACTGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCaa     <  2:687839/150‑1 (MQ=255)
                                                                                                                                                    tGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCACCg   >  2:4405875/1‑150 (MQ=255)
                                                                                                                                                    tGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACg   >  1:1851317/1‑150 (MQ=255)
                                                                                                                                                    tGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACg   >  2:1962773/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTAACCGGTAACCACTATCAGGTTGCATTTGGCGGCAAAGGCGCGAATCAGGCTGTGGCTGCTGGGCGTAGCGGTGCGAATATCGCGTTTATTGCCTGTACGGGTGATGACAGCATTGGTGAGAGCGTTCGCCAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGT  >  CP009273/3930742‑3931040

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: ‑