Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,093,659 C→T 19.1% 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,6590CT19.1% 444.2 / 83.2 173C364Y (TGC→TAC) gnd6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base C (63/77);  new base T (14/19);  total (77/96)
Fisher's exact test for biased strand distribution p-value = 8.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 181 total aligned reads displayed.

GTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAG  >  CP009273/2093513‑2093806
                                                                                                                                                  |                                                                                                                                                   
gTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGca                                                                                                                                                    >  2:764235/1‑148 (MQ=255)
gTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGca                                                                                                                                                    >  2:2622185/1‑148 (MQ=255)
 tACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcag                                                                                                                                                   >  1:261223/1‑148 (MQ=255)
  aCTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGcagc                                                                                                                                                  <  2:248936/148‑1 (MQ=255)
   cTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGtagcc                                                                                                                                                 <  2:1432711/148‑1 (MQ=255)
          aGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACgg                                                                                                                                          >  2:1815/1‑148 (MQ=255)
           gCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGa                                                                                                                                         >  2:1444075/1‑148 (MQ=255)
            cAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGaa                                                                                                                                        >  1:3200861/1‑148 (MQ=255)
            cAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGaa                                                                                                                                        >  1:3198626/1‑148 (MQ=255)
                 aCATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCCGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCt                                                                                                                                   <  2:2252893/148‑1 (MQ=255)
                  cATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCtt                                                                                                                                  <  1:310556/148‑1 (MQ=255)
                   aTCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTc                                                                                                                                 >  1:1448550/1‑148 (MQ=255)
                      aCGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcgcg                                                                                                                              <  1:1255892/148‑1 (MQ=255)
                      aCGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTcgcg                                                                                                                              <  1:2456349/148‑1 (MQ=255)
                              ccTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGcc                                                                                                                      >  1:1786861/1‑148 (MQ=255)
                              ccTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGcc                                                                                                                      <  1:960940/148‑1 (MQ=255)
                               ctgctgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCg                                                                                                                     <  2:2315762/148‑1 (MQ=255)
                                 gctgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTa                                                                                                                   <  1:2829860/148‑1 (MQ=255)
                                 gctgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTa                                                                                                                   >  2:2011788/1‑148 (MQ=255)
                                   tgGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGt                                                                                                                 <  1:2566947/148‑1 (MQ=255)
                                    ggTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGtt                                                                                                                <  1:1949711/148‑1 (MQ=255)
                                      tAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCa                                                                                                              >  1:1053634/1‑148 (MQ=255)
                                        gTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGa                                                                                                            >  1:3015866/1‑148 (MQ=255)
                                        gTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGa                                                                                                            >  1:1509447/1‑148 (MQ=255)
                                         tCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGAt                                                                                                           >  1:2501199/1‑148 (MQ=255)
                                           aTCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATcc                                                                                                         >  1:754850/1‑148 (MQ=255)
                                            tCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATccc                                                                                                        <  2:724410/148‑1 (MQ=255)
                                             cGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCa                                                                                                       >  1:1058819/1‑148 (MQ=255)
                                              ggCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAg                                                                                                      <  2:2085317/148‑1 (MQ=255)
                                                cAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGtt                                                                                                    >  2:3292434/1‑148 (MQ=255)
                                                 aaTTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTg                                                                                                   <  2:2390390/148‑1 (MQ=255)
                                                   tttGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTa                                                                                                 <  1:2306925/148‑1 (MQ=255)
                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                                                                              >  2:1173079/1‑148 (MQ=255)
                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                                                                              <  1:2622185/148‑1 (MQ=255)
                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                                                                              <  1:764235/148‑1 (MQ=255)
                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                                                                              <  1:2465021/148‑1 (MQ=255)
                                                      gCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActc                                                                                              <  1:1448472/148‑1 (MQ=255)
                                                       cTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActct                                                                                             >  1:1691212/1‑148 (MQ=255)
                                                       cTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTActct                                                                                             <  1:1377131/148‑1 (MQ=255)
                                                         tGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTc                                                                                           >  1:1868235/1‑148 (MQ=255)
                                                            aGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCTCGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGa                                                                                        <  2:457116/148‑1 (MQ=255)
                                                            aGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGa                                                                                        <  2:2790181/148‑1 (MQ=255)
                                                             gTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGac                                                                                       <  1:3248770/148‑1 (MQ=255)
                                                              tACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacg                                                                                      >  2:345492/1‑148 (MQ=255)
                                                              tACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacg                                                                                      >  1:2377794/1‑148 (MQ=255)
                                                                      cAGCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgc                                                                              <  1:1670962/148‑1 (MQ=255)
                                                                        gCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcag                                                                            <  1:214563/148‑1 (MQ=255)
                                                                        gCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcag                                                                            <  1:1820365/148‑1 (MQ=255)
                                                                        gCAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcag                                                                            <  2:19237/148‑1 (MQ=255)
                                                                         cAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcagc                                                                           >  2:1744951/1‑148 (MQ=255)
                                                                         cAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacg                                                                               <  2:636828/144‑1 (MQ=255)
                                                                         cAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacg                                                                               <  2:2422439/144‑1 (MQ=255)
                                                                         cAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacg                                                                               >  1:636828/1‑144 (MQ=255)
                                                                         cAACAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGacgcagcacgcagc                                                                           <  2:1114263/148‑1 (MQ=255)
                                                                            cAGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTga                                                                        >  1:46421/1‑148 (MQ=255)
                                                                             aGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgag                                                                       <  2:674099/148‑1 (MQ=255)
                                                                             aGGTTAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgag                                                                       <  1:332884/148‑1 (MQ=255)
                                                                               gTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTgagag                                                                     <  2:1580376/148‑1 (MQ=255)
                                                                                 tAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGaa                                                                   >  1:512220/1‑148 (MQ=255)
                                                                                 tAGCGATCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGaa                                                                   <  1:653800/148‑1 (MQ=255)
                                                                                  aGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAg                                                                  <  1:504781/148‑1 (MQ=255)
                                                                                   gCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGc                                                                 >  2:1004470/1‑148 (MQ=255)
                                                                                    cGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGcc                                                                <  1:3078822/148‑1 (MQ=255)
                                                                                     gATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGccc                                                               <  1:3385959/148‑1 (MQ=255)
                                                                                       tCTGTGGATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTg                                                             >  1:174366/1‑148 (MQ=255)
                                                                                            ggATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGt                                                        <  2:824348/148‑1 (MQ=255)
                                                                                             gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTa                                                       <  2:2085409/148‑1 (MQ=255)
                                                                                             gATTTTAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTa                                                       <  2:2085594/148‑1 (MQ=255)
                                                                                                tttAGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGa                                                    <  1:981710/148‑1 (MQ=255)
                                                                                                 ttaggCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaa                                                   <  2:273559/145‑1 (MQ=255)
                                                                                                 ttaggCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaa                                                   <  1:3340562/145‑1 (MQ=255)
                                                                                                  taggCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaaa                                                  >  2:3145092/3‑148 (MQ=255)
                                                                                                    ggCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACg                                                <  2:1203608/148‑1 (MQ=255)
                                                                                                     gCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGa                                               <  2:311471/148‑1 (MQ=255)
                                                                                                     gCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGa                                               >  1:110137/1‑148 (MQ=255)
                                                                                                      cATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAt                                              <  1:2262429/148‑1 (MQ=255)
                                                                                                       aTAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAtt                                             <  2:1691212/148‑1 (MQ=255)
                                                                                                       aTAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAtt                                             >  2:892656/1‑148 (MQ=255)
                                                                                                       aTAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGAtt                                             <  2:433846/148‑1 (MQ=255)
                                                                                                          aGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTg                                          <  1:3021407/148‑1 (MQ=255)
                                                                                                          aGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTg                                          <  2:2784381/148‑1 (MQ=255)
                                                                                                           gCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGc                                         >  1:1745433/1‑148 (MQ=255)
                                                                                                           gCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGc                                         >  2:2891197/1‑148 (MQ=255)
                                                                                                            cATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGcc                                        >  2:2700033/1‑148 (MQ=255)
                                                                                                              tCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGAGTAAGAAACGATTTTGcaca                                      >  2:953025/1‑148 (MQ=255)
                                                                                                               cGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAg                                     <  2:1000057/148‑1 (MQ=255)
                                                                                                                 gTGATTTTCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGAt                                   >  2:3342005/1‑148 (MQ=255)
                                                                                                                 gTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGAt                                   >  1:2024878/1‑148 (MQ=255)
                                                                                                                 gTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGAt                                   <  2:1890741/148‑1 (MQ=255)
                                                                                                                  tGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATa                                  >  1:1832562/1‑148 (MQ=255)
                                                                                                                      tttCTGCATGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAgc                              >  2:2764082/1‑148 (MQ=255)
                                                                                                                      tttCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAgc                              >  2:2045953/1‑148 (MQ=255)
                                                                                                                        tCTGCAGGAACTGCGCACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAgcgc                            >  2:3091261/1‑148 (MQ=255)
                                                                                                                         cTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCa                           <  1:2011736/148‑1 (MQ=255)
                                                                                                                                 aCTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAAc                   >  2:976238/1‑148 (MQ=255)
                                                                                                                                   tGCGCACGGATGATGCAGCCAGCACGTAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCTTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACtt                 <  1:1815/148‑1 (MQ=255)
                                                                                                                                   tGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACtt                 >  1:859192/1‑148 (MQ=255)
                                                                                                                                      gcACGGATGATGTAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACttttt              <  2:2353329/148‑1 (MQ=255)
                                                                                                                                               atgCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACt     <  1:698904/148‑1 (MQ=255)
                                                                                                                                                tgCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTc    <  2:3200823/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTACTGCATAAGCAACGACATCACGCAGCGCCTGCTGGTAGTCATCGGCAATTTGCTTGAAGTACGGAGCCAGCAACAGGTTAGCGATCTGTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAG  >  CP009273/2093513‑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: ‑