Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,688,006 C→A 100% P299P (CCG→CCT purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,688,0060CA96.1% 513.0 / 12.0 194P299P (CCG→CCTpurLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (91/95);  total (95/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 194 total aligned reads displayed.

TGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACG  >  CP009273/2687863‑2688150
                                                                                                                                               |                                                                                                                                                
tGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggt                                                                                                                                              <  1:1199789/148‑1 (MQ=255)
 gTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtg                                                                                                                                             >  2:1019291/1‑148 (MQ=255)
  ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                            >  1:475655/1‑148 (MQ=255)
   tCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtggt                                                                                                                                           <  1:313459/148‑1 (MQ=255)
   tCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtggt                                                                                                                                           >  2:1721280/1‑148 (MQ=255)
   tCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtggt                                                                                                                                           >  2:1721473/1‑148 (MQ=255)
     gAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtg                                                                                                                                         <  2:1206994/148‑1 (MQ=255)
      aaGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtgt                                                                                                                                        <  2:3556140/148‑1 (MQ=255)
       aGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtgtg                                                                                                                                       >  1:2068700/1‑148 (MQ=255)
       aGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtgtg                                                                                                                                       >  1:2068373/1‑148 (MQ=255)
        gCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGa                                                                                                                                      >  1:2281711/1‑148 (MQ=255)
          cAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGt                                                                                                                                    >  1:112847/1‑148 (MQ=255)
          cAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGt                                                                                                                                    <  1:675779/148‑1 (MQ=255)
           aGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGtt                                                                                                                                   >  2:2589222/1‑148 (MQ=255)
           aGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGtt                                                                                                                                   <  2:1900995/148‑1 (MQ=255)
            ggAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGttt                                                                                                                                  >  1:1779685/1‑148 (MQ=255)
              aaTTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCg                                                                                                                                >  2:1579510/1‑148 (MQ=255)
              aaTTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCg                                                                                                                                >  2:1281373/1‑148 (MQ=255)
               aTTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGa                                                                                                                               <  2:2646303/148‑1 (MQ=255)
               aTTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGa                                                                                                                               <  1:672071/148‑1 (MQ=255)
                ttCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGAc                                                                                                                              >  1:2226641/1‑148 (MQ=255)
                 tCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACt                                                                                                                             >  1:1936897/1‑148 (MQ=255)
                 tCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACt                                                                                                                             >  2:3070877/1‑148 (MQ=255)
                 tCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACt                                                                                                                             >  2:2640595/1‑148 (MQ=255)
                 tCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACt                                                                                                                             >  2:2675400/1‑148 (MQ=255)
                 tCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACt                                                                                                                             >  1:391882/1‑148 (MQ=255)
                    cAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtc                                                                                                                          <  1:1621591/148‑1 (MQ=255)
                    cAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtc                                                                                                                          <  1:1479377/148‑1 (MQ=255)
                      ggTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcat                                                                                                                        <  2:2573187/148‑1 (MQ=255)
                          ggAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGa                                                                                                                    <  1:277227/148‑1 (MQ=255)
                             aaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAata                                                                                                                 <  1:24508/148‑1 (MQ=255)
                              aCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCCCGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAatat                                                                                                                >  1:1662851/1‑148 (MQ=255)
                              aCGGAGAAACCAACCAGACCCGCTTTAGGCATTGCGCCGCGCCCGGTGGCACCTTCATAGCGGATTTCACCGCCGGAAAAGGTCGCCGCGCCCGGCCACCGAGACAACTCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAatat                                                                                                                >  1:1662831/1‑148 (MQ=255)
                               cGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                               <  2:1951132/148‑1 (MQ=255)
                                   gaAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGcc                                                                                                           >  2:1824863/1‑148 (MQ=255)
                                      accaaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGt                                                                                                        <  2:1119509/148‑1 (MQ=255)
                                       ccaaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGtt                                                                                                       <  2:897509/148‑1 (MQ=255)
                                       ccaaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGtt                                                                                                       <  1:739174/148‑1 (MQ=255)
                                       ccaaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGtt                                                                                                       <  1:2159033/148‑1 (MQ=255)
                                           ccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTg                                                                                                   <  1:1956718/148‑1 (MQ=255)
                                            caGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGa                                                                                                  >  1:2593563/1‑148 (MQ=255)
                                              gACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGGGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATg                                                                                                >  1:2523142/1‑148 (MQ=255)
                                               aCCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATgg                                                                                               >  1:2663648/1‑148 (MQ=255)
                                               aCCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATgg                                                                                               >  2:1455584/1‑148 (MQ=255)
                                                cccGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGa                                                                                              <  1:1638934/148‑1 (MQ=255)
                                                 ccGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGaa                                                                                             >  1:3122989/1‑148 (MQ=255)
                                                   gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                           >  1:2639082/1‑148 (MQ=255)
                                                   gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                           <  2:878522/148‑1 (MQ=255)
                                                     tttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCg                                                                                         >  2:198153/1‑148 (MQ=255)
                                                      ttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGt                                                                                        <  1:925588/148‑1 (MQ=255)
                                                        cGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAg                                                                                      <  2:24653/148‑1 (MQ=255)
                                                        cGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAg                                                                                      >  1:1696044/1‑148 (MQ=255)
                                                             ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGcc                                                                                 >  2:2711997/1‑148 (MQ=255)
                                                             ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGcc                                                                                 >  2:702463/1‑148 (MQ=255)
                                                               gcgccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCg                                                                               <  1:3325971/148‑1 (MQ=255)
                                                                 gccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGtt                                                                             <  2:223575/148‑1 (MQ=255)
                                                                  ccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGttt                                                                            <  2:1571859/148‑1 (MQ=255)
                                                                   cgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTc                                                                           <  1:584088/148‑1 (MQ=255)
                                                                   cgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTc                                                                           >  1:2023870/1‑148 (MQ=255)
                                                                    gcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCg                                                                          <  2:953326/148‑1 (MQ=255)
                                                                      gccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTg                                                                        <  1:3204825/148‑1 (MQ=255)
                                                                           gTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAg                                                                   >  1:1791053/1‑148 (MQ=255)
                                                                               cACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCAACAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCaaa                                                               >  1:306528/1‑148 (MQ=255)
                                                                                 ccTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGt                                                             <  1:2611559/148‑1 (MQ=255)
                                                                                    tCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGc                                                          >  1:3104642/1‑148 (MQ=255)
                                                                                    tCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGc                                                          >  1:1329613/1‑148 (MQ=255)
                                                                                       tCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGc                                                       <  1:568590/148‑1 (MQ=255)
                                                                                       tCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGc                                                       <  2:1168159/148‑1 (MQ=255)
                                                                                       tCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGc                                                       <  1:568723/148‑1 (MQ=255)
                                                                                          cgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCa                                                    <  1:3465345/148‑1 (MQ=255)
                                                                                            gATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACt                                                  >  1:1031656/1‑148 (MQ=255)
                                                                                            gATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACt                                                  <  2:333330/148‑1 (MQ=255)
                                                                                             aTTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACtt                                                 >  1:2991669/1‑148 (MQ=255)
                                                                                             aTTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACtt                                                 <  2:2009710/148‑1 (MQ=255)
                                                                                              tttCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTc                                                >  1:1430831/1‑148 (MQ=255)
                                                                                                 cACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGa                                             >  2:620429/1‑148 (MQ=255)
                                                                                                 cACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGa                                             >  2:2411095/1‑148 (MQ=255)
                                                                                                  aCCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGaa                                            >  1:750642/1‑148 (MQ=255)
                                                                                                    cgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAcc                                          >  1:3284608/1‑148 (MQ=255)
                                                                                                    cgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAcc                                          >  1:3407405/1‑148 (MQ=255)
                                                                                                      ccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCtt                                        <  2:1694781/148‑1 (MQ=255)
                                                                                                         gAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCa                                     >  2:2055927/1‑148 (MQ=255)
                                                                                                          aaCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAt                                    <  2:2751540/148‑1 (MQ=255)
                                                                                                           aCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAtt                                   <  1:2675400/148‑1 (MQ=255)
                                                                                                               gTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcgg                               <  1:2711997/148‑1 (MQ=255)
                                                                                                                tCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcggc                              >  1:3251665/1‑148 (MQ=255)
                                                                                                                      cgcCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGtt                        <  1:606264/148‑1 (MQ=255)
                                                                                                                         ccGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATc                     <  1:1444745/148‑1 (MQ=255)
                                                                                                                             ccACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTta                 <  1:1590093/148‑1 (MQ=255)
                                                                                                                                cGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATaa              <  2:2154591/148‑1 (MQ=255)
                                                                                                                                 ggAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATTTGCGCCGGTTCCTGATGGTAATCGTAGCGGCCCGTTTCTTGGTCATCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAg             <  2:1662831/148‑1 (MQ=255)
                                                                                                                                 ggAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAg             <  2:1662851/148‑1 (MQ=255)
                                                                                                                                  gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGTCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc            <  2:391882/148‑1 (MQ=255)
                                                                                                                                  gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc            >  1:493516/1‑148 (MQ=255)
                                                                                                                                   agaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCa           <  1:1721473/148‑1 (MQ=255)
                                                                                                                                   agaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCa           <  1:1721280/148‑1 (MQ=255)
                                                                                                                                     aaaTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCaga         >  2:1179872/1‑148 (MQ=255)
                                                                                                                                        tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag      >  2:3269463/1‑148 (MQ=255)
                                                                                                                                         cGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAgagaga     <  1:1983301/148‑1 (MQ=255)
                                                                                                                                            cGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACg  <  1:1019291/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACG  >  CP009273/2687863‑2688150

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: