Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,397,366 A→T 5.9% S49S (TCT→TCA abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,397,3660AT 5.9% 587.6 / 21.3 187S49S (TCT→TCAabgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base A (90/86);  new base T (5/6);  total (95/92)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 188 total aligned reads displayed.

GTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATC  >  CP009273/1397220‑1397513
                                                                                                                                                  |                                                                                                                                                   
gTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTGGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAg                                                                                                                                                    <  2:938117/148‑1 (MQ=255)
  aGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGAt                                                                                                                                                  >  1:1581767/1‑148 (MQ=255)
   ggTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGAtt                                                                                                                                                 >  2:159909/1‑148 (MQ=255)
      tGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCa                                                                                                                                              <  2:918995/148‑1 (MQ=255)
          cgcAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAgcgc                                                                                                                                          >  2:702389/1‑148 (MQ=255)
          cgcAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAgcgc                                                                                                                                          >  1:2271471/1‑148 (MQ=255)
            cAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCg                                                                                                                                        >  2:863018/1‑148 (MQ=255)
               ccTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAg                                                                                                                                     >  2:1899751/1‑148 (MQ=255)
               ccTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAg                                                                                                                                     >  2:3108826/1‑148 (MQ=255)
               ccTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGAGCCGAAg                                                                                                                                     >  1:2728578/1‑148 (MQ=255)
                cTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGc                                                                                                                                    <  1:159935/148‑1 (MQ=255)
                cTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGc                                                                                                                                    >  2:2792175/1‑148 (MQ=255)
                 tgcttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGcc                                                                                                                                   >  1:2812193/1‑148 (MQ=255)
                  gcttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCa                                                                                                                                  >  1:1676561/1‑148 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                                 >  1:2253796/1‑148 (MQ=255)
                     tgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGAt                                                                                                                               <  1:377735/148‑1 (MQ=255)
                     tgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGAt                                                                                                                               >  2:2488419/1‑148 (MQ=255)
                      gctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATg                                                                                                                              <  1:2037155/148‑1 (MQ=255)
                         gACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTc                                                                                                                           <  2:2977237/148‑1 (MQ=255)
                           cTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCg                                                                                                                         <  1:331393/148‑1 (MQ=255)
                             tAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCt                                                                                                                       <  2:3178787/148‑1 (MQ=255)
                               aaCCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGa                                                                                                                     >  2:2804945/1‑148 (MQ=255)
                               aaCCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGa                                                                                                                     >  2:2805420/1‑148 (MQ=255)
                                 ccTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGAcc                                                                                                                   <  1:2015717/148‑1 (MQ=255)
                                 ccTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGAcc                                                                                                                   >  2:2903137/1‑148 (MQ=255)
                                     ccAGGGCGTCATATTCTCCCAGCAGGGCGAGAACCGGTGTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGaa                                                                                                               <  2:1567496/148‑1 (MQ=255)
                                      cAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAAc                                                                                                              >  1:1410725/1‑148 (MQ=255)
                                       aGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAAct                                                                                                             <  2:1501165/148‑1 (MQ=255)
                                        gggCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAActc                                                                                                            <  2:1530278/148‑1 (MQ=255)
                                        gggCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAActc                                                                                                            >  1:826163/1‑148 (MQ=255)
                                         ggCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAActct                                                                                                           >  1:543458/1‑148 (MQ=255)
                                                  ttCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACgt                                                                                                  <  2:2689389/148‑1 (MQ=255)
                                                      cccAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGttt                                                                                              >  2:1984229/1‑148 (MQ=255)
                                                      cccAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGttt                                                                                              <  2:591601/148‑1 (MQ=255)
                                                      cccAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGttt                                                                                              <  2:2253796/148‑1 (MQ=255)
                                                        cagcagGGCGATAACCGGTGTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCt                                                                                            <  1:2051638/148‑1 (MQ=255)
                                                           cagGGCGATAACCGGTTGGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                         <  2:3050718/148‑1 (MQ=255)
                                                              ggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatga                                                                                      <  1:2498901/148‑1 (MQ=255)
                                                              ggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatga                                                                                      <  1:2293183/148‑1 (MQ=255)
                                                              ggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatga                                                                                      >  1:2012581/1‑148 (MQ=255)
                                                               gCGATTACCGGTTTGCCTTGACCAAACGATGCAATAAAGGCTTTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatgat                                                                                     <  1:1071143/148‑1 (MQ=255)
                                                                  aTAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATccc                                                                                  >  1:2767290/1‑148 (MQ=255)
                                                                    aaCCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAg                                                                                <  1:2558643/148‑1 (MQ=255)
                                                                     aCCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGa                                                                               <  1:127740/148‑1 (MQ=255)
                                                                     aCCGGTGTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGa                                                                               <  1:129431/148‑1 (MQ=255)
                                                                       cGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGAtt                                                                             >  2:2872531/1‑148 (MQ=255)
                                                                       cGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGAtt                                                                             <  1:1899751/148‑1 (MQ=255)
                                                                        ggTTTGCCTTGAACAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCATCGCCGAAGCCAGATGCTCAGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGAttt                                                                            <  2:3285089/148‑1 (MQ=255)
                                                                          tttGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGa                                                                          <  2:2150978/148‑1 (MQ=255)
                                                                           ttGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGAt                                                                         >  2:433201/1‑148 (MQ=255)
                                                                              ccTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCgg                                                                      >  1:2155994/1‑148 (MQ=255)
                                                                              ccTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCgg                                                                      <  1:1704412/148‑1 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    <  1:864931/148‑1 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    >  1:3057822/1‑148 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    <  1:2922618/148‑1 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    >  2:1312418/1‑148 (MQ=255)
                                                                                   aCCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAata                                                                 >  1:540638/1‑148 (MQ=255)
                                                                                   aCCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAata                                                                 <  2:1808708/148‑1 (MQ=255)
                                                                                   aCCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAata                                                                 <  1:2113816/148‑1 (MQ=255)
                                                                                    ccAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAatat                                                                >  2:3136063/1‑148 (MQ=255)
                                                                                    ccAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAatat                                                                <  1:2921066/148‑1 (MQ=255)
                                                                                     cAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATc                                                               <  1:2629770/148‑1 (MQ=255)
                                                                                        actaaGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGt                                                            <  2:2516509/145‑1 (MQ=255)
                                                                                        aCGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGt                                                            >  1:2029307/1‑148 (MQ=255)
                                                                                           aaGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGtata                                                         <  1:3108826/148‑1 (MQ=255)
                                                                                            aGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATaa                                                        >  1:487262/1‑148 (MQ=255)
                                                                                             gCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAAc                                                       >  1:643343/1‑148 (MQ=255)
                                                                                              cAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACg                                                      <  1:415028/148‑1 (MQ=255)
                                                                                               aaTAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGt                                                     <  1:2214373/148‑1 (MQ=255)
                                                                                                  aaaGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTgg                                                  <  1:3028995/148‑1 (MQ=255)
                                                                                                    aGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCg                                                >  1:2591397/1‑148 (MQ=255)
                                                                                                     ggCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGa                                               <  2:1459677/148‑1 (MQ=255)
                                                                                                      gCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGAt                                              <  2:1676561/148‑1 (MQ=255)
                                                                                                        aTTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCg                                            >  2:2362543/1‑148 (MQ=255)
                                                                                                         tttGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCgg                                           >  2:1251139/1‑148 (MQ=255)
                                                                                                          ttGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGc                                          <  2:1473602/148‑1 (MQ=255)
                                                                                                            gggATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCtt                                        >  2:1742944/1‑148 (MQ=255)
                                                                                                            gggATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCtt                                        >  2:359794/1‑148 (MQ=255)
                                                                                                             ggATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTc                                       >  1:2389334/1‑148 (MQ=255)
                                                                                                             ggATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTc                                       >  1:1770385/1‑148 (MQ=255)
                                                                                                              gATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCa                                      >  1:2641118/1‑148 (MQ=255)
                                                                                                              gATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCa                                      >  1:988724/1‑148 (MQ=255)
                                                                                                                tatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAAt                                    <  2:2950808/148‑1 (MQ=255)
                                                                                                                 atTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATc                                   >  1:313170/1‑148 (MQ=255)
                                                                                                                    gCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCa                                <  1:1852139/148‑1 (MQ=255)
                                                                                                                     ccTACGTTGTGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAt                               <  1:101650/148‑1 (MQ=255)
                                                                                                                          gTTGCGTGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAAATCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgtcg                          <  1:2840831/148‑1 (MQ=255)
                                                                                                                          gTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgtcg                          <  2:2282856/148‑1 (MQ=255)
                                                                                                                          gTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgtcg                          <  1:2409526/148‑1 (MQ=255)
                                                                                                                           ttGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGa                         <  1:963462/148‑1 (MQ=255)
                                                                                                                            tGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGAt                        >  2:283134/1‑148 (MQ=255)
                                                                                                                             gCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATa                       >  1:2348040/1‑148 (MQ=255)
                                                                                                                             gCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATa                       <  1:2436962/148‑1 (MQ=255)
                                                                                                                               gggTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaa                     <  2:1182548/148‑1 (MQ=255)
                                                                                                                               gggTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaa                     <  2:530636/148‑1 (MQ=255)
                                                                                                                                    aCGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGAt                <  2:1826615/148‑1 (MQ=255)
                                                                                                                                           aGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTc         <  2:2637598/148‑1 (MQ=255)
                                                                                                                                             ccTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCt       >  2:1232087/1‑148 (MQ=255)
                                                                                                                                               tGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGc     >  2:1950435/1‑148 (MQ=255)
                                                                                                                                                  aGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATc  <  1:1035401/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATC  >  CP009273/1397220‑1397513

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: