Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,394,875 G→T 100% A124A (GCG→GCT)  purA → adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,394,8750GT99.6% 817.2 / ‑4.4 245A124A (GCG→GCT) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (116/128);  total (117/128)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-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 246 total aligned reads displayed.

CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAG  >  CP009273/4394730‑4395019
                                                                                                                                                 |                                                                                                                                                
  ccGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAg                                                                                                                                              <  2:1950199/148‑1 (MQ=255)
  ccGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAg                                                                                                                                              >  1:1180465/1‑148 (MQ=255)
  ccGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAg                                                                                                                                              <  1:2723006/148‑1 (MQ=255)
  ccGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAg                                                                                                                                              <  2:851502/148‑1 (MQ=255)
     gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGa                                                                                                                                           <  1:2063895/148‑1 (MQ=255)
     gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGa                                                                                                                                           <  2:173417/148‑1 (MQ=255)
     gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGa                                                                                                                                           <  2:968535/148‑1 (MQ=255)
         cgcTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCgg                                                                                                                                       <  2:2508380/148‑1 (MQ=255)
             gatgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacc                                                                                                                                   >  2:3345038/1‑148 (MQ=255)
              atgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                  >  2:1754854/1‑148 (MQ=255)
              atgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                  <  2:1539281/148‑1 (MQ=255)
              atgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                  >  2:1147864/1‑148 (MQ=255)
              atgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                  >  1:1947075/1‑148 (MQ=255)
               tgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcaccac                                                                                                                                 >  2:1165566/1‑148 (MQ=255)
                gaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCaccacc                                                                                                                                >  1:79991/1‑148 (MQ=255)
                  aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCgg                                                                                                                              <  2:431075/148‑1 (MQ=255)
                   agagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGt                                                                                                                             <  1:1126953/148‑1 (MQ=255)
                   agagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGt                                                                                                                             >  2:300672/1‑148 (MQ=255)
                      gatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGt                                                                                                                          <  1:3495338/148‑1 (MQ=255)
                         gaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGt                                                                                                                       >  1:696938/1‑148 (MQ=255)
                          aaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTa                                                                                                                      >  2:1478461/1‑148 (MQ=255)
                           aagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTAt                                                                                                                     <  2:1221839/148‑1 (MQ=255)
                            agaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATc                                                                                                                    >  2:2983526/1‑148 (MQ=255)
                                cTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGc                                                                                                                <  1:898812/148‑1 (MQ=255)
                                 tGGTAGAGCGTGGCATCCGCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGGTGATACTTGATTATAACGTTGCGCTGGATAACGCGCTTGAGAAAGCGCGTTGCGCTAAAGCGATCGGCACCAACGTTCGTGTTATCTGGcc                                                                                                               <  2:2790605/148‑1 (MQ=255)
                                  ggAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCt                                                                                                              >  2:2780366/1‑148 (MQ=255)
                                  ggAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCt                                                                                                              >  1:3157899/1‑148 (MQ=255)
                                   gAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTg                                                                                                             <  2:515919/148‑1 (MQ=255)
                                   gAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTg                                                                                                             >  1:3318813/1‑148 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                         >  1:2360261/1‑148 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                         >  2:2260879/1‑148 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                         <  1:2166788/148‑1 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                         >  1:359742/1‑148 (MQ=255)
                                        ccGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTAt                                                                                                        <  1:442711/148‑1 (MQ=255)
                                         cGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATg                                                                                                       <  2:295426/148‑1 (MQ=255)
                                          gTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGa                                                                                                      <  1:2229958/148‑1 (MQ=255)
                                           tGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGaa                                                                                                     >  1:2526650/1‑148 (MQ=255)
                                            ggCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAg                                                                                                    >  1:2776856/1‑148 (MQ=255)
                                             gCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGa                                                                                                   >  2:1317949/1‑148 (MQ=255)
                                              cATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGAt                                                                                                  <  2:1574989/148‑1 (MQ=255)
                                               aTCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATa                                                                                                 <  2:2411909/148‑1 (MQ=255)
                                                tCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaa                                                                                                >  2:3169511/1‑148 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                               >  2:2044014/1‑148 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                               >  2:1290004/1‑148 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                               >  2:3483339/1‑148 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                               >  1:1462510/1‑148 (MQ=255)
                                                   ccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                             <  1:708792/148‑1 (MQ=255)
                                                   ccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                             >  2:1212026/1‑148 (MQ=255)
                                                   ccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                             <  2:1676498/148‑1 (MQ=255)
                                                    cGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTa                                                                                            <  1:2861813/148‑1 (MQ=255)
                                                       tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCa                                                                                         >  2:2650182/1‑148 (MQ=255)
                                                          tGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGt                                                                                      >  2:3024911/1‑148 (MQ=255)
                                                           gAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCATATGAAGATAAAGTAGCACGTc                                                                                     >  1:2944774/1‑148 (MQ=255)
                                                            aGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGAGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcg                                                                                    >  1:516475/1‑148 (MQ=255)
                                                            aGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGAGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcg                                                                                    >  1:2281645/1‑148 (MQ=255)
                                                                 ctgctgctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTc                                                                               >  1:526604/1‑148 (MQ=255)
                                                                      gctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCgt                                                                          <  1:979949/148‑1 (MQ=255)
                                                                         gTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGtt                                                                       >  2:2336372/1‑148 (MQ=255)
                                                                         gTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGtt                                                                       >  2:571845/1‑148 (MQ=255)
                                                                               aGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAc                                                                 <  1:3398528/148‑1 (MQ=255)
                                                                                   tGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCttt                                                             <  1:3239427/148‑1 (MQ=255)
                                                                                    gTCCGCTGATCCTTGATTATCACGTTGCGATGGATAACGCGCGTGATAAAGCGCGTGGCGCTAAAGAGATCGGCACCACCGGTCGTGGTATCGTTCCTGCTTATGAAGATAAAGTAGCACGTCGCGTTCTGCGTGTTGGCGACCtttt                                                            <  1:2709396/148‑1 (MQ=255)
                                                                                     tCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTc                                                           >  1:530278/1‑148 (MQ=255)
                                                                                       cGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGa                                                         <  1:1856298/148‑1 (MQ=255)
                                                                                          tGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACaa                                                      <  1:1165566/148‑1 (MQ=255)
                                                                                            aTCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAg                                                    <  1:1499835/148‑1 (MQ=255)
                                                                                             tCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGa                                                   <  2:2612442/148‑1 (MQ=255)
                                                                                              ccTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaa                                                  >  2:2512649/1‑148 (MQ=255)
                                                                                               cTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaaa                                                 >  2:2099288/1‑148 (MQ=255)
                                                                                               cTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaaa                                                 >  2:2170439/1‑148 (MQ=255)
                                                                                               cTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaaa                                                 <  1:1308686/148‑1 (MQ=255)
                                                                                                 tGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAcc                                               <  2:555143/148‑1 (MQ=255)
                                                                                                 tGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAcc                                               <  1:1969740/148‑1 (MQ=255)
                                                                                                 tGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAcc                                               >  2:756988/1‑148 (MQ=255)
                                                                                                  gATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCt                                              >  1:4941/1‑148 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                             <  1:2614011/148‑1 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                             >  2:2458458/1‑148 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                             >  2:64841/1‑148 (MQ=255)
                                                                                                     tATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCg                                           >  2:45934/1‑148 (MQ=255)
                                                                                                     tATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCg                                           >  1:3001822/1‑148 (MQ=255)
                                                                                                     tATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCg                                           <  1:300575/148‑1 (MQ=255)
                                                                                                     tATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCg                                           <  1:3091830/148‑1 (MQ=255)
                                                                                                          cGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGat                                      <  1:1701657/148‑2 (MQ=255)
                                                                                                             tGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaaa                                   <  1:463664/148‑1 (MQ=255)
                                                                                                                    gATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaag                            <  2:3004049/148‑1 (MQ=255)
                                                                                                                    gATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaag                            >  2:1907620/1‑148 (MQ=255)
                                                                                                                       aaCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaag                         >  2:874421/1‑148 (MQ=255)
                                                                                                                       aaCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaag                         <  2:1425096/148‑1 (MQ=255)
                                                                                                                        aCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                        <  1:1549915/148‑1 (MQ=255)
                                                                                                                        aCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                        <  1:2799866/148‑1 (MQ=255)
                                                                                                                         cgcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTg                       <  1:1824197/148‑1 (MQ=255)
                                                                                                                         cgcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTg                       <  1:47273/148‑1 (MQ=255)
                                                                                                                         cgcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGAGTTGGCGACCTTTTCGACCAAGAAACCTTCGCTGAAAAACTGAAAGAAGTg                       >  2:2241341/1‑148 (MQ=255)
                                                                                                                             cgTGAGAAAGCGAGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATgg                   <  1:1291123/148‑1 (MQ=255)
                                                                                                                                      gcgcGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca          >  1:1627495/1‑148 (MQ=255)
                                                                                                                                      gcgcGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca          <  1:1528282/148‑1 (MQ=255)
                                                                                                                                      gcgcGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca          <  1:418741/148‑1 (MQ=255)
                                                                                                                                       cgcgTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACaa         <  1:1582638/148‑1 (MQ=255)
                                                                                                                                             gcgcTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCa   >  2:2128671/1‑148 (MQ=255)
                                                                                                                                             gcgcTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCa   <  2:993409/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAG  >  CP009273/4394730‑4395019

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: ‑