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,0060CA97.1% 584.7 / 9.0 206P299P (CCG→CCT) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (73/126);  total (76/130)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

GCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAC  >  CP009273/2687861‑2688149
                                                                                                                                                 |                                                                                                                                               
 cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAggg                                                                                                                                              <  2:1556697/148‑1 (MQ=255)
 cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAggg                                                                                                                                              <  2:1557074/148‑1 (MQ=255)
 cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAggg                                                                                                                                              >  2:3444500/1‑148 (MQ=255)
   gTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtg                                                                                                                                            >  2:2710002/1‑148 (MQ=255)
   gTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtg                                                                                                                                            <  1:3546510/148‑1 (MQ=255)
    ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                           <  1:1046082/148‑1 (MQ=255)
    ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                           <  1:2374098/148‑1 (MQ=255)
    ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                           >  2:1019255/1‑148 (MQ=255)
    ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                           <  2:3150029/148‑1 (MQ=255)
     tCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtggt                                                                                                                                          <  2:1635277/148‑1 (MQ=255)
      cGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGtt                                                                                                                                         <  1:3012660/148‑1 (MQ=255)
      cGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGtt                                                                                                                                         <  1:2852880/148‑1 (MQ=255)
      cGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGtt                                                                                                                                         <  1:2400120/148‑1 (MQ=255)
      cGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGtt                                                                                                                                         <  1:2716566/148‑1 (MQ=255)
      cGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGtt                                                                                                                                         <  1:3242393/148‑1 (MQ=255)
        aaGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtgt                                                                                                                                       >  1:1753128/1‑148 (MQ=255)
          gCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGa                                                                                                                                     <  1:144704/148‑1 (MQ=255)
           ccAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAg                                                                                                                                    <  2:411823/148‑1 (MQ=255)
            cAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGt                                                                                                                                   >  2:2779854/1‑148 (MQ=255)
             aGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGGAGGGGGGTTGTGAGtt                                                                                                                                  >  1:3517777/1‑148 (MQ=255)
              ggAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGttt                                                                                                                                 >  2:492332/1‑148 (MQ=255)
               gAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTc                                                                                                                                >  2:284629/1‑148 (MQ=255)
               gAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTc                                                                                                                                >  1:3247669/1‑148 (MQ=255)
                 aTTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGa                                                                                                                              >  1:139488/1‑148 (MQ=255)
                  ttCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCCCGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGAc                                                                                                                             >  1:1623961/1‑148 (MQ=255)
                  ttCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGAc                                                                                                                             >  1:52598/1‑148 (MQ=255)
                  ttCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGAc                                                                                                                             <  1:3398724/148‑1 (MQ=255)
                     gCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACttt                                                                                                                          <  1:3489187/148‑1 (MQ=255)
                      cAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtc                                                                                                                         <  1:1883778/148‑1 (MQ=255)
                             gAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGaa                                                                                                                  >  2:1660391/1‑148 (MQ=255)
                             gAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGaa                                                                                                                  <  2:3236340/148‑1 (MQ=255)
                              aaaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAat                                                                                                                 <  2:288206/148‑1 (MQ=255)
                               aaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAata                                                                                                                <  1:1382922/148‑1 (MQ=255)
                                aCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAatat                                                                                                               >  2:3817610/1‑148 (MQ=255)
                                 cGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                              <  1:2132364/148‑1 (MQ=255)
                                 cGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                              <  2:3286201/148‑1 (MQ=255)
                                   gagaAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATgcg                                                                                                            >  2:1229950/1‑148 (MQ=255)
                                      aaaCCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCg                                                                                                         <  1:151558/148‑1 (MQ=255)
                                      aaaCCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGGGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCg                                                                                                         <  1:2602759/148‑1 (MQ=255)
                                       aaccaaccAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCgg                                                                                                        >  1:2962314/1‑148 (MQ=255)
                                          caaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACCGAGAAATCGCCGTAGGGGGGTTGTGAGTTTCGACTTTAATCAGAATATGCGCCGGTTc                                                                                                     >  1:1632647/1‑148 (MQ=255)
                                           aaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTcc                                                                                                    <  2:2644879/148‑1 (MQ=255)
                                           aaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTcc                                                                                                    <  2:885048/148‑1 (MQ=255)
                                             ccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTg                                                                                                  <  1:316851/148‑1 (MQ=255)
                                                 aCCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATgg                                                                                              >  2:2517827/1‑148 (MQ=255)
                                                    cGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGaaa                                                                                           >  1:71422/1‑148 (MQ=255)
                                                     gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                          <  2:2801844/148‑1 (MQ=255)
                                                       tttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCg                                                                                        >  2:2342381/1‑148 (MQ=255)
                                                       tttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCg                                                                                        >  1:2528776/1‑148 (MQ=255)
                                                       tttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGGAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCg                                                                                        >  1:3285049/1‑148 (MQ=255)
                                                        ttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGt                                                                                       >  1:3278363/1‑148 (MQ=255)
                                                        ttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGt                                                                                       <  1:1286899/148‑1 (MQ=255)
                                                         tCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTa                                                                                      <  1:2905603/148‑1 (MQ=255)
                                                         tCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTa                                                                                      <  2:2185172/148‑1 (MQ=255)
                                                               ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGACATCGTAGCGGcc                                                                                >  1:2108466/1‑148 (MQ=255)
                                                                         cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTgg                                                                      <  2:84800/148‑1 (MQ=255)
                                                                         cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTgg                                                                      <  2:1904916/148‑1 (MQ=255)
                                                                            ggTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCa                                                                   <  1:636597/148‑1 (MQ=255)
                                                                              tGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGc                                                                 >  1:1366253/1‑148 (MQ=255)
                                                                                    cTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTa                                                           <  1:3820461/148‑1 (MQ=255)
                                                                                          cgcgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGcc                                                     <  2:2258424/148‑1 (MQ=255)
                                                                                          cgcgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGcc                                                     >  2:34780/1‑148 (MQ=255)
                                                                                           gcgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGccc                                                    >  1:1876178/1‑148 (MQ=255)
                                                                                             ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCAc                                                  <  2:2567662/148‑1 (MQ=255)
                                                                                               aTTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACtt                                                <  1:1581286/148‑1 (MQ=255)
                                                                                                tttCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTc                                               >  1:2047786/1‑148 (MQ=255)
                                                                                                  tCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAg                                             <  1:152193/148‑1 (MQ=255)
                                                                                                     ccgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAc                                          <  2:1649282/148‑1 (MQ=255)
                                                                                                     ccgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAc                                          <  2:1101681/148‑1 (MQ=255)
                                                                                                      cgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAcc                                         <  2:1444932/148‑1 (MQ=255)
                                                                                                      cgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAcc                                         >  1:1138267/1‑148 (MQ=255)
                                                                                                      cgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACGTCCGAAcc                                         >  1:3043402/1‑148 (MQ=255)
                                                                                                         cgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTc                                      <  2:139118/148‑1 (MQ=255)
                                                                                                         cgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTc                                      <  2:2991474/148‑1 (MQ=255)
                                                                                                            aaCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAt                                   <  1:1680256/148‑1 (MQ=255)
                                                                                                            aaCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAt                                   <  1:916774/148‑1 (MQ=255)
                                                                                                              ccGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTCGCGGCCCACTTCAGAACCTTCCATTa                                 >  2:1920947/1‑148 (MQ=255)
                                                                                                              ccGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTa                                 >  1:181470/1‑148 (MQ=255)
                                                                                                              ccGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTa                                 <  2:3424190/148‑1 (MQ=255)
                                                                                                              ccGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTa                                 >  1:909243/1‑148 (MQ=255)
                                                                                                              ccGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCACCAAAGTAGCGGCCCACTTCAGAACCTTCCATTa                                 <  1:1050920/148‑1 (MQ=255)
                                                                                                                 gTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcgg                              <  2:3303629/148‑1 (MQ=255)
                                                                                                                      cgcgCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGgcggcg                         <  2:2791467/148‑1 (MQ=255)
                                                                                                                           ccGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATc                    <  1:2136489/148‑1 (MQ=255)
                                                                                                                            cGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCt                   <  1:3382340/148‑1 (MQ=255)
                                                                                                                               ccACGGAGAAATCGCCGTAGGGTGGTTGTTATTTTCGACTTTCATCAGTATATGCGCCTGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTta                <  2:746870/148‑1 (MQ=255)
                                                                                                                                 aCGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtata              <  1:470876/148‑1 (MQ=255)
                                                                                                                                    gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc           >  2:2244977/1‑148 (MQ=255)
                                                                                                                                    gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc           <  1:2047273/148‑1 (MQ=255)
                                                                                                                                    gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc           <  2:1623961/148‑1 (MQ=255)
                                                                                                                                    gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc           <  1:2800977/148‑1 (MQ=255)
                                                                                                                                      gaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCag         <  1:2710002/148‑1 (MQ=255)
                                                                                                                                       aaaTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCaga        >  1:1586648/1‑148 (MQ=255)
                                                                                                                                       aaaTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCaga        >  2:57918/1‑148 (MQ=255)
                                                                                                                                        aaTCGCCGTCGGGTGGTTGTGAGGTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagag       >  2:727979/1‑148 (MQ=255)
                                                                                                                                          tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag     >  1:2002862/1‑148 (MQ=255)
                                                                                                                                          tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag     <  1:529426/148‑1 (MQ=255)
                                                                                                                                          tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag     >  1:3700852/1‑148 (MQ=255)
                                                                                                                                           cGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAgagaga    <  1:1233892/148‑1 (MQ=255)
                                                                                                                                             ccGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAc  <  1:1019255/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAC  >  CP009273/2687861‑2688149

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