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.5% 562.5 / 11.6 200P299P (CCG→CCT) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base A (89/102);  total (91/109)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 200 total aligned reads displayed.

GCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACG  >  CP009273/2687861‑2688150
                                                                                                                                                 |                                                                                                                                                
 cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAACTCGCCGTAggg                                                                                                                                               >  2:1997146/1‑148 (MQ=255)
  tGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggt                                                                                                                                              >  1:1575226/1‑148 (MQ=255)
   gTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCGTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtg                                                                                                                                             >  1:1835265/1‑148 (MQ=255)
    ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                            >  1:2091969/1‑148 (MQ=255)
       gAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtg                                                                                                                                         <  2:2191649/148‑1 (MQ=255)
       gAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtg                                                                                                                                         <  2:3551543/148‑1 (MQ=255)
        aaGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtgt                                                                                                                                        >  2:3146772/1‑148 (MQ=255)
         aGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtgtg                                                                                                                                       >  2:340180/1‑148 (MQ=255)
          gCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGa                                                                                                                                      >  2:678738/1‑148 (MQ=255)
             aGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGtt                                                                                                                                   <  2:1794987/148‑1 (MQ=255)
               gAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTc                                                                                                                                 >  1:3238507/1‑148 (MQ=255)
                   tCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACt                                                                                                                             <  2:1642191/148‑1 (MQ=255)
                    cGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACtt                                                                                                                            <  1:1650555/148‑1 (MQ=255)
                     gCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACttt                                                                                                                           <  2:1721307/148‑1 (MQ=255)
                      cAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtc                                                                                                                          <  2:508344/148‑1 (MQ=255)
                        ggTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcat                                                                                                                        <  2:3059875/148‑1 (MQ=255)
                        ggTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcat                                                                                                                        <  2:1289340/148‑1 (MQ=255)
                         gTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcatc                                                                                                                       >  2:2065573/1‑148 (MQ=255)
                         gTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcatc                                                                                                                       >  2:270581/1‑148 (MQ=255)
                          ttGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcatca                                                                                                                      >  2:1493980/1‑148 (MQ=255)
                              aaaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAat                                                                                                                  <  1:61893/148‑1 (MQ=255)
                               aaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAata                                                                                                                 <  1:2326056/148‑1 (MQ=255)
                               aaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAata                                                                                                                 <  1:2876996/148‑1 (MQ=255)
                                 cGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                               >  2:1585240/1‑148 (MQ=255)
                                 cGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                               >  1:1137542/1‑148 (MQ=255)
                                 cGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                               >  2:2549652/1‑148 (MQ=255)
                                   gagaAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATgcg                                                                                                             >  2:1601766/1‑148 (MQ=255)
                                         ccaaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGtt                                                                                                       <  1:1285695/148‑1 (MQ=255)
                                          caaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTc                                                                                                      <  2:3397677/148‑1 (MQ=255)
                                          caaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGGAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTc                                                                                                      >  1:3009653/1‑148 (MQ=255)
                                          caaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGGAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTc                                                                                                      >  1:416945/1‑148 (MQ=255)
                                            accaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCt                                                                                                    <  1:1894651/148‑1 (MQ=255)
                                              caGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGa                                                                                                  <  2:1202824/148‑1 (MQ=255)
                                               aGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGAt                                                                                                 <  1:3412708/148‑1 (MQ=255)
                                                gACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATg                                                                                                >  2:80074/1‑148 (MQ=255)
                                                 aCCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGGTCCTGATgg                                                                                               >  2:2165016/1‑148 (MQ=255)
                                                   ccGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGaa                                                                                             <  2:1028266/148‑1 (MQ=255)
                                                     gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                           >  2:2778542/1‑148 (MQ=255)
                                                        ttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGt                                                                                        >  2:3026647/1‑148 (MQ=255)
                                                         tCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTa                                                                                       <  2:1435145/148‑1 (MQ=255)
                                                            gCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCg                                                                                    >  2:403482/1‑148 (MQ=255)
                                                             cTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCgg                                                                                   >  1:1583713/1‑148 (MQ=255)
                                                             cTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCgg                                                                                   >  2:577279/1‑148 (MQ=255)
                                                               ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGcc                                                                                 >  1:2764750/1‑148 (MQ=255)
                                                               ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGcc                                                                                 >  2:627300/1‑148 (MQ=255)
                                                                  cgccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGt                                                                              <  2:2804734/148‑1 (MQ=255)
                                                                    ccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGttt                                                                            <  1:2690234/148‑1 (MQ=255)
                                                                    ccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGttt                                                                            >  2:2529875/1‑148 (MQ=255)
                                                                    ccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGAATTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGttt                                                                            <  1:24004/148‑1 (MQ=255)
                                                                        gccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTg                                                                        <  1:2908636/148‑1 (MQ=255)
                                                                        gccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTg                                                                        <  1:71236/148‑1 (MQ=255)
                                                                         cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTgg                                                                       <  1:1078285/148‑1 (MQ=255)
                                                                          ccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                      >  1:2033519/1‑148 (MQ=255)
                                                                          ccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                      <  2:3239556/148‑1 (MQ=255)
                                                                          ccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                      <  2:189474/148‑1 (MQ=255)
                                                                          ccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                      <  1:430931/148‑1 (MQ=255)
                                                                              tGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGc                                                                  >  2:2393421/1‑148 (MQ=255)
                                                                                gCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCaa                                                                >  1:2505303/1‑148 (MQ=255)
                                                                                gCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCaa                                                                >  1:3002255/1‑148 (MQ=255)
                                                                                  aCCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAg                                                              <  1:1765207/148‑1 (MQ=255)
                                                                                    cTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTa                                                            >  2:2954448/1‑148 (MQ=255)
                                                                                        aTCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCgg                                                        <  2:2412765/148‑1 (MQ=255)
                                                                                         tCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGTTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGc                                                       <  1:198952/148‑1 (MQ=255)
                                                                                          cgcgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGcc                                                      <  1:1539809/148‑1 (MQ=255)
                                                                                          cgcgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGcc                                                      >  1:3103496/1‑148 (MQ=255)
                                                                                            cgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCa                                                    <  2:3381780/148‑1 (MQ=255)
                                                                                             ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCAc                                                   <  2:3002255/148‑1 (MQ=255)
                                                                                              gATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACt                                                  >  1:2772506/1‑148 (MQ=255)
                                                                                                   cACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGa                                             <  2:920889/148‑1 (MQ=255)
                                                                                                      cgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAcc                                          >  1:615591/1‑148 (MQ=255)
                                                                                                      cgccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCAGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAcc                                          <  2:2642812/148‑1 (MQ=255)
                                                                                                            aaCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAt                                    <  2:1302251/148‑1 (MQ=255)
                                                                                                             aCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAtt                                   >  1:2810711/1‑148 (MQ=255)
                                                                                                              ccGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTTGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTa                                  <  1:927274/148‑1 (MQ=255)
                                                                                                                ggTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcg                                <  2:2864288/148‑1 (MQ=255)
                                                                                                                 gTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcgg                               <  2:2494319/148‑1 (MQ=255)
                                                                                                                   cgccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcggcg                             >  2:1295670/1‑148 (MQ=255)
                                                                                                                    gccgcGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcggcgg                            <  1:2570423/148‑1 (MQ=255)
                                                                                                                     ccgcGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACggcggc                           <  1:2716933/148‑1 (MQ=255)
                                                                                                                      cgcgCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGgcggcg                          <  2:1356960/148‑1 (MQ=255)
                                                                                                                               ccACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTta                 >  1:514982/1‑148 (MQ=255)
                                                                                                                                cACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGTCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtat                <  2:2486853/148‑1 (MQ=255)
                                                                                                                                cACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGTCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtat                <  2:2350605/148‑1 (MQ=255)
                                                                                                                                 aCGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtata               <  2:83997/148‑1 (MQ=255)
                                                                                                                                  cGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATaa              >  2:1164514/1‑148 (MQ=255)
                                                                                                                                   ggAGAAATCGCCGTTGGGTGTTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTTATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAg             <  1:10172/148‑1 (MQ=255)
                                                                                                                                    gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc            >  1:3377981/1‑148 (MQ=255)
                                                                                                                                     agaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGTGGCGTTATCTTTATAAGCa           >  2:2371175/1‑148 (MQ=255)
                                                                                                                                     agaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCa           >  2:1379668/1‑148 (MQ=255)
                                                                                                                                     agaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCa           <  2:2598008/148‑1 (MQ=255)
                                                                                                                                      gaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCag          >  2:3234542/1‑148 (MQ=255)
                                                                                                                                          tCGCCTTAGGGTTGTTGTGAGTTTCGACTTTCATGAGAATATGCGCGGGTTCCTGATGGAAATCGTAGCGGCACGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATATGCagagag      <  2:1781942/148‑1 (MQ=255)
                                                                                                                                          tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAATCagagag      >  1:3297465/1‑148 (MQ=255)
                                                                                                                                          tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag      >  2:1524382/1‑148 (MQ=255)
                                                                                                                                          tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag      >  1:3297126/1‑148 (MQ=255)
                                                                                                                                           cGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAgagaga     <  2:941131/148‑1 (MQ=255)
                                                                                                                                            gCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGaa    <  2:2716264/148‑1 (MQ=255)
                                                                                                                                            gCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGaa    <  2:1835265/148‑1 (MQ=255)
                                                                                                                                              cGTCGGGTGGTTGTGTGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACg  <  1:2620250/148‑1 (MQ=255)
                                                                                                                                              cGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACg  <  1:3146772/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACG  >  CP009273/2687861‑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: ‑