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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,065,1230CT96.0% 639.3 / 13.4 198G287S (GGT→AGT) pgkphosphoglycerate kinase
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (103/87);  total (108/90)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 198 total aligned reads displayed.

AACGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAA  >  CP009273/3064978‑3065266
                                                                                                                                                 |                                                                                                                                               
 aCGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGat                                                                                                                                              <  1:950543/148‑1 (MQ=255)
 aCGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGat                                                                                                                                              <  1:1561473/148‑1 (MQ=255)
  cGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGata                                                                                                                                             <  2:2171700/148‑1 (MQ=255)
   gCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGatat                                                                                                                                            <  2:224106/148‑1 (MQ=255)
    cTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATc                                                                                                                                           <  1:2362022/148‑1 (MQ=255)
    cTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATc                                                                                                                                           >  1:1416478/1‑148 (MQ=255)
         cTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGa                                                                                                                                      >  2:3281987/1‑148 (MQ=255)
             cTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTg                                                                                                                                  >  2:691725/1‑148 (MQ=255)
              tGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGc                                                                                                                                 <  2:2974594/148‑1 (MQ=255)
               gCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCt                                                                                                                                <  1:1370757/148‑1 (MQ=255)
                cGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTc                                                                                                                               >  2:3181139/1‑148 (MQ=255)
                cGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTc                                                                                                                               <  2:725414/148‑1 (MQ=255)
                  aTAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGt                                                                                                                             >  2:2267574/1‑148 (MQ=255)
                  aTAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGt                                                                                                                             >  1:1143627/1‑148 (MQ=255)
                  aTAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGt                                                                                                                             <  2:2345130/148‑1 (MQ=255)
                  aTAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGt                                                                                                                             <  2:3206211/148‑1 (MQ=255)
                   tAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTc                                                                                                                            >  2:1868410/1‑148 (MQ=255)
                       gTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCt                                                                                                                        <  1:3258073/148‑1 (MQ=255)
                        ttAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCtt                                                                                                                       <  2:1183926/148‑1 (MQ=255)
                         tAGCCACGATTTCTGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCttt                                                                                                                      >  1:1133379/1‑148 (MQ=255)
                         tAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCttt                                                                                                                      <  1:361710/148‑1 (MQ=255)
                         tAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCttt                                                                                                                      >  2:2946570/1‑148 (MQ=255)
                         tAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCttt                                                                                                                      >  2:1770122/1‑148 (MQ=255)
                              aCGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACAt                                                                                                                 >  1:462802/1‑148 (MQ=255)
                                 aTTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGGTCTGCTCGTCAGCTTTCACATCGt                                                                                                              >  2:131567/1‑148 (MQ=255)
                                   ttCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTa                                                                                                            <  1:2103496/148‑1 (MQ=255)
                                       gTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAg                                                                                                        <  2:1098416/148‑1 (MQ=255)
                                       gTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAg                                                                                                        <  2:2873260/148‑1 (MQ=255)
                                         aCCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAt                                                                                                      >  1:945225/1‑148 (MQ=255)
                                           cTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAttt                                                                                                    >  2:1672879/1‑148 (MQ=255)
                                           cTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAttt                                                                                                    >  2:2832305/1‑148 (MQ=255)
                                           cTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAttt                                                                                                    >  1:1522022/1‑148 (MQ=255)
                                           cTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAttt                                                                                                    <  2:888718/148‑1 (MQ=255)
                                           cTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAttt                                                                                                    >  2:3537624/1‑148 (MQ=255)
                                            tttGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTc                                                                                                   <  2:244749/148‑1 (MQ=255)
                                            tttGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTc                                                                                                   <  1:942149/148‑1 (MQ=255)
                                             ttGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGAtttga                                                                                                  <  2:453336/148‑3 (MQ=255)
                                             ttGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCa                                                                                                  <  2:3135976/148‑1 (MQ=255)
                                              tGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAg                                                                                                 <  1:1717277/148‑1 (MQ=255)
                                               gCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAgg                                                                                                <  2:649691/148‑1 (MQ=255)
                                                cGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAggg                                                                                               >  1:3149205/1‑148 (MQ=255)
                                                cGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAggg                                                                                               >  1:2235098/1‑148 (MQ=255)
                                                cGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAggg                                                                                               <  2:3326271/148‑1 (MQ=255)
                                                cGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAggg                                                                                               <  2:202602/148‑1 (MQ=255)
                                                   aaGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAg                                                                                            >  1:1906273/1‑148 (MQ=255)
                                                   aaGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAg                                                                                            >  1:1201788/1‑148 (MQ=255)
                                                      ttCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCg                                                                                         >  1:1038176/1‑148 (MQ=255)
                                                           gAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCa                                                                                    >  1:2012682/1‑148 (MQ=255)
                                                             aTTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGt                                                                                  <  1:2966207/148‑1 (MQ=255)
                                                              ttCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGtt                                                                                 >  1:350238/1‑148 (MQ=255)
                                                               tCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGttt                                                                                >  1:2847377/1‑148 (MQ=255)
                                                                  aaCACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCag                                                                             >  1:1081014/1‑148 (MQ=255)
                                                                    cacGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCATGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCagag                                                                           >  2:3253995/1‑148 (MQ=255)
                                                                     acGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAgaga                                                                          >  1:2661339/1‑148 (MQ=255)
                                                                     acGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAgaga                                                                          >  1:880287/1‑148 (MQ=255)
                                                                      cGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGaa                                                                         >  2:865531/1‑148 (MQ=255)
                                                                       gCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAAc                                                                        >  1:717375/1‑148 (MQ=255)
                                                                       gCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAAc                                                                        >  1:1171504/1‑148 (MQ=255)
                                                                         cAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGAGATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTc                                                                      <  2:3159645/148‑1 (MQ=255)
                                                                          aaCCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCg                                                                     >  2:2067294/1‑148 (MQ=255)
                                                                          aaCCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCg                                                                     >  2:2067381/1‑148 (MQ=255)
                                                                           aCCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCgg                                                                    <  2:373388/148‑1 (MQ=255)
                                                                              ggACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTg                                                                 <  1:1745245/148‑1 (MQ=255)
                                                                               gACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGc                                                                <  2:1507524/148‑1 (MQ=255)
                                                                               gACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGc                                                                <  2:1739502/148‑1 (MQ=255)
                                                                               gACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGAGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGATAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGc                                                                >  1:101596/1‑148 (MQ=255)
                                                                                aCCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCt                                                               >  1:1338327/1‑148 (MQ=255)
                                                                                aCCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCt                                                               >  1:3261627/1‑148 (MQ=255)
                                                                                    ttCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTAcgc                                                           >  2:1617076/1‑148 (MQ=255)
                                                                                     tCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTAcgcg                                                          >  1:1247670/1‑148 (MQ=255)
                                                                                     tCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTAcgcg                                                          >  1:1045475/1‑148 (MQ=255)
                                                                                       cacaGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGaa                                                        >  1:3185935/1‑148 (MQ=255)
                                                                                            aaTGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCa                                                   >  2:1632116/1‑148 (MQ=255)
                                                                                             aTGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAg                                                  <  1:1843943/148‑1 (MQ=255)
                                                                                             aTGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAg                                                  >  2:1444465/1‑148 (MQ=255)
                                                                                                      gCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCg                                         <  1:2946570/148‑1 (MQ=255)
                                                                                                       cATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCgg                                        >  2:1035916/1‑148 (MQ=255)
                                                                                                         ttcttcAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGa                                      <  1:2928232/148‑1 (MQ=255)
                                                                                                            ttcAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGt                                   <  1:2267574/148‑1 (MQ=255)
                                                                                                             tcAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGtt                                  >  2:131420/1‑148 (MQ=255)
                                                                                                                ggATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCa                               >  1:2767512/1‑148 (MQ=255)
                                                                                                                ggATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCa                               >  1:1971968/1‑148 (MQ=255)
                                                                                                                 gATTTCAGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAg                              >  1:2285184/1‑148 (MQ=255)
                                                                                                                       aGCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTc                        >  2:3021491/1‑148 (MQ=255)
                                                                                                                        gCCAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTca                       >  2:1127149/1‑148 (MQ=255)
                                                                                                                         ccAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTcag                      >  2:1082900/1‑148 (MQ=255)
                                                                                                                          cAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTcagc                     >  2:1493581/1‑148 (MQ=255)
                                                                                                                          cAGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTcagc                     >  2:1179277/1‑148 (MQ=255)
                                                                                                                           aGTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTcagca                    <  1:571359/148‑1 (MQ=255)
                                                                                                                            gTTCCTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTcagcag                   <  2:959021/148‑1 (MQ=255)
                                                                                                                                cTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGt               <  2:2618021/148‑1 (MQ=255)
                                                                                                                                cTGAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGt               <  1:1691949/148‑1 (MQ=255)
                                                                                                                                  gAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGttt             <  2:5105/148‑1 (MQ=255)
                                                                                                                                  gAGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGttt             >  1:2228422/1‑148 (MQ=255)
                                                                                                                                   aGCGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTa            <  2:1133382/148‑1 (MQ=255)
                                                                                                                                     cGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGc          <  1:984606/148‑1 (MQ=255)
                                                                                                                                     cGGAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTCCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGc          <  1:691725/148‑1 (MQ=255)
                                                                                                                                       gAAGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCtt        >  1:3418075/1‑148 (MQ=255)
                                                                                                                                         aGCATCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCg      >  2:2335217/1‑148 (MQ=255)
                                                                                                                                            aTCACTGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCa   >  2:2115707/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AACGCTTCGCTGTCTGCGATAGCGTTAGCCACGATTTCAGTACCTTTGCGGAAGTTCGGGAATTCGAACACGCCAACCGGACCGTTCCACAGAATGGTTTTCGCATTCTTCAGGATTTCAGCCAGTTCCTGAGCGGAAGCATCACCGATATCCAGGATCTGCTCGTCAGCTTTCACATCGTTAACAGATTTCAGGGTAGCCGGTGCAGTTTCAGAGAACTCGGTTGCTACGCGAACATCAGACGGAACCGGGATGTTGCAGGTGGTCAGCAGACGTTTAGCTTCGTCAA  >  CP009273/3064978‑3065266

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: