Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,194,861 G→T 11.9% G333G (GGC→GGA purD ← phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,194,8610GT11.9% 470.1 / 28.2 145G333G (GGC→GGApurDphosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (61/66);  new base T (11/7);  total (72/73)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

CCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATA  >  CP009273/4194714‑4195008
                                                                                                                                                   |                                                                                                                                                   
cccgccGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACg                                                                                                                                                     >  2:2559249/1‑148 (MQ=255)
   gccgTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCg                                                                                                                                                  <  2:1139202/148‑1 (MQ=255)
       ttGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCgagag                                                                                                                                              >  2:424733/1‑148 (MQ=255)
       ttGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCgagag                                                                                                                                              >  2:1433495/1‑148 (MQ=255)
         ggTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGaa                                                                                                                                            <  1:628022/148‑1 (MQ=255)
            cACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAgcg                                                                                                                                         >  1:111620/1‑148 (MQ=255)
            cACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAgcg                                                                                                                                         >  1:2981447/1‑148 (MQ=255)
              cTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGcgcg                                                                                                                                       >  2:1021041/1‑148 (MQ=255)
              cTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGcgcg                                                                                                                                       >  1:2453005/1‑148 (MQ=255)
               tACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGt                                                                                                                                      >  2:3231180/1‑148 (MQ=255)
                 ccTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTc                                                                                                                                    >  1:155807/1‑148 (MQ=255)
                  cTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCa                                                                                                                                   >  1:3215062/1‑148 (MQ=255)
                   tGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCAt                                                                                                                                  >  1:80286/1‑148 (MQ=255)
                    gCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATc                                                                                                                                 >  1:3099799/1‑148 (MQ=255)
                    gCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATc                                                                                                                                 >  1:3099782/1‑148 (MQ=255)
                     cTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATcc                                                                                                                                >  1:2963231/1‑148 (MQ=255)
                       cGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCa                                                                                                                              <  2:579611/148‑1 (MQ=255)
                         tCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACt                                                                                                                            <  2:1295718/148‑1 (MQ=255)
                          cATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTc                                                                                                                           >  1:1923000/1‑148 (MQ=255)
                            tCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCgg                                                                                                                         >  2:2050739/1‑148 (MQ=255)
                             ccGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGa                                                                                                                        <  1:2422364/148‑1 (MQ=255)
                                ccAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGt                                                                                                                     <  2:3027802/148‑1 (MQ=255)
                                 cAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGtt                                                                                                                    <  2:911413/148‑1 (MQ=255)
                                  aGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGttt                                                                                                                   <  1:2103776/148‑1 (MQ=255)
                                   gTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGtttt                                                                                                                  >  1:1955979/1‑148 (MQ=255)
                                   gTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGtttt                                                                                                                  <  2:238709/148‑1 (MQ=255)
                                         tgCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTAGTc                                                                                                            <  2:2233178/148‑1 (MQ=255)
                                             cgcgTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGTAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGt                                                                                                        <  2:2304197/148‑1 (MQ=255)
                                             cgcgTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGt                                                                                                        <  2:34531/148‑1 (MQ=255)
                                              gcgTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGtt                                                                                                       <  2:938454/148‑1 (MQ=255)
                                               cgTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGttt                                                                                                      >  1:2176435/1‑148 (MQ=255)
                                               cgTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGttt                                                                                                      >  1:2176662/1‑148 (MQ=255)
                                                gTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTg                                                                                                     <  2:1631275/148‑1 (MQ=255)
                                                 tGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGc                                                                                                    <  2:2225512/148‑1 (MQ=255)
                                                  ggAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCt                                                                                                   >  1:1229353/1‑148 (MQ=255)
                                                  ggAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCt                                                                                                   >  1:1229347/1‑148 (MQ=255)
                                                   gAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCtt                                                                                                  >  2:1190945/1‑148 (MQ=255)
                                                    aaCACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCttt                                                                                                 <  1:1921714/148‑1 (MQ=255)
                                                    aaCACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCttt                                                                                                 >  1:396792/1‑148 (MQ=255)
                                                    aaCACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCttt                                                                                                 <  1:1921749/148‑1 (MQ=255)
                                                     acacTTTGCCGCCTGCCACGTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAAGCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTc                                                                                                <  2:2207152/148‑1 (MQ=255)
                                                      cacTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTca                                                                                               >  1:2476741/1‑148 (MQ=255)
                                                       acTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTcac                                                                                              <  2:1600814/148‑1 (MQ=255)
                                                        cTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTcaca                                                                                             >  1:56590/1‑148 (MQ=255)
                                                        cTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTcaca                                                                                             <  1:3210197/148‑1 (MQ=255)
                                                          ttGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAgg                                                                                           >  2:3190546/1‑148 (MQ=255)
                                                               gccTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGgccgcc                                                                                      <  2:2450331/148‑1 (MQ=255)
                                                                    ccACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCa                                                                                 >  1:2724891/1‑148 (MQ=255)
                                                                        ttcttcCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGc                                                                             >  2:446493/1‑148 (MQ=255)
                                                                         tcttcCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCt                                                                            >  2:2416329/1‑148 (MQ=255)
                                                                            tcCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCaa                                                                         <  1:1403879/148‑1 (MQ=255)
                                                                            tcCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCaa                                                                         <  1:1405540/148‑1 (MQ=255)
                                                                                gCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAg                                                                     <  1:2068733/148‑1 (MQ=255)
                                                                                   gCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATc                                                                  <  2:1465190/148‑1 (MQ=255)
                                                                                    cAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCg                                                                 <  1:264167/148‑1 (MQ=255)
                                                                                       gCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGAc                                                              <  2:2702477/148‑1 (MQ=255)
                                                                                          gTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTc                                                           <  1:716872/148‑1 (MQ=255)
                                                                                               tCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACCGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATgcg                                                      >  1:2697543/1‑148 (MQ=255)
                                                                                               tCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATgcg                                                      >  2:504222/1‑148 (MQ=255)
                                                                                                     cACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCAt                                                >  2:1885126/1‑148 (MQ=255)
                                                                                                     cACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCAt                                                >  2:216755/1‑148 (MQ=255)
                                                                                                     cACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTGTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCAt                                                >  2:1441865/1‑148 (MQ=255)
                                                                                                      aCCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATa                                               >  2:780108/1‑148 (MQ=255)
                                                                                                       ccGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATaa                                              <  2:1295664/148‑1 (MQ=255)
                                                                                                       ccGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATaa                                              >  2:1531390/1‑148 (MQ=255)
                                                                                                        cGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAAt                                             >  2:612602/1‑148 (MQ=255)
                                                                                                           tGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCgg                                          <  1:2671863/148‑1 (MQ=255)
                                                                                                            gCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGc                                         <  1:3297634/148‑1 (MQ=255)
                                                                                                                tAATCACCCGGATATCCACCCGCAGCCATCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTggg                                     >  2:1138132/1‑148 (MQ=255)
                                                                                                                  aTCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGtt                                   <  2:2689470/148‑1 (MQ=255)
                                                                                                                  aTCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGtt                                   >  2:1285877/1‑148 (MQ=255)
                                                                                                                  aTCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGAACGGACTTCATGCGCAGCATAATCGGCTGGGtt                                   >  2:505829/1‑148 (MQ=255)
                                                                                                                   tCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGttt                                  <  1:2337176/148‑1 (MQ=255)
                                                                                                                   tCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGttt                                  <  2:3099782/148‑1 (MQ=255)
                                                                                                                    cACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTc                                 >  2:797498/1‑148 (MQ=255)
                                                                                                                    cACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACGCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTc                                 <  1:1004739/148‑1 (MQ=255)
                                                                                                                        cGGATATCCACCCGCAGCCATCACCCCGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCCCAGGCCGCCAGGCAGAGCTCACCCAGATCGGACTGCATGCGCAGCATAATCGGCTGGGTTTCCGGa                             <  2:2777184/148‑1 (MQ=255)
                                                                                                                        cGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGTGTTTCCGGa                             <  2:1739384/148‑1 (MQ=255)
                                                                                                                            tatCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGc                         >  2:198062/1‑148 (MQ=255)
                                                                                                                               ccACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCaa                      <  2:1603563/148‑1 (MQ=255)
                                                                                                                                cACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCaaa                     <  1:142328/148‑1 (MQ=255)
                                                                                                                                 aCCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAg                    <  1:1308913/148‑1 (MQ=255)
                                                                                                                                 aCCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAg                    <  1:716757/148‑1 (MQ=255)
                                                                                                                                 aCCCGAAGCCTTCCCCCCGCCGAGAGAAGCGCGTTCATCCCACTCGGACTTTTTCTCGTCCTGTTTGCTTGCACAGGCCGCCTGGCAGAGGTCAACCAGAGCGGAATTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAg                    <  1:717125/148‑1 (MQ=255)
                                                                                                                                        gCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCACCCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGt             <  1:1046339/148‑1 (MQ=255)
                                                                                                                                        gCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGAGCGGACTTCTTGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGt             <  1:738203/148‑1 (MQ=255)
                                                                                                                                         ccATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGtt            >  2:2593873/1‑148 (MQ=255)
                                                                                                                                         ccATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAAGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGtt            >  2:1132972/1‑148 (MQ=255)
                                                                                                                                         caaTCACCACTCCTAGAGAAGCGCGTTCATCCCACTCGGACTTTTTCTCGTCCAGTTTGCTTACACAGGCCGCCCGGCAGCGCTCAACCCGATCGGACTTCATGCGCAGAATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGtt            <  1:646122/146‑1 (MQ=255)
                                                                                                                                           aTCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTaa          <  2:2176662/148‑1 (MQ=255)
                                                                                                                                           aTCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTaa          <  2:396340/148‑1 (MQ=255)
                                                                                                                                           aTCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCACGTCCTGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCACATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTaa          <  1:424733/148‑1 (MQ=255)
                                                                                                                                           aTCACCACGCCGAGAGAAGCGCGTTCATCCCACACCGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCATATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTaa          <  1:1290969/148‑1 (MQ=255)
                                                                                                                                            tCACCACTCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTaaa         >  1:3150247/1‑148 (MQ=255)
                                                                                                                                             caccacGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAAt        >  2:2918858/1‑148 (MQ=255)
                                                                                                                                                cacGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCCGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCa     <  1:2050305/148‑1 (MQ=255)
                                                                                                                                                cacGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCa     <  1:2050739/148‑1 (MQ=255)
                                                                                                                                                  cGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCACGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGAATTCATGCGCACCATAAACGGCTGGGTTTCCGGATCGCCCAACCGGCAGTTAAAATCAAt   >  1:1176942/1‑148 (MQ=255)
                                                                                                                                                   gCCGAGAGACGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGAATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATa  <  1:80760/148‑1 (MQ=255)
                                                                                                                                                   gCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATa  <  1:490407/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATA  >  CP009273/4194714‑4195008

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: