Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,933,166 C→A 6.1% R421S (CGT→AGT)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,933,1660CA 6.1% 991.2 / 21.7 262R421S (CGT→AGT) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base C (128/118);  new base A (7/9);  total (135/127)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 268 total aligned reads displayed.

AGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTC  >  CP009273/1933021‑1933314
                                                                                                                                                 |                                                                                                                                                    
 gCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTgg                                                                                                                                                  >  1:2490644/1‑149 (MQ=255)
 gCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTgg                                                                                                                                                  <  1:3990519/149‑1 (MQ=255)
 gCTAACCACCTAAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTgg                                                                                                                                                  <  2:2967946/149‑1 (MQ=255)
  cTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGc                                                                                                                                                 >  1:265781/1‑149 (MQ=255)
    aaCCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGt                                                                                                                                               >  1:488794/1‑149 (MQ=255)
     accaccTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGtt                                                                                                                                              >  1:3682168/1‑149 (MQ=255)
      ccaccTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTa                                                                                                                                             >  1:3044627/1‑149 (MQ=255)
       caccTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAc                                                                                                                                            <  1:766735/149‑1 (MQ=255)
           tGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCg                                                                                                                                        >  1:1805343/1‑149 (MQ=255)
             aaaGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGt                                                                                                                                      >  1:2322788/1‑149 (MQ=255)
                   gTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACtt                                                                                                                                >  2:3770550/1‑149 (MQ=255)
                   gTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACtt                                                                                                                                <  2:3140879/149‑1 (MQ=255)
                   gTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACtg                                                                                                                                >  2:911212/1‑148 (MQ=255)
                     tACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTg                                                                                                                              <  2:709452/149‑1 (MQ=255)
                        ggCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATa                                                                                                                           <  2:2585233/149‑1 (MQ=255)
                        ggCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATa                                                                                                                           <  1:2980971/149‑1 (MQ=255)
                         gCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAg                                                                                                                          <  1:2927933/149‑1 (MQ=255)
                           gATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAgcg                                                                                                                        >  2:1581842/1‑149 (MQ=255)
                           gATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAgcg                                                                                                                        >  1:3248518/1‑149 (MQ=255)
                            atcatcACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAgcac                                                                                                                       >  1:1038483/1‑147 (MQ=255)
                             tcatcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTATCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCt                                                                                                                      <  1:649216/149‑1 (MQ=255)
                              catcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTa                                                                                                                     >  1:1412061/1‑149 (MQ=255)
                               atcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTaa                                                                                                                    <  1:3285983/149‑1 (MQ=255)
                                tcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAAt                                                                                                                   <  1:462500/149‑1 (MQ=255)
                                tcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAAt                                                                                                                   >  1:2304105/1‑149 (MQ=255)
                                 caCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATg                                                                                                                  >  2:1139237/1‑149 (MQ=255)
                                  aCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGa                                                                                                                 <  2:2926678/149‑1 (MQ=255)
                                    cATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACg                                                                                                               >  1:3451434/1‑149 (MQ=255)
                                     aTGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAAGGACgg                                                                                                              <  2:2760781/149‑1 (MQ=255)
                                      tGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGc                                                                                                             <  1:709150/149‑1 (MQ=255)
                                       gACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCg                                                                                                            >  1:2747993/1‑149 (MQ=255)
                                       gACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCg                                                                                                            >  1:789508/1‑149 (MQ=255)
                                        aCCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGt                                                                                                           >  1:80805/1‑149 (MQ=255)
                                         ccGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTa                                                                                                          >  1:2642686/1‑149 (MQ=255)
                                            aaTCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTagca                                                                                                       >  1:3821267/1‑149 (MQ=255)
                                            aaTCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTagca                                                                                                       <  2:3552488/149‑1 (MQ=255)
                                            aaTCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTagca                                                                                                       >  1:2609653/1‑149 (MQ=255)
                                                  gTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGcc                                                                                                 <  1:2098526/149‑1 (MQ=255)
                                                        ccGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGaa                                                                                           <  2:1103884/149‑1 (MQ=255)
                                                        ccGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGaa                                                                                           <  2:524381/149‑1 (MQ=255)
                                                          gcgcTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGc                                                                                         <  2:1897474/149‑1 (MQ=255)
                                                              tgatgaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGtt                                                                                     <  1:3129955/149‑1 (MQ=255)
                                                                atgaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTaa                                                                                   <  1:911757/149‑1 (MQ=255)
                                                                 tgaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAAt                                                                                  <  2:3599721/149‑1 (MQ=255)
                                                                  gaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATc                                                                                 <  1:331879/149‑1 (MQ=255)
                                                                   aCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGACGCGGTTAATct                                                                                >  2:1852382/1‑149 (MQ=255)
                                                                     cTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATctgc                                                                              >  2:3848356/1‑149 (MQ=255)
                                                                     cTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATctgc                                                                              >  2:757192/1‑149 (MQ=255)
                                                                     cTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATctgc                                                                              >  1:1608123/1‑149 (MQ=255)
                                                                         cGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTgcg                                                                          >  1:1011463/1‑149 (MQ=255)
                                                                           tATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGcgcg                                                                        >  2:2181067/1‑149 (MQ=255)
                                                                            aTCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGa                                                                       <  2:1433750/149‑1 (MQ=255)
                                                                            aTCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGa                                                                       <  1:898072/149‑1 (MQ=255)
                                                                             tCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGAt                                                                      <  1:2214579/149‑1 (MQ=255)
                                                                             tCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGAt                                                                      <  2:2460665/149‑1 (MQ=255)
                                                                              cAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATa                                                                     >  2:2607606/1‑149 (MQ=255)
                                                                                gCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATTACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATaaa                                                                   <  1:1855936/149‑1 (MQ=255)
                                                                                gCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATaaa                                                                   >  1:323929/1‑149 (MQ=255)
                                                                                 ctctGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAg                                                                  <  2:1301118/149‑1 (MQ=255)
                                                                                  tctGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAgg                                                                 >  2:1162895/1‑149 (MQ=255)
                                                                                   ctGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGt                                                                <  2:1412061/149‑1 (MQ=255)
                                                                                      gTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTAc                                                             <  1:3776809/149‑1 (MQ=255)
                                                                                       tCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGTCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACt                                                            <  2:459068/149‑1 (MQ=255)
                                                                                         tGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTg                                                          <  1:98994/149‑1 (MQ=255)
                                                                                          gCCAATTTTCGCCATGTCGCGCCATGAAAGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGa                                                         >  2:3584914/1‑149 (MQ=255)
                                                                                          gCCAATGTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGTCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGa                                                         <  2:3996878/149‑1 (MQ=255)
                                                                                           ccAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGAt                                                        <  1:2696904/149‑1 (MQ=255)
                                                                                            cAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATg                                                       >  2:2249223/1‑149 (MQ=255)
                                                                                             aaTTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGt                                                      >  2:2007249/1‑149 (MQ=255)
                                                                                              aTTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTc                                                     <  1:2573135/149‑1 (MQ=255)
                                                                                              aTTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTc                                                     >  1:2869576/1‑149 (MQ=255)
                                                                                               ttttCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCt                                                    <  1:394641/149‑1 (MQ=255)
                                                                                                     ccATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGAc                                              <  2:1927623/149‑1 (MQ=255)
                                                                                                     ccATGTCGCGCCATGAACGGACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGAc                                              <  2:2283362/149‑1 (MQ=255)
                                                                                                      cATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGAcc                                             >  2:3181439/1‑149 (MQ=255)
                                                                                                       aTGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCt                                            <  1:3051541/149‑1 (MQ=255)
                                                                                                        tGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTg                                           <  1:3454804/149‑1 (MQ=255)
                                                                                                         gTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTgg                                          <  1:2173711/149‑1 (MQ=255)
                                                                                                          tCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGt                                         <  1:3758323/149‑1 (MQ=255)
                                                                                                           cgcgCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTg                                        >  1:997252/1‑149 (MQ=255)
                                                                                                            gcgcCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGa                                       >  2:2537047/1‑149 (MQ=255)
                                                                                                               ccATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTg                                    >  1:1984223/1‑149 (MQ=255)
                                                                                                                   gAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCAc                                >  1:1593343/1‑149 (MQ=255)
                                                                                                                      cGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCa                             >  1:3835754/1‑149 (MQ=255)
                                                                                                                       gTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAg                            >  1:81245/1‑149 (MQ=255)
                                                                                                                        tACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAgg                           >  2:1023321/1‑149 (MQ=255)
                                                                                                                          cGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGc                         <  2:2747993/149‑1 (MQ=255)
                                                                                                                          cGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGc                         <  2:789498/149‑1 (MQ=255)
                                                                                                                             tGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGAc                      <  1:3584914/149‑1 (MQ=255)
                                                                                                                             tGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGAc                      <  1:1536763/149‑1 (MQ=255)
                                                                                                                                  cTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtgat                 <  2:2323292/149‑1 (MQ=255)
                                                                                                                                   tGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtgatg                >  1:317087/1‑149 (MQ=255)
                                                                                                                                   tGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtgatg                <  2:3412648/149‑1 (MQ=255)
                                                                                                                                   tGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtgatg                >  2:827861/1‑149 (MQ=255)
                                                                                                                                        gCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTAc           <  2:303585/149‑1 (MQ=255)
                                                                                                                                        gCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTAc           <  1:3874217/149‑1 (MQ=255)
                                                                                                                                               aTCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGt    <  2:2575888/149‑1 (MQ=255)
                                                                                                                                                tCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGtt   <  2:4021152/149‑1 (MQ=255)
                                                                                                                                                 cGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTc  <  1:2775075/149‑1 (MQ=255)
                                                                                                                                                 cGTGGCGTTACGCCGGTGCACTTTGATAGCGCGAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTc  <  1:2775105/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGTTC  >  CP009273/1933021‑1933314

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