Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,093,612 C→A 29.1% E380* (GAA→TAA)  gnd ← 6‑phosphogluconate dehydrogenase, decarboxylating

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,093,6120CA29.1% 443.6 / 152.1 219E380* (GAA→TAA) gnd6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base C (85/70);  new base A (33/31);  total (118/101)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

TCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAAC  >  CP009273/2093464‑2093759
                                                                                                                                                    |                                                                                                                                                   
     aTAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCAt                                                                                                                                                >  2:3834702/1‑149 (MQ=255)
       aGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATaa                                                                                                                                              <  1:845139/149‑1 (MQ=255)
            aCCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATc                                                                                                                                         >  1:2605435/1‑149 (MQ=255)
            aCCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATc                                                                                                                                         <  1:3879195/149‑1 (MQ=255)
             ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCg                                                                                                                                        >  2:2778361/1‑149 (MQ=255)
             ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCg                                                                                                                                        >  1:1451154/1‑149 (MQ=255)
             ccGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCg                                                                                                                                        <  2:2347004/149‑1 (MQ=255)
                 tGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGa                                                                                                                                    >  2:2281740/1‑149 (MQ=255)
                  gCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGAt                                                                                                                                   <  1:230785/149‑1 (MQ=255)
                   cGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGAtt                                                                                                                                  <  2:3082680/149‑1 (MQ=255)
                   cGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGAtt                                                                                                                                  <  2:1400801/149‑1 (MQ=255)
                      agaAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTc                                                                                                                               <  1:1025956/149‑1 (MQ=255)
                      agaAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTc                                                                                                                               <  2:1952509/149‑1 (MQ=255)
                      agaAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTc                                                                                                                               <  2:1321030/149‑1 (MQ=255)
                      agaAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTc                                                                                                                               <  1:2855428/149‑1 (MQ=255)
                          ggTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCa                                                                                                                           <  2:120116/149‑1 (MQ=255)
                          ggTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCa                                                                                                                           <  1:1944955/149‑1 (MQ=255)
                              ggAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGaa                                                                                                                       <  1:993225/149‑1 (MQ=255)
                               gAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAAc                                                                                                                      <  1:1003821/149‑1 (MQ=255)
                                 aCCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTg                                                                                                                    >  1:3800775/1‑149 (MQ=255)
                                  ccGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATGTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTgc                                                                                                                   <  1:2242171/149‑1 (MQ=255)
                                    ggAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTgcgc                                                                                                                 <  1:3043749/149‑1 (MQ=255)
                                    ggAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTgcgc                                                                                                                 <  1:2950076/149‑1 (MQ=255)
                                     gAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCa                                                                                                                <  1:3634339/149‑1 (MQ=255)
                                      aaTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCAc                                                                                                               <  1:1435956/149‑1 (MQ=255)
                                       aTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACg                                                                                                              >  2:429392/1‑149 (MQ=255)
                                       aTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACg                                                                                                              >  2:2527744/1‑149 (MQ=255)
                                       aTACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACg                                                                                                              >  2:3333583/1‑149 (MQ=255)
                                       aTACAGTTCTGTACTGCATAAGCAACGACATCACGCAGAGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGAGCACg                                                                                                              >  1:3355855/1‑149 (MQ=255)
                                         aCCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGga                                                                                                            <  1:1397099/149‑1 (MQ=255)
                                            gTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatga                                                                                                         >  1:2672628/1‑149 (MQ=255)
                                            gTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatga                                                                                                         <  2:1524884/149‑1 (MQ=255)
                                            gTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatga                                                                                                         >  2:3052880/1‑149 (MQ=255)
                                             ttCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGgatgat                                                                                                        <  1:1093247/149‑1 (MQ=255)
                                               cTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGc                                                                                                      >  2:2701676/1‑149 (MQ=255)
                                                tGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGatgatgta                                                                                                     >  2:3651358/1‑147 (MQ=255)
                                                tGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGca                                                                                                     >  2:1369562/1‑149 (MQ=255)
                                                      gCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGtagccagc                                                                                               <  2:92670/149‑1 (MQ=255)
                                                      gCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGAGGATGtagccagc                                                                                               >  1:3129597/1‑149 (MQ=255)
                                                        aTAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCAc                                                                                             >  1:94712/1‑149 (MQ=255)
                                                           aGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGa                                                                                          >  1:493665/1‑149 (MQ=255)
                                                           aGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGa                                                                                          >  2:3235734/1‑149 (MQ=255)
                                                              aaCGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGaaaa                                                                                       >  1:2862828/1‑149 (MQ=255)
                                                                 gACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCt                                                                                    >  1:2561857/1‑149 (MQ=255)
                                                                 gACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCt                                                                                    >  2:2473298/1‑149 (MQ=255)
                                                                    aTCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcg                                                                                 >  2:1276408/1‑149 (MQ=255)
                                                                     tCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcgc                                                                                >  1:4092040/1‑149 (MQ=255)
                                                                       aCGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGa                                                                              <  2:727112/149‑1 (MQ=255)
                                                                        cGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGAt                                                                             <  2:3436632/149‑1 (MQ=255)
                                                                        cGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAt                                                                             >  1:3724530/1‑149 (MQ=255)
                                                                         gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGAtt                                                                            >  2:1896318/1‑149 (MQ=255)
                                                                         gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGAtt                                                                            >  2:3888001/1‑149 (MQ=255)
                                                                         gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAtt                                                                            <  1:2701676/149‑1 (MQ=255)
                                                                          cAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAttt                                                                           <  1:3680907/149‑1 (MQ=255)
                                                                           aGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTc                                                                          >  1:881614/1‑149 (MQ=255)
                                                                           aGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTc                                                                          >  2:1688640/1‑149 (MQ=255)
                                                                                ctgctgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGt                                                                     >  1:2465068/1‑149 (MQ=255)
                                                                                ctgctgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGt                                                                     >  1:2120546/1‑149 (MQ=255)
                                                                                  gctgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCTCCGTAg                                                                   >  1:2443910/1‑149 (MQ=255)
                                                                                    tgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGAAGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGtt                                                                 >  2:2250534/1‑149 (MQ=255)
                                                                                     ggTAGTCTTCGGCAATTTGCTTGAAGTACGGAGCGAGAAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGTTGATTTTCTGCAGGAACTGCGCACGGATGATGAAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTc                                                                <  2:2899708/149‑1 (MQ=255)
                                                                                      gTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCa                                                               >  1:664965/1‑149 (MQ=255)
                                                                                      gTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCa                                                               >  2:765003/1‑149 (MQ=255)
                                                                                      gTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCa                                                               >  2:240006/1‑149 (MQ=255)
                                                                                         gTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGAt                                                            >  1:615751/1‑149 (MQ=255)
                                                                                           cATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATcc                                                          >  2:1290187/1‑149 (MQ=255)
                                                                                            aTCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATccc                                                         <  2:3486309/149‑1 (MQ=255)
                                                                                            aTCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATccc                                                         >  2:2918900/1‑149 (MQ=255)
                                                                                                   aTTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTa                                                  >  1:2723907/1‑149 (MQ=255)
                                                                                                   aTTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTa                                                  >  1:2723792/1‑149 (MQ=255)
                                                                                                      tGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                               <  2:1009276/149‑1 (MQ=255)
                                                                                                       gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActct                                              >  2:2830090/1‑149 (MQ=255)
                                                                                                         ttGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTc                                            <  1:727415/149‑1 (MQ=255)
                                                                                                           gAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAg                                          <  2:491856/149‑1 (MQ=255)
                                                                                                           gAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAg                                          <  2:1141910/149‑1 (MQ=255)
                                                                                                            aaGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGa                                         >  1:2589110/1‑149 (MQ=255)
                                                                                                              gTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacg                                       >  2:1331581/1‑149 (MQ=255)
                                                                                                              gTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacg                                       <  2:3728892/149‑1 (MQ=255)
                                                                                                              gTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacg                                       <  2:3724489/149‑1 (MQ=255)
                                                                                                               tACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgc                                      <  2:3392217/149‑1 (MQ=255)
                                                                                                                aCGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgca                                     >  2:1150659/1‑149 (MQ=255)
                                                                                                                aCGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgca                                     >  2:175370/1‑149 (MQ=255)
                                                                                                                   gAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagca                                  <  1:3938218/149‑1 (MQ=255)
                                                                                                                   gAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTAGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagca                                  >  1:204800/1‑149 (MQ=255)
                                                                                                                       cAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgca                              >  1:1797837/1‑149 (MQ=255)
                                                                                                                            aCAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTga                         <  1:3588205/149‑1 (MQ=255)
                                                                                                                               ggTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgagag                      <  1:1628159/149‑1 (MQ=255)
                                                                                                                               ggTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgagag                      <  1:3729075/149‑1 (MQ=255)
                                                                                                                               ggTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgagag                      <  1:346005/149‑1 (MQ=255)
                                                                                                                                   aGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGc                  <  1:395712/149‑1 (MQ=255)
                                                                                                                                    gCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGcc                 >  1:3981210/1‑149 (MQ=255)
                                                                                                                                      gATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCt               >  1:3245687/1‑149 (MQ=255)
                                                                                                                                       aTCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTg              <  1:1150682/149‑1 (MQ=255)
                                                                                                                                       aTCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTg              >  2:1794209/1‑149 (MQ=255)
                                                                                                                                       aTCTGTAGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTg              <  1:1150361/149‑1 (MQ=255)
                                                                                                                                            tgGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTAGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGt         <  1:2800680/149‑1 (MQ=255)
                                                                                                                                              gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTaa       >  1:1895131/1‑149 (MQ=255)
                                                                                                                                                 tttCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaa    >  1:1387998/1‑149 (MQ=255)
                                                                                                                                                  ttCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaaa   <  2:2103888/149‑1 (MQ=255)
                                                                                                                                                  ttCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaaa   >  2:3144605/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGTAATAGGCAACCGCTGCGGAGAAGGTCGGAACCGGAATACCGTTCTGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAAC  >  CP009273/2093464‑2093759

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: