Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,093,659 C→T 20.0% C364Y (TGC→TAC)  gnd ← 6‑phosphogluconate dehydrogenase, decarboxylating

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,093,6590CT20.0% 522.1 / 111.8 210C364Y (TGC→TAC) gnd6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base C (85/83);  new base T (19/23);  total (104/106)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

TGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAG  >  CP009273/2093512‑2093806
                                                                                                                                                   |                                                                                                                                                   
tGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGc                                                                                                                                                     <  1:3732430/148‑1 (MQ=255)
tGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGc                                                                                                                                                     <  2:3440569/148‑1 (MQ=255)
   aCTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagc                                                                                                                                                  <  2:644956/148‑1 (MQ=255)
      gCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagccag                                                                                                                                               >  2:1016998/1‑148 (MQ=255)
       cATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagccagc                                                                                                                                              >  1:3538533/1‑148 (MQ=255)
         tAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCAc                                                                                                                                            >  2:763790/1‑148 (MQ=255)
          aaGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACg                                                                                                                                           >  1:1519277/1‑148 (MQ=255)
               aCGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGaaaa                                                                                                                                      <  1:1079235/148‑1 (MQ=255)
                cGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAAt                                                                                                                                     >  2:542786/1‑148 (MQ=255)
                  aCATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCt                                                                                                                                   <  2:3808526/148‑1 (MQ=255)
                  aCATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCt                                                                                                                                   >  2:825680/1‑148 (MQ=255)
                  aCATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCt                                                                                                                                   <  1:534366/148‑1 (MQ=255)
                   cATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCtt                                                                                                                                  >  1:426392/1‑148 (MQ=255)
                      cACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcgc                                                                                                                               >  2:1851476/1‑148 (MQ=255)
                      cACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcgc                                                                                                                               >  2:1351881/1‑148 (MQ=255)
                       aCGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcgcg                                                                                                                              <  1:2448631/148‑1 (MQ=255)
                        cGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGa                                                                                                                             <  1:365653/148‑1 (MQ=255)
                         gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAt                                                                                                                            <  2:736977/148‑1 (MQ=255)
                         gCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAt                                                                                                                            <  1:680473/148‑1 (MQ=255)
                          cAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGAtt                                                                                                                           <  1:1954479/148‑1 (MQ=255)
                           aGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGAttt                                                                                                                          <  1:1810452/148‑1 (MQ=255)
                                cttctgGTATTCAACGGCAATTTGCTTGAATTAAGGAGCCAGCAACAGGTTAGCGATATTTGGATTTTAGGCATAAGCATAGGTGATTTTATTCACGAACTGCGCACGGATGATGTAGCAAGCACCGAAAATCTTCGCGATTTAGCCg                                                                                                                     >  2:2328802/4‑148 (MQ=255)
                                    tgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGt                                                                                                                 <  1:126916/148‑1 (MQ=255)
                                     ggTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGtt                                                                                                                <  1:36095/148‑1 (MQ=255)
                                          tCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGAt                                                                                                           >  1:2343822/1‑148 (MQ=255)
                                           cATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATc                                                                                                          >  2:3252170/1‑148 (MQ=255)
                                            aTCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATcc                                                                                                         <  2:551949/148‑1 (MQ=255)
                                            aTCGGCAATGTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATcc                                                                                                         <  2:2659537/148‑1 (MQ=255)
                                             tCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATccc                                                                                                        <  2:2251211/148‑1 (MQ=255)
                                                gCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGt                                                                                                     >  2:1435916/1‑148 (MQ=255)
                                                    tttGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTa                                                                                                 <  1:1851476/148‑1 (MQ=255)
                                                    tttGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTa                                                                                                 >  2:481357/1‑148 (MQ=255)
                                                     ttGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTAGCCGTAGTTCAGATCCCAGTTGTAc                                                                                                >  1:3309434/1‑148 (MQ=255)
                                                     ttGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTAGCCGTAGTTCAGATCCCAGTTGTAc                                                                                                >  1:2892409/1‑148 (MQ=255)
                                                         ttGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCtt                                                                                            <  1:951023/148‑1 (MQ=255)
                                                          tGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTc                                                                                           <  2:3014431/148‑1 (MQ=255)
                                                           gAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCa                                                                                          <  1:20484/148‑1 (MQ=255)
                                                           gAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCa                                                                                          <  1:1952768/148‑1 (MQ=255)
                                                           gAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCa                                                                                          <  1:3747468/148‑1 (MQ=255)
                                                             aGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGa                                                                                        <  1:3587424/148‑1 (MQ=255)
                                                                aCGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgc                                                                                     <  1:769580/148‑1 (MQ=255)
                                                                    aGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTAGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagca                                                                                 >  2:2537199/1‑148 (MQ=255)
                                                                            aCAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTg                                                                         >  1:1916694/1‑148 (MQ=255)
                                                                            aCAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTg                                                                         <  2:3751454/148‑1 (MQ=255)
                                                                            aCAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTg                                                                         <  1:1704618/148‑1 (MQ=255)
                                                                              aGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgag                                                                       <  2:3678965/148‑1 (MQ=255)
                                                                                    gCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGc                                                                 <  2:1519287/148‑1 (MQ=255)
                                                                                       aTCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCt                                                              >  2:2958076/1‑148 (MQ=255)
                                                                                        tCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGTAGCCCTg                                                             <  1:2413919/148‑1 (MQ=255)
                                                                                        tCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTg                                                             >  1:3299700/1‑148 (MQ=255)
                                                                                           gtgGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGc                                                          >  1:1302635/1‑148 (MQ=255)
                                                                                            tgGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCg                                                         <  2:2377913/148‑1 (MQ=255)
                                                                                            tgGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCg                                                         <  2:455400/148‑1 (MQ=255)
                                                                                              gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTa                                                       <  1:2531763/148‑1 (MQ=255)
                                                                                              gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTa                                                       <  1:2531701/148‑1 (MQ=255)
                                                                                               atttgaggCATTAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTaa                                                      <  2:3595771/142‑1 (MQ=255)
                                                                                               aTTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTaa                                                      >  1:1378919/1‑148 (MQ=255)
                                                                                                 tttCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGa                                                    <  2:1141255/148‑1 (MQ=255)
                                                                                                 tttCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGa                                                    <  2:1822817/148‑1 (MQ=255)
                                                                                                 tttAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGa                                                    >  1:565747/1‑148 (MQ=255)
                                                                                                 tttAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGa                                                    <  1:1145716/148‑1 (MQ=255)
                                                                                                   taggCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaaa                                                  <  2:3137629/146‑1 (MQ=255)
                                                                                                      gCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGa                                               >  1:2135463/1‑148 (MQ=255)
                                                                                                      gCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTAGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGa                                               <  1:3277253/148‑1 (MQ=255)
                                                                                                       cATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAt                                              <  1:3404723/148‑1 (MQ=255)
                                                                                                         tAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAttt                                            >  2:3081609/1‑148 (MQ=255)
                                                                                                         tAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAttt                                            >  1:905776/1‑148 (MQ=255)
                                                                                                          aaGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTGAGAAACGAtttt                                           >  2:643940/1‑148 (MQ=255)
                                                                                                          aaGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAtttt                                           >  1:3670408/1‑148 (MQ=255)
                                                                                                          aaGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAtttt                                           >  2:3044465/1‑148 (MQ=255)
                                                                                                           aGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTg                                          <  2:1629237/148‑1 (MQ=255)
                                                                                                           aGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTg                                          >  1:706057/1‑148 (MQ=255)
                                                                                                             cATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGcc                                        >  1:2115134/1‑148 (MQ=255)
                                                                                                                cGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCGAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCTGACGCAGCACGCAGCTGAGAGATGCCCTGGGCGTAAGAAACGATTTTGCCCAg                                     <  1:1016996/148‑1 (MQ=255)
                                                                                                                cGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAg                                     <  2:3538533/148‑1 (MQ=255)
                                                                                                                   tGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATa                                  >  2:3267622/1‑148 (MQ=255)
                                                                                                                     aTTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACa                                >  2:2826824/1‑148 (MQ=255)
                                                                                                                       tttCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAgc                              >  2:3490093/1‑148 (MQ=255)
                                                                                                                            tcAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCacg                         <  2:3496579/147‑1 (MQ=255)
                                                                                                                            gCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCacg                         <  1:2265357/148‑1 (MQ=255)
                                                                                                                            gCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCacg                         >  1:3078895/1‑148 (MQ=255)
                                                                                                                            gCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCacg                         >  1:3093642/1‑148 (MQ=255)
                                                                                                                            gCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCacg                         >  1:844539/1‑148 (MQ=255)
                                                                                                                             cAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCacga                        >  2:734525/1‑148 (MQ=255)
                                                                                                                               ggAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCacgacg                      >  2:1886245/1‑148 (MQ=255)
                                                                                                                                 aaCTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGaa                    <  2:3352876/148‑1 (MQ=255)
                                                                                                                                 aaCTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGaa                    <  1:2959264/148‑1 (MQ=255)
                                                                                                                                   cTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACt                  <  2:1493753/148‑1 (MQ=255)
                                                                                                                                    tGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACtt                 <  1:183918/148‑1 (MQ=255)
                                                                                                                                       gcACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACttttt              <  2:2873382/148‑1 (MQ=255)
                                                                                                                                       gcACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACttttt              <  1:2931778/148‑1 (MQ=255)
                                                                                                                                       gcACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACttttt              <  1:3452757/148‑1 (MQ=255)
                                                                                                                                         aCGGATGATGCAGCCAGCACTGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCg            <  2:1275410/148‑1 (MQ=255)
                                                                                                                                         aCGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCg            <  1:3376950/148‑1 (MQ=255)
                                                                                                                                          cGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGa           >  2:3402687/1‑148 (MQ=255)
                                                                                                                                             atgatgTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGa        >  2:1116347/1‑148 (MQ=255)
                                                                                                                                             atgatgTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGa        <  2:1936304/148‑1 (MQ=255)
                                                                                                                                              tgatgTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGaa       <  2:3361904/148‑1 (MQ=255)
                                                                                                                                              tgatgCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGaa       <  1:2366913/148‑1 (MQ=255)
                                                                                                                                                   cagccagcACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAg  >  1:3709312/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAG  >  CP009273/2093512‑2093806

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