Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,336,599 T→C 14.5% L141P (CTT→CCT)  pyrF → orotidine‑5'‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,336,5990TC14.5% 516.5 / 78.8 200L141P (CTT→CCT) pyrForotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (93/78);  new base C (20/9);  total (113/87)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

Only 100 of 202 total aligned reads displayed.

CGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCGGG  >  CP009273/1336453‑1336746
                                                                                                                                                  |                                                                                                                                                   
cGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCtt                                                                                                                                                    <  2:493237/148‑1 (MQ=255)
   gtgGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGc                                                                                                                                                 <  2:3526186/148‑1 (MQ=255)
            gAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTg                                                                                                                                        <  1:1694167/148‑1 (MQ=255)
             aaTGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGt                                                                                                                                       >  1:2446683/1‑148 (MQ=255)
              tttttCTTGCGTCTGGTGGGGCGTGTTTGTTGACCGCTGCTCGTGTGGCACTTTTTGAGTTTGGCAATGATGCACGGCTTTTGTTTGTTGTTACAGTTTTGACCTGATTGGATGCAAGCGACCTGGTCGATCTTGGCATGACACTGTc                                                                                                                                      <  1:349018/145‑1 (MQ=255)
              aTGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTc                                                                                                                                      <  1:2142418/148‑1 (MQ=255)
               tGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCa                                                                                                                                     >  1:2807322/1‑148 (MQ=255)
                gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCAc                                                                                                                                    >  1:1315212/1‑148 (MQ=255)
                gTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCAc                                                                                                                                    >  2:2204773/1‑148 (MQ=255)
                  tCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCt                                                                                                                                  <  2:103220/148‑1 (MQ=255)
                    aTGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGc                                                                                                                                >  2:886688/1‑148 (MQ=255)
                       ccTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCTCTGGTTCGGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTTTTGACAAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGa                                                                                                                             <  1:2597389/148‑1 (MQ=255)
                       ccTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCACGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGa                                                                                                                             >  1:315316/1‑148 (MQ=255)
                        ctctGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGAt                                                                                                                            <  1:3411025/148‑1 (MQ=255)
                         tctGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGAtt                                                                                                                           <  2:2033116/148‑1 (MQ=255)
                               ggggCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAg                                                                                                                     <  2:924767/148‑1 (MQ=255)
                                 ggCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGaa                                                                                                                   <  1:1946845/148‑1 (MQ=255)
                                    gcgTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGt                                                                                                                <  1:2529447/148‑1 (MQ=255)
                                     cgTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTTGCATGACACTGTCACCTGCAGATTATGCAGAACGTc                                                                                                               <  2:138707/148‑1 (MQ=255)
                                     cgTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTGGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTc                                                                                                               <  1:2524424/148‑1 (MQ=255)
                                       tATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTg                                                                                                             >  2:1316553/1‑148 (MQ=255)
                                       tATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTg                                                                                                             <  1:1719456/148‑1 (MQ=255)
                                       tATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTg                                                                                                             >  2:536151/1‑148 (MQ=255)
                                        atgatgACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTgg                                                                                                            >  2:1754499/1‑148 (MQ=255)
                                          gatgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcg                                                                                                          >  1:898110/1‑148 (MQ=255)
                                           atgaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTggcgg                                                                                                         >  2:2069066/1‑148 (MQ=255)
                                             gaCCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCa                                                                                                       <  1:2320438/148‑1 (MQ=255)
                                               ccGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACt                                                                                                     >  1:2884948/1‑148 (MQ=255)
                                               ccGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACt                                                                                                     >  1:2884958/1‑148 (MQ=255)
                                                cGCAGCGCGTGCGGCACTGTTTCCGTTTGGCAAAGATGCACCGCTTTTGATTTCTTTGACATTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTg                                                                                                    <  1:2984956/148‑1 (MQ=255)
                                                cGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTg                                                                                                    <  1:714888/148‑1 (MQ=255)
                                                 gCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGa                                                                                                   <  1:1468309/148‑1 (MQ=255)
                                                   aGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACg                                                                                                 <  2:574633/148‑1 (MQ=255)
                                                    gcgcGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGc                                                                                                >  1:2945436/1‑148 (MQ=255)
                                                        gTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCaaaa                                                                                            <  1:420612/148‑1 (MQ=255)
                                                         tGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCaaaaa                                                                                           <  1:1535836/148‑1 (MQ=255)
                                                         tGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCaaaaa                                                                                           <  2:2088012/148‑1 (MQ=255)
                                                            ggCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAACtgt                                                                                        >  2:2519003/1‑148 (MQ=255)
                                                              cACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTgg                                                                                      >  2:3179243/1‑148 (MQ=255)
                                                                cTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGcc                                                                                    >  2:372000/1‑148 (MQ=255)
                                                                cTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGcc                                                                                    <  1:3024511/148‑1 (MQ=255)
                                                                cTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGcc                                                                                    >  2:197716/1‑148 (MQ=255)
                                                                  ggTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCtt                                                                                  >  2:1092373/1‑148 (MQ=255)
                                                                  ggTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGGCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCtt                                                                                  <  2:154814/148‑1 (MQ=255)
                                                                     tCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAt                                                                               >  2:634660/1‑148 (MQ=255)
                                                                      ccGTTTGGCACAGTTGCACCGATTTTGATTGCTGGGACAGTGTTGAGCAGCATGGAAGCCAGCGACCTGGTGGATCTTGGCATGACACTGTCACCTGCTGATTATGCAGAAAGTCTGGCGGCACTGACGCAAAAATTTGGCCTTGAtg                                                                              <  2:544592/148‑1 (MQ=255)
                                                                        gTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggt                                                                            >  2:2369643/1‑148 (MQ=255)
                                                                         tttGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtg                                                                           <  1:1316543/148‑1 (MQ=255)
                                                                              cAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCAGCAGATTATTCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGAtggtgtggtg                                                                      >  1:3277843/1‑148 (MQ=255)
                                                                               aaaGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATggtgtggtgt                                                                     >  2:2861414/1‑148 (MQ=255)
                                                                                taGATGCACCGCTTTTGATTGCTGTGACTGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGgtgtggtgtg                                                                    <  2:1500678/147‑1 (MQ=255)
                                                                                     gCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTg                                                               >  2:3228864/1‑148 (MQ=255)
                                                                                     gCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTg                                                               >  2:2619333/1‑148 (MQ=255)
                                                                                     gCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTg                                                               >  2:1425699/1‑148 (MQ=255)
                                                                                       aCCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCt                                                             >  2:3203554/1‑148 (MQ=255)
                                                                                          gCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAg                                                          >  2:1074430/1‑148 (MQ=255)
                                                                                           cTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAgg                                                         <  2:3340102/148‑1 (MQ=255)
                                                                                             tttGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGaa                                                       >  2:63350/1‑148 (MQ=255)
                                                                                                   tGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTTACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCtgtg                                                 >  1:3375797/1‑148 (MQ=255)
                                                                                                     cTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTgcg                                               >  1:3346072/1‑148 (MQ=255)
                                                                                                       gtgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCt                                             <  1:436701/148‑1 (MQ=255)
                                                                                                        tgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCtt                                            <  1:606445/148‑1 (MQ=255)
                                                                                                        tgACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCtt                                            <  1:610088/148‑1 (MQ=255)
                                                                                                         gACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCttt                                           >  1:1364381/1‑148 (MQ=255)
                                                                                                         gACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCttt                                           <  1:1018846/148‑1 (MQ=255)
                                                                                                              tgtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACa                                      >  1:3011733/1‑148 (MQ=255)
                                                                                                               gtTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAg                                     >  1:3081850/1‑148 (MQ=255)
                                                                                                                 tGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGt                                   >  2:1068918/1‑148 (MQ=255)
                                                                                                                 tGACCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGt                                   >  2:1068909/1‑148 (MQ=255)
                                                                                                                  gACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTa                                  <  1:1235971/148‑1 (MQ=255)
                                                                                                                  gACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTa                                  <  1:853216/148‑1 (MQ=255)
                                                                                                                   aCCAGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTAt                                 >  2:1498441/1‑148 (MQ=255)
                                                                                                                      aGCATGGAAGCCAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCg                              <  2:2446683/148‑1 (MQ=255)
                                                                                                                       gCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTTATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCgg                             <  2:2670818/148‑1 (MQ=255)
                                                                                                                       gCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCgg                             <  1:226586/148‑1 (MQ=255)
                                                                                                                       gCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCgg                             <  1:226430/148‑1 (MQ=255)
                                                                                                                        cATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGt                            <  1:1754499/148‑1 (MQ=255)
                                                                                                                          tGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCa                          <  1:404677/148‑1 (MQ=255)
                                                                                                                            gAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAgg                        >  2:3108114/1‑148 (MQ=255)
                                                                                                                               gCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGt                     <  1:2609803/148‑1 (MQ=255)
                                                                                                                                ccAGCGACCTGGTCGATCCTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGtt                    <  2:1175918/148‑1 (MQ=255)
                                                                                                                                 cAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTc                   >  1:1162641/1‑148 (MQ=255)
                                                                                                                                 cAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTc                   >  1:334463/1‑148 (MQ=255)
                                                                                                                                   gCGACCTGGTCGATCTTGGCATTACACTGTAACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaa                 >  1:2272830/1‑148 (MQ=255)
                                                                                                                                   gCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaa                 <  2:1928164/148‑1 (MQ=255)
                                                                                                                                   gCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaa                 >  2:1792370/1‑148 (MQ=255)
                                                                                                                                    cGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa                >  2:842639/1‑148 (MQ=255)
                                                                                                                                    cGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCaaa                >  1:2202771/1‑148 (MQ=255)
                                                                                                                                     gACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAAc               <  2:2751505/148‑1 (MQ=255)
                                                                                                                                     gACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAAc               <  1:765903/148‑1 (MQ=255)
                                                                                                                                       ccTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTg             >  1:1155859/1‑148 (MQ=255)
                                                                                                                                        cTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTgg            <  2:1262007/148‑1 (MQ=255)
                                                                                                                                        cTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTgg            >  2:2287904/1‑148 (MQ=255)
                                                                                                                                           gTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTa         <  2:78169/148‑1 (MQ=255)
                                                                                                                                           gTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTa         >  2:2158859/1‑148 (MQ=255)
                                                                                                                                                tCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCg    >  1:1479376/1‑148 (MQ=255)
                                                                                                                                                 cTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCgg   >  1:1248366/1‑148 (MQ=255)
                                                                                                                                                 cTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCgg   >  1:2383876/1‑148 (MQ=255)
                                                                                                                                                 cTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCgg   >  2:884131/1‑148 (MQ=255)
                                                                                                                                                 cTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCgg   >  2:884981/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGCATGGAAGCCAGCGACCTGGTCGATCTTGGCATGACACTGTCACCTGCAGATTATGCAGAACGTCTGGCGGCACTGACGCAAAAATGTGGCCTTGATGGTGTGGTGTGTTCTGCTCAGGAAGCTGTGCGCTTTAAACAGGTATTCGGTCAGGAGTTCAAACTGGTTACGCCGGG  >  CP009273/1336453‑1336746

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: