Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,902,302 C→A 30.2% D242Y (GAT→TAT)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,902,3020CA30.2% 346.7 / 138.0 189D242Y (GAT→TAT) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (75/57);  new base A (32/25);  total (107/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 189 total aligned reads displayed.

GACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGA  >  CP009273/2902156‑2902446
                                                                                                                                                  |                                                                                                                                                
gACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCg                                                                                                                                                 <  1:2046253/148‑1 (MQ=255)
 aCTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGa                                                                                                                                                <  1:1007490/148‑1 (MQ=255)
 aCTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGa                                                                                                                                                <  1:2185014/148‑1 (MQ=255)
  cTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGAc                                                                                                                                               >  2:1387034/1‑148 (MQ=255)
   ttCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACg                                                                                                                                              >  2:1186925/1‑148 (MQ=255)
   ttCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACg                                                                                                                                              >  2:3438371/1‑148 (MQ=255)
     cACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTc                                                                                                                                            >  2:1316396/1‑148 (MQ=255)
        ttACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCttt                                                                                                                                         <  1:2787400/148‑1 (MQ=255)
        ttACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCttt                                                                                                                                         >  1:1977097/1‑148 (MQ=255)
           ccGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCag                                                                                                                                      >  1:463554/1‑148 (MQ=255)
            cGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCaga                                                                                                                                     <  1:763197/148‑1 (MQ=255)
            cGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATAGTCCAGCCCCTGAGATTTCAACAGGCCCGGAAATTTATAGATGGAATCGACGTCTTTCaca                                                                                                                                     >  1:1816717/1‑146 (MQ=255)
               gTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGaa                                                                                                                                  <  2:2297016/148‑1 (MQ=255)
               gTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGaa                                                                                                                                  <  1:2839560/148‑1 (MQ=255)
                    cttcttCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAAc                                                                                                                             >  1:1572613/1‑148 (MQ=255)
                       cttcGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGc                                                                                                                          <  2:3497800/148‑1 (MQ=255)
                        ttcGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGTAATAGACGTCTTTCAGAGAAATAACCGCt                                                                                                                         >  2:1448512/1‑148 (MQ=255)
                          cGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCttt                                                                                                                       >  2:2873280/1‑148 (MQ=255)
                            aaGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCttttt                                                                                                                     >  1:307778/1‑148 (MQ=255)
                            aaGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCttttt                                                                                                                     >  1:237811/1‑148 (MQ=255)
                              gATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTcc                                                                                                                   >  1:2045960/1‑148 (MQ=255)
                              gATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTcc                                                                                                                   >  2:3002860/1‑148 (MQ=255)
                              gATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTcc                                                                                                                   >  2:2385475/1‑148 (MQ=255)
                               aTAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCCGAGAAATAACCGCTTTTTCCg                                                                                                                  >  2:2353592/1‑148 (MQ=255)
                                 aaCCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGa                                                                                                                >  2:1009741/1‑148 (MQ=255)
                                        tCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTAc                                                                                                         <  1:1814437/148‑1 (MQ=255)
                                           cATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGa                                                                                                      <  2:3139470/148‑1 (MQ=255)
                                              tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGCTTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTACAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   >  2:2435844/1‑148 (MQ=255)
                                              tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   <  1:2539882/148‑1 (MQ=255)
                                              tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   <  1:107600/148‑1 (MQ=255)
                                              tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   <  1:2358060/148‑1 (MQ=255)
                                              tCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACa                                                                                                   <  1:2357674/148‑1 (MQ=255)
                                               cGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACaa                                                                                                  <  2:839518/148‑1 (MQ=255)
                                                ggACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAAt                                                                                                 >  2:3256952/1‑148 (MQ=255)
                                                ggACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGACCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCATTCAGAGAAATAACCGCTTTTTCCGGACCATTACAGAACAAt                                                                                                 >  2:3259046/1‑148 (MQ=255)
                                                    aGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCaa                                                                                             <  2:2691367/148‑1 (MQ=255)
                                                        tCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCtt                                                                                         <  2:247380/148‑1 (MQ=255)
                                                          gCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCg                                                                                       <  2:2140730/148‑1 (MQ=255)
                                                          gCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCg                                                                                       <  2:2071764/148‑1 (MQ=255)
                                                          gCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCg                                                                                       >  2:532596/1‑148 (MQ=255)
                                                            ttCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCa                                                                                     >  2:1160618/1‑148 (MQ=255)
                                                             tCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAc                                                                                    >  1:1436778/1‑148 (MQ=255)
                                                               cGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgt                                                                                  <  2:1781991/148‑1 (MQ=255)
                                                                  gCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACCGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcg                                                                               >  1:3097564/1‑148 (MQ=255)
                                                                  gCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcg                                                                               <  2:1653133/148‑1 (MQ=255)
                                                                  gCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcg                                                                               <  2:1312914/148‑1 (MQ=255)
                                                                    aGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCAcgttcgtt                                                                             <  1:1052819/148‑1 (MQ=255)
                                                                     gTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACgttcgttc                                                                            <  2:613008/148‑1 (MQ=255)
                                                                           gCTGAATCGTTTACAAATATAATCTTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTAAGATAAATAACCGCTTTTTACTGAACATTAAAGAACAATGCAATATTCGCACGTTCGTTCGCAGGa                                                                      >  1:1361766/1‑148 (MQ=255)
                                                                            cTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGaa                                                                     >  1:2853538/1‑148 (MQ=255)
                                                                                    tttACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGa                                                             >  2:986855/1‑148 (MQ=255)
                                                                                    tttACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACAGTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGa                                                             <  1:3099581/148‑1 (MQ=255)
                                                                                    tttACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTGAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGa                                                             <  2:2289816/148‑1 (MQ=255)
                                                                                    tttACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGa                                                             >  1:1528349/1‑148 (MQ=255)
                                                                                     ttACAAATATACTCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGAt                                                            <  2:396788/148‑1 (MQ=255)
                                                                                     ttACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAAAATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGAt                                                            >  1:1317729/1‑148 (MQ=255)
                                                                                      tACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATc                                                           <  1:2967688/148‑1 (MQ=255)
                                                                                        cAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTg                                                         >  2:69779/1‑148 (MQ=255)
                                                                                          aaTATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGaa                                                       >  1:3190516/1‑148 (MQ=255)
                                                                                          aaTATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGaa                                                       >  1:3190190/1‑148 (MQ=255)
                                                                                            tataATCGTCAAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACg                                                     >  1:852895/1‑148 (MQ=255)
                                                                                             ataATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGa                                                    <  2:1572613/148‑1 (MQ=255)
                                                                                              taATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGAc                                                   <  1:3189910/148‑1 (MQ=255)
                                                                                                 tCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACaaa                                                <  1:3320842/148‑1 (MQ=255)
                                                                                                 tCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACaaa                                                >  2:3284340/1‑148 (MQ=255)
                                                                                                 tCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTATATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCACTCATATCACGTTCGTTCGACGGAACAGCGCGATCTGAACGACaaa                                                >  1:592889/1‑148 (MQ=255)
                                                                                                     ccAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAg                                            >  1:3372035/1‑148 (MQ=255)
                                                                                                      cAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAgg                                           <  2:2617581/148‑1 (MQ=255)
                                                                                                      cAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAgg                                           <  2:1977097/148‑1 (MQ=255)
                                                                                                      cAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAgg                                           <  1:805876/148‑1 (MQ=255)
                                                                                                       aGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGa                                          <  1:2911833/148‑1 (MQ=255)
                                                                                                       aGCCCCTGACAATTCAACAGGCCCGGAATTTAATAGATGTAATCGACGACTTTCAGAGAAATAACCGCTTTTTCCTGAACAGTACAGAACAATGCAATCTTCGCACGATCGATAGCCGGAACAGCGAGATCTGAACGACAAATCAGGa                                          >  2:24542/1‑148 (MQ=255)
                                                                                                        gCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGAt                                         >  2:3426826/1‑148 (MQ=255)
                                                                                                              gagaTTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAgg                                   >  1:1722801/1‑148 (MQ=255)
                                                                                                              gagaTTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAgg                                   >  1:1144299/1‑148 (MQ=255)
                                                                                                              gagaTTTCAACAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAgg                                   >  2:283998/1‑148 (MQ=255)
                                                                                                               agaTTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGc                                  <  1:2610394/148‑1 (MQ=255)
                                                                                                               agaTTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGc                                  >  1:458635/1‑148 (MQ=255)
                                                                                                                  tttCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTgg                               <  1:1316396/148‑1 (MQ=255)
                                                                                                                  tttCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTgg                               >  1:126110/1‑148 (MQ=255)
                                                                                                                     cAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATa                            >  1:3243017/1‑148 (MQ=255)
                                                                                                                       aCAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATAcc                          >  2:2663822/1‑148 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCg                         <  1:535018/148‑1 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCg                         <  1:2323082/148‑1 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCg                         >  2:1052640/1‑148 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCg                         >  2:939792/1‑148 (MQ=255)
                                                                                                                        cAGGCCCGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCg                         >  1:1648755/1‑148 (MQ=255)
                                                                                                                          ggCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGAt                       <  2:126110/148‑1 (MQ=255)
                                                                                                                          ggCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGAt                       <  2:223577/148‑1 (MQ=255)
                                                                                                                           gCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATg                      >  1:3539801/1‑148 (MQ=255)
                                                                                                                             ccGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGga                    <  1:1340573/148‑1 (MQ=255)
                                                                                                                             ccGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGga                    <  2:210724/148‑1 (MQ=255)
                                                                                                                              cGGAATTTTATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGgag                   >  2:1766967/1‑148 (MQ=255)
                                                                                                                                gAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGagag                 >  1:369517/1‑148 (MQ=255)
                                                                                                                                   ttttATAGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGagcag              <  1:3438371/148‑1 (MQ=255)
                                                                                                                                    tttATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGagcagc             >  1:1832308/1‑148 (MQ=255)
                                                                                                                                        taGATGGAATAGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCtt         >  2:1642403/1‑148 (MQ=255)
                                                                                                                                           aTGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTa      <  2:3371130/148‑1 (MQ=255)
                                                                                                                                            tGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTAc     >  1:2695366/1‑148 (MQ=255)
                                                                                                                                            tGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTAc     >  2:1707620/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACAGATTCGCTTCCGGGCAGTTTAAGCTGAATCGTTTACAAATATAATCGTCCAGCCCCTGAGATTTCAACAGGCCCGGAATTTTATAGATGGAATCGACGTCTTTCAGAGAAATAACCGCTTTTTCCGGAACATTACAGAACAATGCAATCTTCGCACGTTCGTTCGCCGGAACAGCGCGATCTGAACGACAAATCAGGATGTCAGGCTGGATACCGATGGAGAGCAGCTCTTTTACAGA  >  CP009273/2902156‑2902446

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: