Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,065,123 C→T 93.4% G287S (GGT→AGT)  pgk ← phosphoglycerate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,065,1230CT93.4% 622.3 / 28.0 197G287S (GGT→AGT) pgkphosphoglycerate kinase
Reads supporting (aligned to +/- strand):  ref base C (8/5);  new base T (88/96);  total (96/101)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 198 total aligned reads displayed.

AACGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAG  >  CP009273/3064978‑3065270
                                                                                                                                                 |                                                                                                                                                   
aaCGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGa                                                                                                                                                   <  2:318891/148‑1 (MQ=255)
  cGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGata                                                                                                                                                 <  2:1589946/148‑1 (MQ=255)
  cGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGata                                                                                                                                                 >  1:2154772/1‑148 (MQ=255)
  cGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGata                                                                                                                                                 >  1:2031806/1‑148 (MQ=255)
    cTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCTGAATTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATc                                                                                                                                               <  2:1643521/148‑1 (MQ=255)
    cTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATc                                                                                                                                               >  1:588724/1‑148 (MQ=255)
     ttCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATcc                                                                                                                                              <  2:1838255/148‑1 (MQ=255)
     ttCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATcc                                                                                                                                              <  1:1585985/148‑1 (MQ=255)
      tCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCa                                                                                                                                             >  2:2226343/1‑148 (MQ=255)
      tCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCa                                                                                                                                             <  1:196203/148‑1 (MQ=255)
        gCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAgg                                                                                                                                           <  2:1655325/148‑1 (MQ=255)
           gTCTGCGATAGCGTTAGCCACGATTTCTGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATc                                                                                                                                        <  1:1238139/148‑1 (MQ=255)
             cTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTg                                                                                                                                      <  2:2367908/148‑1 (MQ=255)
               gCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCt                                                                                                                                    >  2:461902/1‑148 (MQ=255)
                     gCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAg                                                                                                                              >  2:3151913/1‑148 (MQ=255)
                       gTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCt                                                                                                                            <  1:1752587/148‑1 (MQ=255)
                        ttAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCtt                                                                                                                           >  1:521523/1‑148 (MQ=255)
                        ttAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCtt                                                                                                                           >  1:1276530/1‑148 (MQ=255)
                        ttAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCtt                                                                                                                           >  1:1276431/1‑148 (MQ=255)
                         tAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCttt                                                                                                                          >  2:3072854/1‑148 (MQ=255)
                         tAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGCAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCttt                                                                                                                          >  2:2029207/1‑148 (MQ=255)
                          aGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTc                                                                                                                         <  2:1281424/148‑1 (MQ=255)
                           gCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTca                                                                                                                        <  2:536974/148‑1 (MQ=255)
                            ccACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTcac                                                                                                                       <  2:1532273/148‑1 (MQ=255)
                                gATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCg                                                                                                                   <  1:1296669/148‑1 (MQ=255)
                                 agttCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGt                                                                                                                  <  2:1978181/146‑1 (MQ=255)
                                 aTTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGt                                                                                                                  <  1:1569238/148‑1 (MQ=255)
                                    tCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTaa                                                                                                               <  2:359289/148‑1 (MQ=255)
                                        tACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGa                                                                                                           >  1:1399646/1‑148 (MQ=255)
                                         aCCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAt                                                                                                          <  2:2058057/148‑1 (MQ=255)
                                            tttGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTc                                                                                                       <  1:3180674/148‑1 (MQ=255)
                                            tttGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTc                                                                                                       >  1:2583052/1‑148 (MQ=255)
                                              tGCTGAATTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAg                                                                                                     <  2:601792/148‑1 (MQ=255)
                                              tGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAg                                                                                                     <  2:2538103/148‑1 (MQ=255)
                                              tGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAg                                                                                                     <  2:3136300/148‑1 (MQ=255)
                                              tGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAg                                                                                                     <  1:3166687/148‑1 (MQ=255)
                                               gCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAgg                                                                                                    <  1:2263695/148‑1 (MQ=255)
                                                 ggAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGt                                                                                                  <  2:453244/148‑1 (MQ=255)
                                                  gAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTa                                                                                                 >  2:432477/1‑148 (MQ=255)
                                                   aaGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAg                                                                                                >  1:1842650/1‑148 (MQ=255)
                                                     gTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGcc                                                                                              >  1:337729/1‑148 (MQ=255)
                                                      ttCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCg                                                                                             >  1:2921005/1‑148 (MQ=255)
                                                       tCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCgg                                                                                            >  2:1917188/1‑148 (MQ=255)
                                                         gggAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTg                                                                                          <  2:568600/148‑1 (MQ=255)
                                                             aTTCGAACACGCCCACCGGACCGGTCCACAGAATGGTTTTCGCAGTCTTCCGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGt                                                                                      <  1:2326627/148‑1 (MQ=255)
                                                             aTTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGt                                                                                      <  1:668285/148‑1 (MQ=255)
                                                              ttCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGtt                                                                                     >  2:700054/1‑148 (MQ=255)
                                                               tCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGttt                                                                                    <  2:2515780/148‑1 (MQ=255)
                                                               tCGAACACGCCAACCGGACCGGTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGttt                                                                                    <  2:2515806/148‑1 (MQ=255)
                                                                cGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTc                                                                                   >  2:1744217/1‑148 (MQ=255)
                                                                  aaCACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCag                                                                                 <  1:2079528/148‑1 (MQ=255)
                                                                   acacGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCaga                                                                                >  1:1692587/1‑148 (MQ=255)
                                                                   acacGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCaga                                                                                >  1:1692591/1‑148 (MQ=255)
                                                                    cacGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAATCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCAAATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCagag                                                                               >  2:1467413/1‑148 (MQ=255)
                                                                    cacGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCagag                                                                               <  2:2109715/148‑1 (MQ=255)
                                                                    cacGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCagag                                                                               >  2:2422183/1‑148 (MQ=255)
                                                                     acGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAgaga                                                                              <  2:3065962/148‑1 (MQ=255)
                                                                          aaCCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCg                                                                         <  1:3237633/148‑1 (MQ=255)
                                                                            ccGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGt                                                                       >  2:2876636/1‑148 (MQ=255)
                                                                            ccGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGt                                                                       >  1:3269612/1‑148 (MQ=255)
                                                                            ccGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGt                                                                       >  1:1332080/1‑148 (MQ=255)
                                                                                aCCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCt                                                                   >  1:2389428/1‑148 (MQ=255)
                                                                                 ccGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTa                                                                  <  2:1276431/148‑1 (MQ=255)
                                                                                 ccGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTa                                                                  <  2:1276530/148‑1 (MQ=255)
                                                                                   tttCCTCAGAATGGTTTTCGCATTCTTCAGGATTTCCGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTGGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTAcg                                                                <  2:3215465/147‑1 (MQ=255)
                                                                                   gTTCCACAGAATGGTTTTCTCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTAcg                                                                <  2:3011853/148‑1 (MQ=255)
                                                                                           gAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATc                                                        >  1:360836/1‑148 (MQ=255)
                                                                                             aTGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAg                                                      >  1:3064915/1‑148 (MQ=255)
                                                                                              tGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGa                                                     <  1:2731083/148‑1 (MQ=255)
                                                                                              tGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGa                                                     <  1:559726/148‑1 (MQ=255)
                                                                                                   ttCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGaa                                                <  1:2510464/148‑1 (MQ=255)
                                                                                                       cATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCgg                                            >  2:8135/1‑148 (MQ=255)
                                                                                                        aTTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCggg                                           >  2:251055/1‑148 (MQ=255)
                                                                                                        aTTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCggg                                           >  2:3278544/1‑148 (MQ=255)
                                                                                                         ttcttcAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGa                                          <  2:518210/148‑1 (MQ=255)
                                                                                                         ttcttcAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGa                                          <  2:2505704/148‑1 (MQ=255)
                                                                                                             tcAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGtt                                      >  1:863367/1‑148 (MQ=255)
                                                                                                             tcAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGtt                                      >  1:409817/1‑148 (MQ=255)
                                                                                                               aGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGc                                    >  2:1011969/1‑148 (MQ=255)
                                                                                                                ggATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCa                                   <  2:1679843/148‑1 (MQ=255)
                                                                                                                   tttCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAggt                                >  2:2610786/1‑148 (MQ=255)
                                                                                                                   tttCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAggt                                >  1:635945/1‑148 (MQ=255)
                                                                                                                   tttCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAggt                                <  1:837214/148‑1 (MQ=255)
                                                                                                                    ttCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAggtg                               >  2:493974/1‑148 (MQ=255)
                                                                                                                     tCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAggtgg                              >  2:1648659/1‑148 (MQ=255)
                                                                                                                       aGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTc                            <  2:3104224/148‑1 (MQ=255)
                                                                                                                        gCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTca                           >  2:1106316/1‑148 (MQ=255)
                                                                                                                        gCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTca                           >  2:2758249/1‑148 (MQ=255)
                                                                                                                         ccAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTcag                          <  1:1541098/148‑1 (MQ=255)
                                                                                                                             ttCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGa                      >  1:1627895/1‑148 (MQ=255)
                                                                                                                              tCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGAc                     >  2:3150404/1‑148 (MQ=255)
                                                                                                                               ccTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACg                    >  2:632672/1‑148 (MQ=255)
                                                                                                                                 tGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGtt                  <  2:1228840/148‑1 (MQ=255)
                                                                                                                                 tGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGtt                  >  2:721659/1‑148 (MQ=255)
                                                                                                                                  gAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGttt                 <  2:2921005/148‑1 (MQ=255)
                                                                                                                                       gAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCtt            <  2:1815474/148‑1 (MQ=255)
                                                                                                                                           cATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTc        >  1:1037463/1‑148 (MQ=255)
                                                                                                                                             tCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCaa      >  2:1361/1‑148 (MQ=255)
                                                                                                                                             tCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCaa      >  1:1585035/1‑148 (MQ=255)
                                                                                                                                                 cGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAg  <  1:1094833/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AACGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAACCAG  >  CP009273/3064978‑3065270

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: