Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,783,091 G→T 100% E134* (GAG→TAG)  aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,783,0910GT97.9% 606.5 / 3.3 195E134* (GAG→TAG) aroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (114/77);  total (117/78)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-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 195 total aligned reads displayed.

GGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTC  >  CP009273/1782949‑1783232
                                                                                                                                              |                                                                                                                                             
    aTCGTAATGCTCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACCGTCTGGAGCGGGCGCGCAAATTACTTTTACAGTTACATTAGcagtgct                                                                                                                                      >  2:3082554/1‑142 (MQ=255)
    aTCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCg                                                                                                                                      >  1:2869084/1‑148 (MQ=255)
      cGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTc                                                                                                                                    <  2:2781689/148‑1 (MQ=255)
       gTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTcc                                                                                                                                   <  2:1098096/148‑1 (MQ=255)
         aaTGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCa                                                                                                                                 <  1:1502184/148‑1 (MQ=255)
          aTGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTcccca                                                                                                                                >  2:2391103/1‑146 (MQ=255)
           tGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAAc                                                                                                                               >  1:374569/1‑148 (MQ=255)
           tGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAAc                                                                                                                               >  2:1412499/1‑148 (MQ=255)
           tGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAAc                                                                                                                               >  2:1894389/1‑148 (MQ=255)
           tGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAAc                                                                                                                               >  2:1894321/1‑148 (MQ=255)
                aCCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGa                                                                                                                          <  1:2089500/148‑1 (MQ=255)
                 ccTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGAc                                                                                                                         <  2:1883127/148‑1 (MQ=255)
                  cTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGAcc                                                                                                                        >  2:3096193/1‑148 (MQ=255)
                    aTTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGa                                                                                                                      >  1:215421/1‑148 (MQ=255)
                    aTTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGa                                                                                                                      >  2:1631033/1‑148 (MQ=255)
                          aaaaaCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCt                                                                                                                >  1:2469064/1‑148 (MQ=255)
                            aaaCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCg                                                                                                              >  2:682438/1‑148 (MQ=255)
                            aaaCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCg                                                                                                              >  1:3279602/1‑148 (MQ=255)
                             aaCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGa                                                                                                             <  1:1916159/148‑1 (MQ=255)
                             aaCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGa                                                                                                             >  1:914061/1‑148 (MQ=255)
                              aCCACGACCTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGat                                                                                                            <  1:2622722/148‑1 (MQ=255)
                              aCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGat                                                                                                            <  1:706403/148‑1 (MQ=255)
                              aCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGat                                                                                                            <  2:813610/148‑1 (MQ=255)
                               ccACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGata                                                                                                           >  1:1631464/1‑148 (MQ=255)
                                cccGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGatat                                                                                                          <  1:105663/146‑1 (MQ=255)
                                cACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGatat                                                                                                          <  1:230857/148‑1 (MQ=255)
                                    aaCTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTg                                                                                                      <  1:1681890/148‑1 (MQ=255)
                                    aaCTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTg                                                                                                      >  1:1290391/1‑148 (MQ=255)
                                    aaCTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTg                                                                                                      <  1:3576594/148‑1 (MQ=255)
                                      cTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGAc                                                                                                    >  2:1216479/1‑148 (MQ=255)
                                      cTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGAc                                                                                                    >  2:3795170/1‑148 (MQ=255)
                                       tgttgtCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGAcc                                                                                                   >  2:2683079/1‑148 (MQ=255)
                                         ttgtCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCgg                                                                                                 >  2:2070772/1‑148 (MQ=255)
                                            tCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCa                                                                                              <  2:3703242/148‑1 (MQ=255)
                                                 tGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTa                                                                                         <  1:3012149/148‑1 (MQ=255)
                                                  ggAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTAt                                                                                        <  1:1830381/148‑1 (MQ=255)
                                                  ggAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTAt                                                                                        <  2:1190183/148‑1 (MQ=255)
                                                  ggAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTAt                                                                                        <  1:3703162/148‑1 (MQ=255)
                                                  ggAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTAt                                                                                        <  2:1560807/148‑1 (MQ=255)
                                                     aaGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTCTTGc                                                                                     >  1:687150/1‑148 (MQ=255)
                                                     aaGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGc                                                                                     >  1:687137/1‑148 (MQ=255)
                                                      aGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCt                                                                                    >  2:2816107/1‑148 (MQ=255)
                                                        gACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGa                                                                                  >  2:1889239/1‑148 (MQ=255)
                                                        gACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGa                                                                                  >  2:2495850/1‑148 (MQ=255)
                                                        gACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGa                                                                                  <  2:3665882/148‑1 (MQ=255)
                                                        gACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGa                                                                                  <  2:1872100/148‑1 (MQ=255)
                                                         aCTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATCTGGTGACCGGTCAGTGTATTGCTGAt                                                                                 >  1:2757407/1‑148 (MQ=255)
                                                            aaTCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTa                                                                              >  1:424720/1‑148 (MQ=255)
                                                            aaTCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTa                                                                              >  2:1498695/1‑148 (MQ=255)
                                                            aaTCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTa                                                                              >  2:1498361/1‑148 (MQ=255)
                                                            aaTCTCCGATCAATATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTa                                                                              <  1:61517/148‑1 (MQ=255)
                                                              tctcCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAAt                                                                            >  1:1872106/1‑148 (MQ=255)
                                                               ctcCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATc                                                                           <  2:844050/148‑1 (MQ=255)
                                                               ctcCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATc                                                                           >  1:1189297/1‑148 (MQ=255)
                                                                tcCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATAGGGTGACCGGTCAGTTTATTGCTGAGTTAATCa                                                                          >  1:3175695/1‑148 (MQ=255)
                                                                    aTCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTg                                                                      >  1:2415900/1‑148 (MQ=255)
                                                                            ttAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGa                                                              <  1:350200/148‑1 (MQ=255)
                                                                                aCGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTgg                                                          >  1:2222085/1‑148 (MQ=255)
                                                                                 cGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGAGTGgc                                                         >  2:3037957/1‑148 (MQ=255)
                                                                                   gCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGCTATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGgcgc                                                       >  2:3376847/1‑148 (MQ=255)
                                                                                   gCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGgcgc                                                       >  2:3376869/1‑148 (MQ=255)
                                                                                    cAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCa                                                      >  2:3628469/1‑148 (MQ=255)
                                                                                      gCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCGTCGATATGGTGACCGGTCCGTTTATGGCTGATTTAATCAGTTGGGGGGCGATTGGCGCACg                                                    >  2:3168064/1‑148 (MQ=255)
                                                                                       cTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGt                                                   >  2:3716596/1‑148 (MQ=255)
                                                                                       cTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGt                                                   >  1:3099813/1‑148 (MQ=255)
                                                                                       cTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGt                                                   <  2:72696/148‑1 (MQ=255)
                                                                                         aTCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtac                                                 >  2:443977/1‑148 (MQ=255)
                                                                                         aTCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtac                                                 >  1:1115071/1‑148 (MQ=255)
                                                                                         aTCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtac                                                 >  2:2831864/1‑148 (MQ=255)
                                                                                          tCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtact                                                >  1:1968445/1‑148 (MQ=255)
                                                                                          tCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGGTCCTCGATCTGGTGACCGGTCAGTTTATTGCTGATTTACTCAGTTGGGGCGCGATTGGCGCACGtact                                                >  1:1964688/1‑148 (MQ=255)
                                                                                           cGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtacta                                               <  2:2163978/148‑1 (MQ=255)
                                                                                           cGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtacta                                               <  2:1487269/148‑1 (MQ=255)
                                                                                               tAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGa                                           >  1:1546294/1‑148 (MQ=255)
                                                                                                aaaTCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGaa                                          <  1:3217429/148‑1 (MQ=255)
                                                                                                aaaTCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGaa                                          <  1:790978/148‑1 (MQ=255)
                                                                                                   tCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGt                                       <  1:868216/148‑1 (MQ=255)
                                                                                                    cACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTc                                      <  2:1877551/148‑1 (MQ=255)
                                                                                                    cACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTc                                      >  2:1076591/1‑148 (MQ=255)
                                                                                                    cACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTc                                      >  1:523227/1‑148 (MQ=255)
                                                                                                    cACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGAGATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTc                                      >  2:2951878/1‑148 (MQ=255)
                                                                                                      cGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAg                                    >  2:17871/1‑148 (MQ=255)
                                                                                                            ggAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCAc                              <  2:1444225/148‑1 (MQ=255)
                                                                                                             caGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCAcc                             <  2:3523396/147‑1 (MQ=255)
                                                                                                                cTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACcgcg                          <  2:2869084/148‑1 (MQ=255)
                                                                                                                  ggCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGaa                        >  2:3193046/1‑148 (MQ=255)
                                                                                                                   gcgcgcAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGaaa                       <  1:2441576/148‑1 (MQ=255)
                                                                                                                   gcgcgcAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGCGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGaaa                       >  2:359725/1‑148 (MQ=255)
                                                                                                                         aaaTTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCtt                 <  2:3289470/148‑1 (MQ=255)
                                                                                                                           aTTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCg               >  1:3357723/1‑148 (MQ=255)
                                                                                                                           aTTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCg               >  1:3411587/1‑148 (MQ=255)
                                                                                                                            ttACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCgg              >  2:3261004/1‑148 (MQ=255)
                                                                                                                              aCTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCa            <  2:1631873/148‑1 (MQ=255)
                                                                                                                               cTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCAc           >  1:2826732/1‑148 (MQ=255)
                                                                                                                                ttttACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCAct          >  1:2112188/1‑148 (MQ=255)
                                                                                                                                ttttACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCAct          >  1:182685/1‑148 (MQ=255)
                                                                                                                                    aCAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTcc      <  2:2405435/148‑1 (MQ=255)
                                                                                                                                       ggTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGt   <  2:2404190/148‑1 (MQ=255)
                                                                                                                                       ggTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGt   >  2:542568/1‑148 (MQ=255)
                                                                                                                                        gTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTc  >  2:3181786/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTC  >  CP009273/1782949‑1783232

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: