Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,422,841 C→A 100% D293Y (GAT→TAT)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,422,8410CA99.6% 893.6 / ‑3.9 249D293Y (GAT→TAT) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (111/137);  total (112/137)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 249 total aligned reads displayed.

GCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCAC  >  CP009273/2422694‑2422987
                                                                                                                                                   |                                                                                                                                                  
gctgctGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATAtt                                                                                                                                                  <  1:1189021/150‑1 (MQ=255)
  tgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTcc                                                                                                                                                >  2:2780343/1‑150 (MQ=255)
    ctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCa                                                                                                                                              >  2:1838588/1‑150 (MQ=255)
        cccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTcg                                                                                                                                          >  2:4015143/1‑150 (MQ=255)
         ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTcgc                                                                                                                                         <  2:1013193/150‑1 (MQ=255)
              aTGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCa                                                                                                                                    <  2:1530976/150‑1 (MQ=255)
              aTGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCa                                                                                                                                    <  2:2477056/150‑1 (MQ=255)
              aTGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCa                                                                                                                                    <  2:3475111/150‑1 (MQ=255)
               tGATAAAGGTGCGGCCAACATAGCGTTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGAGGGCaa                                                                                                                                   >  1:1766986/1‑150 (MQ=255)
                gATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAt                                                                                                                                  <  1:1406795/150‑1 (MQ=255)
                gATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAt                                                                                                                                  <  2:962309/150‑1 (MQ=255)
                 aTAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAtt                                                                                                                                 <  2:3399027/150‑1 (MQ=255)
                   aaaGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAAtttt                                                                                                                               <  2:1098811/150‑1 (MQ=255)
                    aaGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtc                                                                                                                              <  1:2674929/150‑1 (MQ=255)
                     aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtct                                                                                                                             <  2:796347/150‑1 (MQ=255)
                          cGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGccc                                                                                                                        <  2:1979878/150‑1 (MQ=255)
                           ggCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCa                                                                                                                       <  2:3319205/150‑1 (MQ=255)
                           ggCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCa                                                                                                                       <  2:2116915/150‑1 (MQ=255)
                           ggCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCa                                                                                                                       <  2:1672279/150‑1 (MQ=255)
                              cAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGCTCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGtt                                                                                                                    >  2:2411598/1‑150 (MQ=255)
                                aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTc                                                                                                                  >  1:3241431/1‑150 (MQ=255)
                                aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTc                                                                                                                  <  2:1355600/150‑1 (MQ=255)
                                 cATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCg                                                                                                                 >  1:254093/1‑150 (MQ=255)
                                 cATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCg                                                                                                                 >  2:1209256/1‑150 (MQ=255)
                                           tttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATAtt                                                                                                       >  1:3403417/1‑150 (MQ=255)
                                            ttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTc                                                                                                      >  1:2289271/1‑150 (MQ=255)
                                               aCGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTcaca                                                                                                   >  1:1707120/1‑150 (MQ=255)
                                                cGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCacac                                                                                                  >  2:1764198/1‑150 (MQ=255)
                                                 gAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcg                                                                                                 >  2:2622713/1‑150 (MQ=255)
                                                   aGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgcg                                                                                               >  1:1353046/1‑150 (MQ=255)
                                                    gCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgcgc                                                                                              <  1:1024090/150‑1 (MQ=255)
                                                      ccTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACgcgcgc                                                                                            <  2:3802458/150‑1 (MQ=255)
                                                            ggTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTaaa                                                                                      >  1:698646/1‑150 (MQ=255)
                                                             gTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAAc                                                                                     >  1:165301/1‑150 (MQ=255)
                                                                 ggTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaa                                                                                 >  1:3839958/1‑150 (MQ=255)
                                                                 ggTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaa                                                                                 >  1:1501219/1‑150 (MQ=255)
                                                                  gTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaaa                                                                                <  2:1994195/150‑1 (MQ=255)
                                                                  gTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaaa                                                                                >  1:1637924/1‑150 (MQ=255)
                                                                  gTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaaa                                                                                >  1:517012/1‑150 (MQ=255)
                                                                      gCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAAtttt                                                                            <  2:1282694/150‑1 (MQ=255)
                                                                      gCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAAttgt                                                                            >  2:2080857/1‑148 (MQ=255)
                                                                       cccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTg                                                                           <  1:3715987/150‑1 (MQ=255)
                                                                        ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                          >  1:2265125/1‑150 (MQ=255)
                                                                        ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACAGAAATTTTGt                                                                          <  2:473523/150‑1 (MQ=255)
                                                                        ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACAGAAATTTTGt                                                                          <  2:2886751/150‑1 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                         <  1:353676/150‑1 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                         <  2:1313931/150‑1 (MQ=255)
                                                                         cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTc                                                                         >  2:676707/1‑150 (MQ=255)
                                                                           gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAANTTTGTCGa                                                                       <  1:3075543/150‑1 (MQ=255)
                                                                             aTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATa                                                                     <  1:1328081/150‑1 (MQ=255)
                                                                              tACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCCCACGCGCGCTGTAAACGGAAATTTTGTCGATaa                                                                    <  2:3043405/150‑1 (MQ=255)
                                                                              tACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATaa                                                                    >  2:3549173/1‑150 (MQ=255)
                                                                               aCGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATaaa                                                                   <  2:204894/150‑1 (MQ=255)
                                                                                cGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAAc                                                                  <  2:3716146/150‑1 (MQ=255)
                                                                                   gCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAg                                                               >  2:2520992/1‑150 (MQ=255)
                                                                                   gCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAg                                                               >  2:1055173/1‑150 (MQ=255)
                                                                                   gCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAg                                                               >  1:4112958/1‑150 (MQ=255)
                                                                                    cAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGt                                                              >  1:2888783/1‑150 (MQ=255)
                                                                                     aaTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCTGGCAATTTTCTCGCCCAGTTTCGTGCCCATCTTCACACTCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTc                                                             <  1:4100851/150‑1 (MQ=255)
                                                                                      aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTcc                                                            <  2:934553/150‑1 (MQ=255)
                                                                                      aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTcc                                                            <  2:934707/150‑1 (MQ=255)
                                                                                          ccAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggc                                                        >  1:423927/1‑150 (MQ=255)
                                                                                           cAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcg                                                       <  2:1412462/150‑1 (MQ=255)
                                                                                            aGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcgg                                                      >  2:401354/1‑150 (MQ=255)
                                                                                            aGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcgg                                                      >  2:3742546/1‑150 (MQ=255)
                                                                                             gcgcGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcggg                                                     <  2:2544194/150‑1 (MQ=255)
                                                                                              cgcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGTTATTCCCATTCGCGGGCAATTTTCTCGCCGAGTTTCGGGCCCATCGTCAGAGGCGCGCTGTACTCGGAAAGTTTGTCGAGGAACGTGTCCgggcgggc                                                    >  2:1889219/1‑150 (MQ=255)
                                                                                               gcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCa                                                   <  1:1079365/150‑1 (MQ=255)
                                                                                                  atatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAg                                                >  2:3066132/1‑150 (MQ=255)
                                                                                                  atatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAg                                                >  1:2741791/1‑150 (MQ=255)
                                                                                                    atCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGta                                              <  1:3886344/150‑1 (MQ=255)
                                                                                                     tCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtat                                             <  1:482392/150‑1 (MQ=255)
                                                                                                     tCACACGAGGTGTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtat                                             <  2:2513523/150‑1 (MQ=255)
                                                                                                      cacaCGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtata                                            >  1:853949/1‑150 (MQ=255)
                                                                                                         acGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacat                                         <  2:1593959/150‑1 (MQ=255)
                                                                                                         acGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacat                                         <  1:2621061/150‑1 (MQ=255)
                                                                                                          cGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCAGGGCGGGCAAAGTatacata                                        >  1:4169216/1‑150 (MQ=255)
                                                                                                             ggTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTc                                     >  2:3628585/1‑150 (MQ=255)
                                                                                                             ggTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTc                                     <  2:1568273/150‑1 (MQ=255)
                                                                                                             ggTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTc                                     <  2:3258426/150‑1 (MQ=255)
                                                                                                             ggTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTc                                     >  1:2590877/1‑150 (MQ=255)
                                                                                                               tttCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCaa                                   >  1:1594272/1‑150 (MQ=255)
                                                                                                                ttCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCaaa                                  >  2:671129/1‑150 (MQ=255)
                                                                                                                ttCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCaaa                                  >  2:493695/1‑150 (MQ=255)
                                                                                                                       atcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCAc                           <  2:2345013/150‑1 (MQ=255)
                                                                                                                       atcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCAc                           <  2:761298/150‑1 (MQ=255)
                                                                                                                         cgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCAAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACgg                         >  1:1916526/1‑150 (MQ=255)
                                                                                                                                   acGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCACATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGAcc               <  2:1501742/150‑1 (MQ=255)
                                                                                                                                     gTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCgg             <  2:165281/150‑1 (MQ=255)
                                                                                                                                       cGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAt           >  1:3633266/1‑150 (MQ=255)
                                                                                                                                        gATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGAtt          <  1:3325528/150‑1 (MQ=255)
                                                                                                                                           atCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTc       >  2:973219/1‑150 (MQ=255)
                                                                                                                                           atCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTc       <  2:2856391/150‑1 (MQ=255)
                                                                                                                                            tCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCa      >  2:3856792/1‑150 (MQ=255)
                                                                                                                                            tCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCAATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCa      >  1:2967824/1‑150 (MQ=255)
                                                                                                                                             ccAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg     >  2:3766866/1‑150 (MQ=255)
                                                                                                                                             ccAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg     <  2:190064/150‑1 (MQ=255)
                                                                                                                                             ccAGATATTCCCAGTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAg     <  2:2071827/150‑1 (MQ=255)
                                                                                                                                              cAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGc    <  2:3267893/150‑1 (MQ=255)
                                                                                                                                                gATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACACAAACAGGCACGGATTGCTGACCGGATTGTCAGcac  >  1:714691/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCAC  >  CP009273/2422694‑2422987

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: