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,8750GT98.6% 961.2 / 2.8 285A124A (GCG→GCT) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base T (113/168);  total (114/171)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
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 287 total aligned reads displayed.

CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAGT  >  CP009273/4394730‑4395020
                                                                                                                                                 |                                                                                                                                                 
   cGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCg                                                                                                                                             <  1:114480/149‑1 (MQ=255)
    ggCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGa                                                                                                                                            <  1:4506023/149‑1 (MQ=255)
     gCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGAt                                                                                                                                           <  1:2821595/149‑1 (MQ=255)
       cgcgCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCg                                                                                                                                         <  2:3917866/149‑1 (MQ=255)
                 aaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCgg                                                                                                                               <  2:3202308/149‑1 (MQ=255)
                  aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGt                                                                                                                              <  2:3676003/149‑1 (MQ=255)
                  aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGt                                                                                                                              <  2:269321/149‑1 (MQ=255)
                  aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGt                                                                                                                              >  1:3492409/1‑149 (MQ=255)
                   agagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTc                                                                                                                             >  1:593319/1‑149 (MQ=255)
                    gagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCg                                                                                                                            >  2:605766/1‑149 (MQ=255)
                       atgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTgg                                                                                                                         >  1:1029765/1‑149 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGt                                                                                                                        >  2:4383316/1‑149 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGt                                                                                                                        >  1:2652138/1‑149 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGt                                                                                                                        >  1:259058/1‑149 (MQ=255)
                        tgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGt                                                                                                                        <  1:3817838/149‑1 (MQ=255)
                           aagaACTGGAAGACCGTGGAATCCACTTTAGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATc                                                                                                                     <  1:4542166/149‑1 (MQ=255)
                            agaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCg                                                                                                                    <  2:2042839/149‑1 (MQ=255)
                             gaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTCTCgg                                                                                                                   >  1:1047440/1‑149 (MQ=255)
                              aaCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCggg                                                                                                                  >  2:3485496/1‑149 (MQ=255)
                               aCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGc                                                                                                                 >  2:202337/1‑149 (MQ=255)
                               aCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGc                                                                                                                 >  2:2650248/1‑149 (MQ=255)
                               aCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGc                                                                                                                 >  1:2891790/1‑149 (MQ=255)
                                 tGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCt                                                                                                               >  1:628925/1‑149 (MQ=255)
                                    aaGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCt                                                                                                            <  2:3365837/149‑1 (MQ=255)
                                      gACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                          >  1:441639/1‑149 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTAt                                                                                                         >  2:725346/1‑149 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTAt                                                                                                         <  2:2318224/149‑1 (MQ=255)
                                        ccGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATg                                                                                                        <  2:2622386/149‑1 (MQ=255)
                                            ggCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGa                                                                                                    >  1:2087210/1‑149 (MQ=255)
                                              cATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATa                                                                                                  >  1:3624000/1‑149 (MQ=255)
                                                tCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                                <  1:766056/149‑1 (MQ=255)
                                                tCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                                <  1:800817/149‑1 (MQ=255)
                                                tCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                                <  1:3359042/149‑1 (MQ=255)
                                                tCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                                >  1:4395846/1‑149 (MQ=255)
                                                tCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                                <  1:425045/149‑1 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAg                                                                                               <  1:908003/149‑1 (MQ=255)
                                                  cccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                              <  2:4119203/149‑1 (MQ=255)
                                                    cGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAg                                                                                            <  1:2613836/149‑1 (MQ=255)
                                                       tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCAc                                                                                         <  1:2251095/149‑1 (MQ=255)
                                                             gCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcgcg                                                                                   >  1:3792855/1‑149 (MQ=255)
                                                             gCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcgcg                                                                                   <  1:1225152/149‑1 (MQ=255)
                                                              cGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCgg                                                                                  <  2:1623124/149‑1 (MQ=255)
                                                                tCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTc                                                                                <  1:31299/149‑1 (MQ=255)
                                                                 ctgctgctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCt                                                                               <  2:3036665/149‑1 (MQ=255)
                                                                       ctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCgtgt                                                                         <  1:3859144/149‑1 (MQ=255)
                                                                          tCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTgg                                                                      <  2:2157182/149‑1 (MQ=255)
                                                                            tGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCg                                                                    >  1:4122159/1‑149 (MQ=255)
                                                                              aaGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAc                                                                  >  1:2987078/1‑149 (MQ=255)
                                                                               aGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAcc                                                                 <  1:3668232/149‑1 (MQ=255)
                                                                                   tGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtttt                                                             >  1:4013435/1‑149 (MQ=255)
                                                                                   tGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtttt                                                             <  2:4275360/149‑1 (MQ=255)
                                                                                      ccGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGa                                                          <  1:1869706/149‑1 (MQ=255)
                                                                                       cGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGAc                                                         <  2:1591503/149‑1 (MQ=255)
                                                                                          tGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACaaa                                                      <  1:3441878/149‑1 (MQ=255)
                                                                                            aTCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGa                                                    >  2:1540003/1‑149 (MQ=255)
                                                                                              ccTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaaa                                                  <  2:560727/149‑1 (MQ=255)
                                                                                                  gATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                              >  1:1699863/1‑149 (MQ=255)
                                                                                                  gATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTATTCGACAAAGAAACCtt                                              >  1:4088746/1‑149 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTc                                             <  1:3963086/149‑1 (MQ=255)
                                                                                                     tATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGc                                           <  2:3107421/149‑1 (MQ=255)
                                                                                                      aTCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCt                                          <  1:3168755/149‑1 (MQ=255)
                                                                                                      aTCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCt                                          <  2:4608871/149‑1 (MQ=255)
                                                                                                      aTCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCt                                          <  2:1144652/149‑1 (MQ=255)
                                                                                                      aTCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCt                                          >  2:3978486/1‑149 (MQ=255)
                                                                                                      aTCACGTTGCGCGGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCt                                          <  2:4012939/149‑1 (MQ=255)
                                                                                                         aCGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaa                                       >  2:3331329/1‑149 (MQ=255)
                                                                                                           gTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaa                                     >  1:4075204/1‑149 (MQ=255)
                                                                                                            ttGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaaa                                    >  2:2849412/1‑149 (MQ=255)
                                                                                                            ttGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaaa                                    <  2:959135/149‑1 (MQ=255)
                                                                                                            ttGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaaaaa                                    >  2:1919586/1‑149 (MQ=255)
                                                                                                              gcgcTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACt                                  <  1:2439048/149‑1 (MQ=255)
                                                                                                                 cTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGaa                               <  2:259058/149‑1 (MQ=255)
                                                                                                                   ggATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaag                             >  1:505159/1‑149 (MQ=255)
                                                                                                                    gATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaaga                            >  2:2452518/1‑149 (MQ=255)
                                                                                                                     aTAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGTTATCGGTCCTGCTTATGTAGATAAAGTAGCACGTCGCGGTCTGAGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagac                           <  1:68617/149‑2 (MQ=255)
                                                                                                                      tAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaag                          <  1:2312141/149‑1 (MQ=255)
                                                                                                                        aCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTg                        >  1:1416178/1‑149 (MQ=255)
                                                                                                                         cgcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGa                       <  1:1941080/149‑1 (MQ=255)
                                                                                                                          gcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGAt                      <  1:1574018/149‑1 (MQ=255)
                                                                                                                          gcgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGAt                      <  1:1727338/149‑1 (MQ=255)
                                                                                                                           cgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATg                     <  1:635614/149‑1 (MQ=255)
                                                                                                                           cgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATg                     <  1:99672/149‑1 (MQ=255)
                                                                                                                            gcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATgg                    >  2:4242800/1‑149 (MQ=255)
                                                                                                                            gcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATgg                    <  1:4138251/149‑1 (MQ=255)
                                                                                                                               tGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAat                 <  2:1853415/149‑1 (MQ=255)
                                                                                                                               tGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAat                 >  1:123898/1‑149 (MQ=255)
                                                                                                                                gagaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAata                <  2:869878/149‑1 (MQ=255)
                                                                                                                                  gaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACTTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAACTTAGCTGAAAAACTGAAATAAGTGATGGAATATc              >  2:447319/1‑149 (MQ=255)
                                                                                                                                   aaaGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATca             >  2:2805098/1‑149 (MQ=255)
                                                                                                                                    aaGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcac            >  1:1601847/1‑149 (MQ=255)
                                                                                                                                     aGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca           >  2:381028/1‑149 (MQ=255)
                                                                                                                                     aGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca           <  1:2756999/149‑1 (MQ=255)
                                                                                                                                     aGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca           >  1:2122140/1‑149 (MQ=255)
                                                                                                                                           tGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACACCTTcc     <  2:3528823/149‑1 (MQ=255)
                                                                                                                                           tGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTcc     <  2:1421276/149‑1 (MQ=255)
                                                                                                                                           tGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTcc     >  1:2041200/1‑149 (MQ=255)
                                                                                                                                            ggCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCa    <  2:1599253/149‑1 (MQ=255)
                                                                                                                                            ggCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCa    <  1:1589936/149‑1 (MQ=255)
                                                                                                                                             gcgcTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAg   <  1:3284121/149‑1 (MQ=255)
                                                                                                                                             gcgcTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAg   >  1:825187/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCCAGT  >  CP009273/4394730‑4395020

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