Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,688,006 C→A 91.6% 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,0060CA91.6% 423.1 / 35.3 176P299P (CCG→CCTpurLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (11/4);  new base A (72/85);  total (83/93)
Fisher's exact test for biased strand distribution p-value = 5.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

Only 100 of 177 total aligned reads displayed.

CGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACG  >  CP009273/2687859‑2688150
                                                                                                                                                   |                                                                                                                                                
cGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTc                                                                                                                                                  >  2:1524836/1‑148 (MQ=255)
  gCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAgg                                                                                                                                                >  1:2798772/1‑148 (MQ=255)
   cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCTCTTTCGGCTTTGCGACGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCAGGACACGGAGAAATCGCCGTAggg                                                                                                                                               >  2:1627774/1‑148 (MQ=255)
   cTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAggg                                                                                                                                               <  1:538594/148‑1 (MQ=255)
     gTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtg                                                                                                                                             >  1:1617918/1‑148 (MQ=255)
     gTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtg                                                                                                                                             >  1:2218009/1‑148 (MQ=255)
     gTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtg                                                                                                                                             >  1:1617963/1‑148 (MQ=255)
      ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGggtgg                                                                                                                                            >  2:797025/1‑148 (MQ=255)
      ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                            >  2:1341088/1‑148 (MQ=255)
      ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                            >  2:1644221/1‑148 (MQ=255)
      ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                            <  2:2830427/148‑1 (MQ=255)
      ttCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtgg                                                                                                                                            <  2:2829944/148‑1 (MQ=255)
       tCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGggtggt                                                                                                                                           <  1:918848/148‑1 (MQ=255)
          aaGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTtgt                                                                                                                                        >  1:312810/1‑148 (MQ=255)
            gCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGa                                                                                                                                      >  1:2542236/1‑148 (MQ=255)
             ccAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAg                                                                                                                                     <  1:1710315/148‑1 (MQ=255)
              cAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTTAGt                                                                                                                                    <  2:273724/148‑1 (MQ=255)
              cAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGt                                                                                                                                    <  2:3018527/148‑1 (MQ=255)
              cAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGt                                                                                                                                    <  2:530603/148‑1 (MQ=255)
                ggAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGttt                                                                                                                                  <  2:1779053/148‑1 (MQ=255)
                ggAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGttt                                                                                                                                  <  2:2840987/148‑1 (MQ=255)
                ggAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGttt                                                                                                                                  >  1:374609/1‑148 (MQ=255)
                ggAAGTCGCAGGGTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGttt                                                                                                                                  <  2:2840940/148‑1 (MQ=255)
                   aTTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGa                                                                                                                               >  1:1269031/1‑148 (MQ=255)
                   aTTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGa                                                                                                                               <  2:521594/148‑1 (MQ=255)
                         aGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTtca                                                                                                                         <  2:2035065/148‑1 (MQ=255)
                          ggTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGGTGTGAGTTTCGACTTccat                                                                                                                        >  2:730281/1‑148 (MQ=255)
                                 aaCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAata                                                                                                                 <  2:2432389/148‑1 (MQ=255)
                                  aCGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAatat                                                                                                                >  2:469353/1‑148 (MQ=255)
                                   cGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATg                                                                                                               >  2:2072190/1‑148 (MQ=255)
                                      agaAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATgcgc                                                                                                            >  2:2055588/1‑148 (MQ=255)
                                        aaaCCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCg                                                                                                          >  1:1270892/1‑148 (MQ=255)
                                          accaaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGt                                                                                                        <  2:1759938/148‑1 (MQ=255)
                                          accaaccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGt                                                                                                        <  2:942639/148‑1 (MQ=255)
                                              accaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGACTATGCGCCGGTTCCt                                                                                                    >  2:2576791/1‑148 (MQ=255)
                                               ccaGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTg                                                                                                   <  2:1095268/148‑1 (MQ=255)
                                                caGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGa                                                                                                  >  1:210657/1‑148 (MQ=255)
                                                caGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGa                                                                                                  >  1:210656/1‑148 (MQ=255)
                                                caGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGa                                                                                                  <  1:1422486/148‑1 (MQ=255)
                                                   aCCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATgg                                                                                               >  2:660617/1‑148 (MQ=255)
                                                    cccGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGa                                                                                              >  1:240621/1‑148 (MQ=255)
                                                     ccGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGaa                                                                                             <  2:1125796/148‑1 (MQ=255)
                                                     ccGCTTTCGGCTGTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGaa                                                                                             <  1:2780227/148‑1 (MQ=255)
                                                       gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                           <  1:1500803/148‑1 (MQ=255)
                                                       gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                           <  1:1444268/148‑1 (MQ=255)
                                                       gCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAAt                                                                                           <  1:1443815/148‑1 (MQ=255)
                                                         tttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCg                                                                                         >  1:605853/1‑148 (MQ=255)
                                                          ttCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGt                                                                                        <  1:2370683/148‑1 (MQ=255)
                                                              gCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCg                                                                                    >  1:860023/1‑148 (MQ=255)
                                                              gCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCg                                                                                    >  1:860852/1‑148 (MQ=255)
                                                               cTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCgg                                                                                   >  2:357295/1‑148 (MQ=255)
                                                                 ttGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGcc                                                                                 >  1:1042225/1‑148 (MQ=255)
                                                                        gcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCg                                                                          <  2:2844674/148‑1 (MQ=255)
                                                                        gcgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCg                                                                          <  1:2170861/148‑1 (MQ=255)
                                                                         cgccCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGt                                                                         >  2:2258657/1‑148 (MQ=255)
                                                                           cccGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTgg                                                                       >  1:967111/1‑148 (MQ=255)
                                                                             cGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTc                                                                     >  2:924585/1‑148 (MQ=255)
                                                                             cGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTc                                                                     >  1:2851164/1‑148 (MQ=255)
                                                                                  gCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACGTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCac                                                                >  2:2988169/1‑147 (MQ=255)
                                                                                     ccTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTAAGCAAAGt                                                             >  2:1783423/1‑148 (MQ=255)
                                                                                     ccTTCATCGCGGATTTCACCGCCGGAACCGGGCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGt                                                             >  2:1088813/1‑148 (MQ=255)
                                                                                      cTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTa                                                            <  2:2865042/148‑1 (MQ=255)
                                                                                      cTTAATCGAGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTAGTGGACAGCAAAgaa                                                            >  2:129538/1‑146 (MQ=255)
                                                                                       ttCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAg                                                           <  2:2326987/148‑1 (MQ=255)
                                                                                        tCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGc                                                          <  1:3122750/148‑1 (MQ=255)
                                                                                        tCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGc                                                          <  1:1480571/148‑1 (MQ=255)
                                                                                            cgcgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGcc                                                      <  1:3096331/148‑1 (MQ=255)
                                                                                            cgcgGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGcc                                                      <  2:2926072/148‑1 (MQ=255)
                                                                                                 aTTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACtt                                                 <  1:1185814/148‑1 (MQ=255)
                                                                                                  tttCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTc                                                <  1:2072190/148‑1 (MQ=255)
                                                                                                          ccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCtt                                        <  1:77976/148‑1 (MQ=255)
                                                                                                          ccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCtt                                        <  1:1754241/148‑1 (MQ=255)
                                                                                                          ccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCtt                                        <  2:765465/148‑1 (MQ=255)
                                                                                                          ccgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCtt                                        <  2:3282059/148‑1 (MQ=255)
                                                                                                           cgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGACATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTc                                       >  2:2143657/1‑148 (MQ=255)
                                                                                                           cgGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTc                                       <  2:2724010/148‑1 (MQ=255)
                                                                                                            ggAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTACATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTcc                                      <  1:842587/148‑1 (MQ=255)
                                                                                                             gAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCa                                     >  1:3023411/1‑148 (MQ=255)
                                                                                                              aaCCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCAt                                    <  1:2885996/148‑1 (MQ=255)
                                                                                                                ccGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTa                                  <  2:456673/148‑1 (MQ=255)
                                                                                                                       ccgcGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACggcggc                           >  1:30308/1‑148 (MQ=255)
                                                                                                                          cgcCCGGCCACGGAGAAATCGCCGTAGGGTGTTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGtt                        <  1:2074319/148‑1 (MQ=255)
                                                                                                                           gcCCGGCCACGGAGAAATCGCCGTTGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTa                       <  1:1412703/148‑1 (MQ=255)
                                                                                                                            cccGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTCCGGCGGCGTTAt                      >  2:835661/1‑148 (MQ=255)
                                                                                                                            cccGGCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTAt                      >  2:1496518/1‑148 (MQ=255)
                                                                                                                                gCCACGGAGAAATCGCCGTTGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCttt                  <  1:1935306/148‑1 (MQ=255)
                                                                                                                                gCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGTGGCGTTATCttt                  <  1:1600006/148‑1 (MQ=255)
                                                                                                                                gCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCttt                  <  1:265844/148‑1 (MQ=255)
                                                                                                                                gCCACGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCttt                  <  1:1935670/148‑1 (MQ=255)
                                                                                                                                   aCGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtata               <  2:312810/148‑1 (MQ=255)
                                                                                                                                   aCGGAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTtata               <  2:1980382/148‑1 (MQ=255)
                                                                                                                                     ggAGAAATCGCCGTAGGGTGGTTGTGATTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAg             <  1:2721308/148‑1 (MQ=255)
                                                                                                                                     ggAGAAATCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAg             <  1:1068214/148‑1 (MQ=255)
                                                                                                                                         aaaTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCaga         >  1:3083511/1‑148 (MQ=255)
                                                                                                                                          aaTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagag        <  2:1911619/148‑1 (MQ=255)
                                                                                                                                           aTCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCTGCGGTATCTTTATAAGCataga       >  2:2201989/1‑148 (MQ=255)
                                                                                                                                           aTCGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCagaga       >  2:440955/1‑148 (MQ=255)
                                                                                                                                             cGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAgagaga     <  1:660618/148‑1 (MQ=255)
                                                                                                                                             cGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTCCCATTACGGCGGCGTTATCTTTATAAGCAgagaga     <  1:660617/148‑1 (MQ=255)
                                                                                                                                             cGCCGTAGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAgagaga     <  1:1024316/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGGCTGTTCGAAGCCAGGAATTCGCAGGTTGGAAACGGAGAAACCAACCAGACCCGCTTTCGGCTTTGCGCCGCGCCCGGTGGCACCTTCATCGCGGATTTCACCGCCGGAACCGGTCGCCGCGCCCGGCCACGGAGAAATCGCCGTCGGGTGGTTGTGAGTTTCGACTTTCATCAGAATATGCGCCGGTTCCTGATGGAAATCGTAGCGGCCCGTTTCGTGGTCAGCAAAGTAGCGGCCCACTTCAGAACCTTCCATTACGGCGGCGTTATCTTTATAAGCAGAGAGAACG  >  CP009273/2687859‑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: