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,8750GT97.8% 719.8 / 5.8 228A124A (GCG→GCTpurAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (104/119);  total (108/120)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-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 228 total aligned reads displayed.

CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCC  >  CP009273/4394730‑4395017
                                                                                                                                                 |                                                                                                                                              
ctcCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTaa                                                                                                                                              >  1:2422361/1‑148 (MQ=255)
  ccGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAg                                                                                                                                            <  1:3505938/148‑1 (MQ=255)
   cGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGc                                                                                                                                           >  2:1996330/1‑148 (MQ=255)
      ccGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGAt                                                                                                                                        <  2:724713/148‑1 (MQ=255)
      ccGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGAt                                                                                                                                        <  1:3301417/148‑1 (MQ=255)
         cgcTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCgg                                                                                                                                     >  2:1897265/1‑148 (MQ=255)
           cTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGca                                                                                                                                   <  1:852895/148‑1 (MQ=255)
             gatgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacc                                                                                                                                 >  2:3132346/1‑148 (MQ=255)
              atgaaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGcacca                                                                                                                                <  2:3285547/148‑1 (MQ=255)
                 aaagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCg                                                                                                                             <  2:686989/148‑1 (MQ=255)
                  aagagatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCgg                                                                                                                            >  1:2960052/1‑148 (MQ=255)
                      gatgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGt                                                                                                                        <  2:668161/148‑1 (MQ=255)
                       atgaaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTg                                                                                                                       >  1:595303/1‑148 (MQ=255)
                          aaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTa                                                                                                                    >  1:1369499/1‑148 (MQ=255)
                          aaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTa                                                                                                                    >  2:3015749/1‑148 (MQ=255)
                          aaagaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTa                                                                                                                    >  1:1724327/1‑148 (MQ=255)
                             gaACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCg                                                                                                                 <  2:2901378/148‑1 (MQ=255)
                              aaCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCgg                                                                                                                >  1:3176573/1‑148 (MQ=255)
                               aCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCggg                                                                                                               <  1:1197690/148‑1 (MQ=255)
                               aCTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCggg                                                                                                               <  1:625085/148‑1 (MQ=255)
                                       aCCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTa                                                                                                       <  2:1054422/148‑1 (MQ=255)
                                          gTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGa                                                                                                    <  2:1145604/148‑1 (MQ=255)
                                           tGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGaa                                                                                                   >  2:127730/1‑148 (MQ=255)
                                            ggCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAg                                                                                                  >  2:2927644/1‑148 (MQ=255)
                                              cATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGAt                                                                                                <  2:3183872/148‑1 (MQ=255)
                                               aTCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATa                                                                                               <  1:1897265/148‑1 (MQ=255)
                                                 ccccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATaaa                                                                                             >  2:1590615/1‑148 (MQ=255)
                                                   ccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                           <  2:1013130/148‑1 (MQ=255)
                                                   ccGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGt                                                                                           <  1:623002/148‑1 (MQ=255)
                                                    cGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATc                                                                                                                                       >  1:371589/1‑103 (MQ=255)
                                                      ttCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCTTGAGAAAGCGCGTGTCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGTGCCTGCTTATGAAGATAAAGTAGc                                                                                        <  2:2535241/148‑1 (MQ=255)
                                                      ttCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGc                                                                                        <  2:1452382/148‑1 (MQ=255)
                                                       tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCa                                                                                       >  2:2845009/1‑148 (MQ=255)
                                                       tCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCa                                                                                       >  2:636818/1‑148 (MQ=255)
                                                         gTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACg                                                                                     <  1:309184/148‑1 (MQ=255)
                                                           gAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTc                                                                                   >  1:1073466/1‑148 (MQ=255)
                                                              cGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcgcg                                                                                <  1:697256/148‑1 (MQ=255)
                                                              cGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTcgcg                                                                                <  1:2748139/148‑1 (MQ=255)
                                                                tCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGt                                                                              >  1:3334195/1‑148 (MQ=255)
                                                                tCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGt                                                                              <  1:234063/148‑1 (MQ=255)
                                                                 ctgctgctgTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTc                                                                             >  1:2358360/1‑148 (MQ=255)
                                                                         gTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGtt                                                                     <  2:1343355/148‑1 (MQ=255)
                                                                             gAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCg                                                                 <  1:2745944/148‑1 (MQ=255)
                                                                              aaGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGa                                                                >  1:1116673/1‑148 (MQ=255)
                                                                                gCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAcc                                                              >  2:1758416/1‑148 (MQ=255)
                                                                                gCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAcc                                                              <  1:3095249/148‑1 (MQ=255)
                                                                                gCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGAcc                                                              <  1:2884222/148‑1 (MQ=255)
                                                                                  aTGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCtt                                                            >  1:1955729/1‑148 (MQ=255)
                                                                                   tGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCttt                                                           <  2:8221/148‑1 (MQ=255)
                                                                                   tGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCttt                                                           >  1:741151/1‑148 (MQ=255)
                                                                                   tGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCttt                                                           >  1:533825/1‑148 (MQ=255)
                                                                                      ccGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCg                                                        <  2:2298435/148‑1 (MQ=255)
                                                                                         cTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACa                                                     >  2:2748081/1‑148 (MQ=255)
                                                                                           gATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACaaa                                                   <  1:1814497/148‑1 (MQ=255)
                                                                                            aTCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAg                                                  >  1:496505/1‑148 (MQ=255)
                                                                                            aTCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAg                                                  >  1:1231783/1‑148 (MQ=255)
                                                                                             tCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACAGGTCGTGGTATCGGGCCAGCTTATGAAGATAAAGTAGCACGTCGAGCTCTGCGTGTTGGCGACCTTTTAGACAAAGa                                                 >  2:2585373/1‑148 (MQ=255)
                                                                                              ccTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaa                                                >  1:2809346/1‑148 (MQ=255)
                                                                                              ccTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaa                                                <  2:2237073/148‑1 (MQ=255)
                                                                                              ccTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaa                                                <  2:7985/148‑1 (MQ=255)
                                                                                               cTTGATTCTCACGTGGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGaaa                                               <  2:2504984/148‑1 (MQ=255)
                                                                                               cTTGATTATCACGTTGCTCTTGATATCGCGCGTGAGATAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCAAGTCGCGGTATGCGTGTTGGCGACCTTTTCGACTAAGtaa                                               <  2:1376162/148‑1 (MQ=255)
                                                                                                ttGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAc                                              >  2:3077313/1‑148 (MQ=255)
                                                                                                ttGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAc                                              <  2:1723414/148‑1 (MQ=255)
                                                                                                ttGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAc                                              >  1:2749543/1‑148 (MQ=255)
                                                                                                 tGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAAcc                                             <  1:518508/148‑1 (MQ=255)
                                                                                                  gATTATCACGTTGCTCTGGATAAGGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGTCCTGCTTATGAAGATAATGTAGCACGTCGCGGTCTGCGTGTTGTCGACCTTTTCGACAAAGAAACCt                                            <  1:1648154/148‑1 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                           <  2:588471/148‑1 (MQ=255)
                                                                                                   aTTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCtt                                           <  1:1575825/148‑1 (MQ=255)
                                                                                                    ttATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTc                                          >  1:1594048/1‑148 (MQ=255)
                                                                                                     tATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCg                                         <  1:2667670/148‑1 (MQ=255)
                                                                                                       tCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCt                                       <  2:2422361/148‑1 (MQ=255)
                                                                                                        cACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACTGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTg                                      >  1:31317/1‑148 (MQ=255)
                                                                                                          cGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGaa                                    <  2:894844/148‑1 (MQ=255)
                                                                                                                gcTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTg                              <  1:2134580/148‑1 (MQ=255)
                                                                                                                gcTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTg                              <  1:3204499/148‑1 (MQ=255)
                                                                                                                gcTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTg                              <  1:338939/148‑1 (MQ=255)
                                                                                                                gcTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTg                              <  2:457838/148‑1 (MQ=255)
                                                                                                                gcTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTg                              <  1:2586821/148‑1 (MQ=255)
                                                                                                                 cTGGATAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGa                             >  1:2310399/1‑148 (MQ=255)
                                                                                                                     aTAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaaga                         <  2:1073466/148‑1 (MQ=255)
                                                                                                                     aTAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaaga                         <  2:3532171/148‑1 (MQ=255)
                                                                                                                      tAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaa                        <  1:1885524/148‑1 (MQ=255)
                                                                                                                      tAACGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAaagaa                        <  2:3179023/148‑1 (MQ=255)
                                                                                                                        aCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                      <  2:1699324/148‑1 (MQ=255)
                                                                                                                        aCGCGCGTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGt                      >  1:1039193/1‑148 (MQ=255)
                                                                                                                           cgcgTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGAt                   <  1:3155654/148‑1 (MQ=255)
                                                                                                                              gTGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGa                <  1:809276/148‑1 (MQ=255)
                                                                                                                               tGAGAAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGaa               >  1:344557/1‑148 (MQ=255)
                                                                                                                                 agaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAata             >  1:1244431/1‑148 (MQ=255)
                                                                                                                                 agaAAGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAata             <  1:1835893/148‑1 (MQ=255)
                                                                                                                                   aaaGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATc           <  1:127730/148‑1 (MQ=255)
                                                                                                                                   aaaGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATc           >  1:2870100/1‑148 (MQ=255)
                                                                                                                                   aaaGCGCGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATc           >  1:1754826/1‑148 (MQ=255)
                                                                                                                                      gcgcGTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATcaca        <  1:250861/148‑1 (MQ=255)
                                                                                                                                        gcgTGGCGCTAAAGAGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAAc      <  2:2946258/148‑1 (MQ=255)
                                                                                                                                         cgTGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGACTATCACAACt     >  2:1833970/1‑148 (MQ=255)
                                                                                                                                           tGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTc   <  2:3030562/148‑1 (MQ=255)
                                                                                                                                           tGGCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACATc   >  1:2842202/1‑148 (MQ=255)
                                                                                                                                            ggCGCTAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTcc  <  1:474658/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTCCGGCCGCGCTGATGAAAGAGATGAAAGAACTGGAAGACCGTGGCATCCCCGTTCGTGAGCGTCTGCTGCTGTCTGAAGCATGTCCGCTGATCCTTGATTATCACGTTGCGCTGGATAACGCGCGTGAGAAAGCGCGTGGCGCGAAAGCGATCGGCACCACCGGTCGTGGTATCGGGCCTGCTTATGAAGATAAAGTAGCACGTCGCGGTCTGCGTGTTGGCGACCTTTTCGACAAAGAAACCTTCGCTGAAAAACTGAAAGAAGTGATGGAATATCACAACTTCC  >  CP009273/4394730‑4395017

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: