Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 137,915 G→T 100% intergenic (‑203/‑3) gcd ← / → hpt glucose dehydrogenase/hypoxanthine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273137,9150GT97.7% 542.7 / 2.0 172intergenic (‑203/‑3)gcd/hptglucose dehydrogenase/hypoxanthine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (89/79);  total (92/80)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 177 total aligned reads displayed.

ATGTGTGATCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGAGATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCATTTATGTTT  >  CP009273/137771‑138061
                                                                                                                                                |                                                                                                                                                  
aTGTGTGATCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGAtata                                                                                                                                                 <  1:1355643/148‑1 (MQ=255)
 tgtgtgATCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGatatat                                                                                                                                                <  2:71645/148‑1 (MQ=255)
    gtgATCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGaa                                                                                                                                             >  1:1445127/1‑148 (MQ=255)
    gtgATCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGaa                                                                                                                                             >  1:720753/1‑148 (MQ=255)
     tgATCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGaaa                                                                                                                                            <  2:684413/148‑1 (MQ=255)
        tCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACAt                                                                                                                                         >  2:3115311/1‑148 (MQ=255)
        tCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACAt                                                                                                                                         <  2:1912970/148‑1 (MQ=255)
          gTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATAc                                                                                                                                       <  1:1371481/148‑1 (MQ=255)
           tcatcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACt                                                                                                                                      >  1:1752484/1‑148 (MQ=255)
             atcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGt                                                                                                                                    <  2:320752/148‑1 (MQ=255)
             atcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGt                                                                                                                                    <  1:2234001/148‑1 (MQ=255)
             atcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGt                                                                                                                                    >  1:661852/1‑148 (MQ=255)
             atcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGt                                                                                                                                    <  1:1045220/148‑1 (MQ=255)
             atcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGt                                                                                                                                    <  1:1045429/148‑1 (MQ=255)
              tcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTa                                                                                                                                   >  2:485165/1‑148 (MQ=255)
              tcaCAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTa                                                                                                                                   >  1:2497479/1‑148 (MQ=255)
               cacaATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAg                                                                                                                                  <  2:884676/148‑1 (MQ=255)
               cacaATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAg                                                                                                                                  >  1:493261/1‑148 (MQ=255)
                acaATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGa                                                                                                                                 >  1:204871/1‑148 (MQ=255)
                acaATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGAGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATTAAACATACTGTAGa                                                                                                                                 >  1:2000297/1‑148 (MQ=255)
                    ttCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTa                                                                                                                             <  1:3160709/148‑1 (MQ=255)
                     tCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTaa                                                                                                                            >  1:189847/1‑148 (MQ=255)
                       gAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATg                                                                                                                          <  2:3177915/148‑1 (MQ=255)
                       gAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATg                                                                                                                          <  2:3177910/148‑1 (MQ=255)
                        aGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGa                                                                                                                         <  2:1761925/148‑1 (MQ=255)
                        aGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGa                                                                                                                         <  2:3475764/148‑1 (MQ=255)
                          cTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATc                                                                                                                       <  1:920765/148‑1 (MQ=255)
                           tttATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATcc                                                                                                                      <  2:1476840/148‑1 (MQ=255)
                            ttattaACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATccc                                                                                                                     >  2:1739078/1‑148 (MQ=255)
                                taACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGaa                                                                                                                 >  1:3082040/1‑148 (MQ=255)
                                   cAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCg                                                                                                              >  2:1384665/1‑148 (MQ=255)
                                     gATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGga                                                                                                            >  2:17334/1‑148 (MQ=255)
                                      aTTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGgag                                                                                                           >  2:653482/1‑148 (MQ=255)
                                        tCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGAt                                                                                                         >  2:810899/1‑148 (MQ=255)
                                        tCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGAt                                                                                                         >  2:810922/1‑148 (MQ=255)
                                        tCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTACCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGAGATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGAt                                                                                                         <  2:3086857/148‑1 (MQ=255)
                                            cgAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTaaa                                                                                                     <  2:1536305/148‑1 (MQ=255)
                                            cgAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTaaa                                                                                                     <  2:2031588/148‑1 (MQ=255)
                                            cgAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTaaa                                                                                                     <  2:2607585/148‑1 (MQ=255)
                                            cgAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTaaa                                                                                                     <  2:2608331/148‑1 (MQ=255)
                                              aatgaatAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAgc                                                                                                   >  2:2640056/1‑148 (MQ=255)
                                              aatgaatAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAgc                                                                                                   >  2:1187729/1‑148 (MQ=255)
                                                tgaatAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAgcgc                                                                                                 >  1:2927075/1‑148 (MQ=255)
                                                    tAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTAt                                                                                             <  2:2325463/148‑1 (MQ=255)
                                                    tAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTAt                                                                                             <  2:1250328/148‑1 (MQ=255)
                                                     aGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATc                                                                                            >  2:1616340/1‑148 (MQ=255)
                                                      gTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCg                                                                                           <  1:158852/148‑1 (MQ=255)
                                                       ttttACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGc                                                                                          >  2:2771245/1‑148 (MQ=255)
                                                        tttACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGcc                                                                                         >  1:644259/1‑148 (MQ=255)
                                                               gTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTgg                                                                                  <  1:1410826/148‑1 (MQ=255)
                                                                tataCTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTggg                                                                                 <  2:3323050/148‑1 (MQ=255)
                                                                  taCTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGgtc                                                                               >  1:828403/1‑148 (MQ=255)
                                                                    cTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGgtcgt                                                                             >  2:3100309/1‑148 (MQ=255)
                                                                     tGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGgtcgtc                                                                            <  1:1103998/148‑1 (MQ=255)
                                                                     tGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGgtcgtc                                                                            <  2:1160661/148‑1 (MQ=255)
                                                                            cTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTAc                                                                     >  1:2274762/1‑148 (MQ=255)
                                                                               gcgcTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGAGATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGa                                                                  >  1:2908386/1‑148 (MQ=255)
                                                                                 gcTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGc                                                                >  1:1603766/1‑148 (MQ=255)
                                                                                  cTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCg                                                               >  2:1382303/1‑148 (MQ=255)
                                                                                         gCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACaaa                                                        >  1:151329/1‑148 (MQ=255)
                                                                                         gCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACaaa                                                        <  1:1639973/148‑1 (MQ=255)
                                                                                          cGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAg                                                       <  1:1508405/148‑1 (MQ=255)
                                                                                           ggTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGa                                                      <  2:2147519/148‑1 (MQ=255)
                                                                                            gTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGAc                                                     <  2:2910089/148‑1 (MQ=255)
                                                                                              gCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAg                                                   <  2:2179540/148‑1 (MQ=255)
                                                                                              gCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAg                                                   <  2:2179568/148‑1 (MQ=255)
                                                                                              gCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAg                                                   <  1:1707629/148‑1 (MQ=255)
                                                                                               ccAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGc                                                  <  2:262148/148‑1 (MQ=255)
                                                                                                   aaCCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCa                                              >  2:2090812/1‑148 (MQ=255)
                                                                                                       tGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATTATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGa                                          <  1:2610774/148‑1 (MQ=255)
                                                                                                       tGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGa                                          >  2:396942/1‑148 (MQ=255)
                                                                                                        gCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGat                                         >  1:1791870/1‑148 (MQ=255)
                                                                                                         cccGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGata                                        <  2:1939207/148‑1 (MQ=255)
                                                                                                          ccGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGatat                                       <  2:3047354/148‑1 (MQ=255)
                                                                                                           cGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATg                                      <  2:3448430/148‑1 (MQ=255)
                                                                                                            gTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATgg                                     >  1:3116092/1‑148 (MQ=255)
                                                                                                             tGCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGt                                    >  1:57834/1‑148 (MQ=255)
                                                                                                              gCGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTg                                   >  2:1066793/1‑148 (MQ=255)
                                                                                                               cGAAGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGc                                  <  1:3412163/148‑1 (MQ=255)
                                                                                                                 aaGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCtg                                >  1:298150/1‑148 (MQ=255)
                                                                                                                  aGTGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCtgg                               >  1:979194/1‑148 (MQ=255)
                                                                                                                    tGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCtggtg                             <  2:2484442/148‑1 (MQ=255)
                                                                                                                    tGATTTGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCtggtg                             <  2:1770566/148‑1 (MQ=255)
                                                                                                                        ttGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTc                         >  2:3082655/1‑148 (MQ=255)
                                                                                                                         tGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTct                        <  1:3293421/148‑1 (MQ=255)
                                                                                                                         tGTTTTTAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTct                        <  1:2839386/148‑1 (MQ=255)
                                                                                                                            ttttAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGAGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTctgct                     >  2:3057479/1‑148 (MQ=255)
                                                                                                                             tttAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTctgctg                    <  2:481138/148‑1 (MQ=255)
                                                                                                                              ttAAATCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCtgctgc                   >  2:566144/1‑148 (MQ=255)
                                                                                                                                 aaTCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTg                >  1:1183522/1‑148 (MQ=255)
                                                                                                                                 aaTCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTg                >  1:2029598/1‑148 (MQ=255)
                                                                                                                                  aTCATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTgg               >  2:2855790/1‑148 (MQ=255)
                                                                                                                                    cATATGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCt             >  2:67894/1‑148 (MQ=255)
                                                                                                                                       atGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGGGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCAt          >  1:2313092/1‑148 (MQ=255)
                                                                                                                                       atGGTTAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCAt          <  1:2713552/148‑1 (MQ=255)
                                                                                                                                       atGGTTAGAGATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCAt          >  1:3067725/1‑148 (MQ=255)
                                                                                                                                           ttAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCAtttat      >  2:188439/1‑148 (MQ=255)
                                                                                                                                            tAGATATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCAtttatg     >  2:794215/1‑148 (MQ=255)
                                                                                                                                             agatATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCAtttatgt    >  2:12202/1‑148 (MQ=255)
                                                                                                                                               agaTATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCAtttatgttt  <  2:1013608/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ATGTGTGATCGTCATCACAATTCGAGCTTTATTAACAGATTCCGCGAATGAATAGTTTTACTGGTATACTGCGTGTCTTGCGCTTTGTTGCGGTGCCAAAACCTGCCCGTGCGAAGTGATTTGTTTTTAAATCATATGGTTAGAGATATGAAACATACTGTAGAAGTAATGATCCCCGAAGCGGAGATTAAAGCGCGTATCGCCGAACTGGGTCGTCAGATTACTGAGCGTTACAAAGACAGCGGCAGCGATATGGTGCTGGTGGGTCTGCTGCGTGGCTCATTTATGTTT  >  CP009273/137771‑138061

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