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,0060CA95.6% 715.5 / 19.8 251P299P (CCG→CCTpurLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (3/8);  new base A (119/120);  total (122/129)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 251 total aligned reads displayed.

GGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAC  >  CP009273/2687860‑2688149
                                                                                                                                                  |                                                                                                                                               
ggCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAgg                                                                                                                                               >  2:2415999/1‑149 (MQ=255)
  cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggt                                                                                                                                             <  1:3177186/149‑1 (MQ=255)
  cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggt                                                                                                                                             <  1:220080/149‑1 (MQ=255)
     ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtggt                                                                                                                                          <  1:2467029/149‑1 (MQ=255)
      tCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGtt                                                                                                                                         <  2:3714928/149‑1 (MQ=255)
      tCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGtt                                                                                                                                         <  2:1564572/149‑1 (MQ=255)
          aGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGa                                                                                                                                     >  2:3393774/1‑149 (MQ=255)
          aGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGa                                                                                                                                     <  2:3464719/149‑1 (MQ=255)
          aGCCAGGAATTCGCACGTAGGAAACGGAGAAACCAACCATACCCGATTACGGCTTAGCGCCTCGCCCGGTGGCACCTTCAACGCGGATTTCAACGCCGGAACCGGTCGCCCCGCCCCGCCACCGAGAAATCTCCCAAGGGAGGTTGTGa                                                                                                                                     >  2:1344827/1‑149 (MQ=255)
           gCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAg                                                                                                                                    >  1:344454/1‑149 (MQ=255)
             cAGGAAGTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGtt                                                                                                                                  <  2:3186248/149‑1 (MQ=255)
              aGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGttt                                                                                                                                 >  2:2187308/1‑149 (MQ=255)
                gAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCg                                                                                                                               >  1:3595538/1‑149 (MQ=255)
                gAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGGAGGGTGGTTGTGAGTTTCg                                                                                                                               >  1:1073681/1‑149 (MQ=255)
                     cGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACttt                                                                                                                          <  1:1801035/149‑1 (MQ=255)
                      gCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGACCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtc                                                                                                                         >  2:200799/1‑149 (MQ=255)
                        aGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcat                                                                                                                       >  1:1010770/1‑149 (MQ=255)
                         ggTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcatc                                                                                                                      <  1:2333537/149‑1 (MQ=255)
                          gTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcatca                                                                                                                     <  2:3757826/149‑1 (MQ=255)
                          gTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtcatca                                                                                                                     <  1:2772464/149‑1 (MQ=255)
                           ttGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAg                                                                                                                    >  1:4049899/1‑149 (MQ=255)
                           ttGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAg                                                                                                                    >  1:2936597/1‑149 (MQ=255)
                            tGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGa                                                                                                                   >  1:1130309/1‑149 (MQ=255)
                             ggAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGaa                                                                                                                  <  1:3882507/149‑1 (MQ=255)
                               aaaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGTGTTTCGACTTTCATCAGAata                                                                                                                >  1:839582/1‑149 (MQ=255)
                                aaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAACTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAatat                                                                                                               >  2:1500029/1‑149 (MQ=255)
                                 aCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                              <  1:934963/149‑1 (MQ=255)
                                       aaaCCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCgg                                                                                                        <  2:2900769/149‑1 (MQ=255)
                                       aaaCCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCgg                                                                                                        <  2:1584431/149‑1 (MQ=255)
                                        aaccaaccAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGt                                                                                                       >  2:4061287/1‑149 (MQ=255)
                                        aaccaaccAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGt                                                                                                       <  2:2547513/149‑1 (MQ=255)
                                        aaccaaccAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGt                                                                                                       <  2:478327/149‑1 (MQ=255)
                                        aaccaaccAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGACGCGCCCGGCCACGGAGAAATCGCCGGAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGt                                                                                                       >  2:585439/1‑149 (MQ=255)
                                           caaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGGAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTcc                                                                                                    >  1:654670/1‑149 (MQ=255)
                                            ccccaGCCCCGCTTTCGGCATGGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCt                                                                                                   <  2:2399354/147‑1 (MQ=255)
                                             accaGACCCGCTTTCGGCTTTTCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTAACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGGAGGGGGGGTGTGAGTTTCGACTTTAATAAGAATATGCGCCGGTTccag                                                                                                  >  1:1049585/1‑147 (MQ=255)
                                             accaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTg                                                                                                  <  1:2692526/149‑1 (MQ=255)
                                                aGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCCCGGAGAAATCGCCGTAGGGTGGTGGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATg                                                                                               >  2:3780159/1‑149 (MQ=255)
                                                aGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATg                                                                                               >  1:1646584/1‑149 (MQ=255)
                                                aGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAAACGCCGTAGGGTGGTTGTGAGTTTCGAATTTCATCAGAATATGCGCCGGTTCCTGATg                                                                                               >  1:1646599/1‑149 (MQ=255)
                                                 gACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCCCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATgg                                                                                              >  1:193044/1‑149 (MQ=255)
                                                  aCCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGa                                                                                             <  2:3525891/149‑1 (MQ=255)
                                                   cccGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGaa                                                                                            >  1:589224/1‑149 (MQ=255)
                                                    ccGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGaaa                                                                                           <  2:766261/149‑1 (MQ=255)
                                                     cGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                          >  2:2773243/1‑149 (MQ=255)
                                                      gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATc                                                                                         <  2:1659594/149‑1 (MQ=255)
                                                        tttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGt                                                                                       >  2:2968667/1‑149 (MQ=255)
                                                         ttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTa                                                                                      <  2:3215381/149‑1 (MQ=255)
                                                               tttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGcc                                                                                <  2:1425249/149‑1 (MQ=255)
                                                                     ccgcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTc                                                                          <  2:3655210/149‑1 (MQ=255)
                                                                        cgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTg                                                                       <  1:720914/149‑1 (MQ=255)
                                                                          cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                     <  2:1016232/149‑1 (MQ=255)
                                                                          cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                     <  1:685337/149‑1 (MQ=255)
                                                                          cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                     >  1:2175563/1‑149 (MQ=255)
                                                                          cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGt                                                                     <  1:3573188/149‑1 (MQ=255)
                                                                           ccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTc                                                                    >  1:2560820/1‑149 (MQ=255)
                                                                              gTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGc                                                                 >  1:1737151/1‑149 (MQ=255)
                                                                                ggCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCaa                                                               <  2:216799/149‑1 (MQ=255)
                                                                                 gCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCaaa                                                              <  1:1122993/149‑1 (MQ=255)
                                                                                   aCCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGt                                                            >  2:1693218/1‑149 (MQ=255)
                                                                                   aCCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGt                                                            >  1:4096759/1‑149 (MQ=255)
                                                                                     cTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAg                                                          <  2:3265944/149‑1 (MQ=255)
                                                                                     cTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAg                                                          >  1:3513519/1‑149 (MQ=255)
                                                                                       tCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCg                                                        <  2:2050736/149‑1 (MQ=255)
                                                                                         aTCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGc                                                      >  2:2635203/1‑149 (MQ=255)
                                                                                          tCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGcc                                                     >  1:2552852/1‑149 (MQ=255)
                                                                                             cgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCAc                                                  <  1:3860050/149‑1 (MQ=255)
                                                                                              ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACt                                                 >  2:854419/1‑149 (MQ=255)
                                                                                              ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACt                                                 >  2:93712/1‑149 (MQ=255)
                                                                                              ggATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACt                                                 >  1:1057181/1‑149 (MQ=255)
                                                                                                    cACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGaa                                           <  2:176909/149‑1 (MQ=255)
                                                                                                    cACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGaa                                           <  2:3019981/149‑1 (MQ=255)
                                                                                                     aCCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAAc                                          <  2:2709327/149‑1 (MQ=255)
                                                                                                        gccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCtt                                       >  1:2707218/1‑149 (MQ=255)
                                                                                                           ggAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCa                                    <  2:3938339/149‑1 (MQ=255)
                                                                                                            gAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAt                                   <  1:1257855/149‑1 (MQ=255)
                                                                                                            gAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGCGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAt                                   <  2:1105434/149‑1 (MQ=255)
                                                                                                             aaCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAtt                                  <  1:2472404/149‑1 (MQ=255)
                                                                                                                 ggTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTAcgg                              <  1:2611387/149‑1 (MQ=255)
                                                                                                                       cgcgCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGt                        <  2:2343194/149‑1 (MQ=255)
                                                                                                                             cGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCtt                  >  2:2307628/1‑149 (MQ=255)
                                                                                                                             cGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCtt                  <  2:2034795/149‑1 (MQ=255)
                                                                                                                              ggCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGTAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCAACTTCAGAACCTTCCATTACGGAGGCGTTATCttt                 >  2:266170/1‑149 (MQ=255)
                                                                                                                              ggCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCttt                 <  2:1883644/149‑1 (MQ=255)
                                                                                                                                ccACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCCTTACGGCGGCGTTATCTTtat               <  1:3621186/149‑1 (MQ=255)
                                                                                                                                 cACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtata              >  1:3059107/1‑149 (MQ=255)
                                                                                                                                 cACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtata              <  1:2019843/149‑1 (MQ=255)
                                                                                                                                    ggAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCGCTTTTCGTTGTCAGCATAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGc           <  1:2356421/149‑1 (MQ=255)
                                                                                                                                     gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCa          <  1:2714911/149‑1 (MQ=255)
                                                                                                                                     gagaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCa          >  2:3083315/1‑149 (MQ=255)
                                                                                                                                      agaAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCag         <  1:476372/149‑1 (MQ=255)
                                                                                                                                         aaTCGGCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGTGCGGGTTCGTGATGGAAATAGTAGCGGCGCGTTTCGTGGTGAGCATAGTAGCGGCCCACTTCAGAACCTTCGATTACGGCGGCGTTATATTTTTTAGCagaga      <  1:1344827/149‑1 (MQ=255)
                                                                                                                                         aaTCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagaga      >  2:2556392/1‑149 (MQ=255)
                                                                                                                                         aaTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagaga      <  1:3393774/149‑1 (MQ=255)
                                                                                                                                         aaTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagaga      <  1:693768/149‑1 (MQ=255)
                                                                                                                                          aTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAACTCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGCCGGCGTTATATTTATAAGCagagag     >  2:3185467/1‑149 (MQ=255)
                                                                                                                                          aTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagagag     >  2:3185190/1‑149 (MQ=255)
                                                                                                                                           tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAgagaga    <  2:3434714/149‑1 (MQ=255)
                                                                                                                                           tCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAgagaga    >  1:2649473/1‑149 (MQ=255)
                                                                                                                                             gCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAc  >  2:404510/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAAC  >  CP009273/2687860‑2688149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: