Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,190,061 T→A 37.3% intergenic (+39/‑148) yjaG → / → hupA DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,190,0610TA37.3% 184.2 / 229.2 193intergenic (+39/‑148)yjaG/hupADUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (53/68);  new base A (30/42);  total (83/110)
Fisher's exact test for biased strand distribution p-value = 8.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 194 total aligned reads displayed.

CCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAAC  >  CP009273/4189914‑4190206
                                                                                                                                                   |                                                                                                                                                 
ccGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCt                                                                                                                                                   <  2:3173871/148‑1 (MQ=255)
     ttttAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGt                                                                                                                                              <  2:1271169/148‑1 (MQ=255)
     ttttAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGAGTTGTCGTACCTGGAGt                                                                                                                                              >  2:1922910/1‑148 (MQ=255)
       ttAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCt                                                                                                                                            >  1:647952/1‑148 (MQ=255)
         aGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTc                                                                                                                                          <  2:419063/148‑1 (MQ=255)
              aGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCttt                                                                                                                                     >  1:1072934/1‑148 (MQ=255)
                  cGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcc                                                                                                                                 <  2:1201598/148‑1 (MQ=255)
                   gaagaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGccc                                                                                                                                <  1:777366/148‑1 (MQ=255)
                    aagaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGcccc                                                                                                                               <  2:2531479/148‑1 (MQ=255)
                     agaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAg                                                                                                                                               >  2:392348/1‑131 (MQ=255)
                      gaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcccccc                                                                                                                             >  1:2137705/1‑148 (MQ=255)
                      gaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGcccccc                                                                                                                             >  1:2717638/1‑148 (MQ=255)
                        aCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGt                                                                                                                           >  2:1635864/1‑148 (MQ=255)
                          gcgATCTCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCt                                                                                                                         <  2:604522/148‑1 (MQ=255)
                           cgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTg                                                                                                                        >  1:2510144/1‑148 (MQ=255)
                           cgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTg                                                                                                                        >  2:1318774/1‑148 (MQ=255)
                            gATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTgg                                                                                                                       <  1:730119/148‑1 (MQ=255)
                             atatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGt                                                                                                                      >  1:3340695/1‑148 (MQ=255)
                              tatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTc                                                                                                                     >  1:3361954/1‑148 (MQ=255)
                                 cGAGCTGATAAAATGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTAc                                                                                                                  <  2:1559324/148‑1 (MQ=255)
                                 cGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTAc                                                                                                                  <  1:1094360/148‑1 (MQ=255)
                                    gCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAtt                                                                                                               <  1:3534075/148‑1 (MQ=255)
                                    gCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAtt                                                                                                               <  2:1985325/148‑1 (MQ=255)
                                     cTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAttt                                                                                                              >  1:900988/1‑148 (MQ=255)
                                     cTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAttt                                                                                                              >  2:1835881/1‑148 (MQ=255)
                                      tGATAACAGGCCTTAGGGCAGACCTTCGTGAGGCGGGTGAGAGCAATATTGTTATAATTTTTCAGCAATAAGACCAGAAATCGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTg                                                                                                             <  1:2647228/148‑1 (MQ=255)
                                      tGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTg                                                                                                             <  1:2647305/148‑1 (MQ=255)
                                      tGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTg                                                                                                             <  1:456149/148‑1 (MQ=255)
                                         tAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGGGAGAGCAATATTGGTATAATTTTTCAGCACTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgggg                                                                                                          >  2:2836371/1‑148 (MQ=255)
                                              ggCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGa                                                                                                     >  1:1401243/1‑148 (MQ=255)
                                               gCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                                                                                    >  1:2978101/1‑148 (MQ=255)
                                               gCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                                                                                    >  2:2858068/1‑148 (MQ=255)
                                                    aGGGCAGACCGGCGTGAGGCGGGTGAGAGCAATATTGGTATAATGTTTCAGCAATAAGACCAGAAAACGTGCTTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaaaa                                                                                               <  2:2841640/148‑1 (MQ=255)
                                                      ggCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGt                                                                                             <  2:1037777/148‑1 (MQ=255)
                                                       gCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTg                                                                                            <  1:2263016/148‑1 (MQ=255)
                                                           aCCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTa                                                                                        >  2:2135669/1‑148 (MQ=255)
                                                            ccTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt                                                                                       >  1:1584535/1‑148 (MQ=255)
                                                              tGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCg                                                                                     <  2:1413634/148‑1 (MQ=255)
                                                                   gAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCg                                                                                <  2:203306/148‑1 (MQ=255)
                                                                    aGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgt                                                                               <  2:219075/148‑1 (MQ=255)
                                                                    aGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgt                                                                               <  1:1325905/148‑1 (MQ=255)
                                                                    aGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgt                                                                               <  2:2015539/148‑1 (MQ=255)
                                                                      gCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt                                                                             >  1:3096770/1‑148 (MQ=255)
                                                                           tGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGNAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa                                                                        <  1:1805966/148‑1 (MQ=255)
                                                                           tGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa                                                                        >  2:1005400/1‑148 (MQ=255)
                                                                            gagaGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAg                                                                       <  1:3201053/148‑1 (MQ=255)
                                                                            gagaGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAg                                                                       <  1:3200736/148‑1 (MQ=255)
                                                                               agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                                    >  1:292130/1‑148 (MQ=255)
                                                                               agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                                    >  2:2496577/1‑148 (MQ=255)
                                                                                gCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgaga                                                                   <  2:1835338/148‑1 (MQ=255)
                                                                                gCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgaga                                                                   <  2:1835308/148‑1 (MQ=255)
                                                                                gCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgaga                                                                   <  1:1099875/148‑1 (MQ=255)
                                                                                 cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                  >  2:1600028/1‑148 (MQ=255)
                                                                                 cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                  >  1:916840/1‑148 (MQ=255)
                                                                                   atatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTg                                                                >  2:3110350/1‑148 (MQ=255)
                                                                                      ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTa                                                             <  2:1647393/148‑1 (MQ=255)
                                                                                      ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTa                                                             >  2:882260/1‑148 (MQ=255)
                                                                                      gtGGTATAATTTTACAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTa                                                             <  2:270462/147‑1 (MQ=255)
                                                                                        ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAtt                                                           <  1:136238/148‑1 (MQ=255)
                                                                                         gTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTc                                                          <  2:2299355/148‑1 (MQ=255)
                                                                                           ataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTg                                                        >  1:1848212/1‑148 (MQ=255)
                                                                                           ataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTg                                                        <  2:492501/148‑1 (MQ=255)
                                                                                                ttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAtt                                                   <  1:1498110/148‑1 (MQ=255)
                                                                                                ttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAtt                                                   >  1:2835074/1‑148 (MQ=255)
                                                                                                 tttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAttt                                                  <  1:3466858/148‑1 (MQ=255)
                                                                                                 attCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCGGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAttt                                                  <  1:2022215/147‑1 (MQ=255)
                                                                                                    cAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCTCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCg                                               >  1:2061051/1‑148 (MQ=255)
                                                                                                    cAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCg                                               >  1:2955744/1‑148 (MQ=255)
                                                                                                     aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGt                                              >  2:2431965/1‑148 (MQ=255)
                                                                                                     aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGt                                              >  1:3443883/1‑148 (MQ=255)
                                                                                                         aTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCa                                          <  1:871487/148‑1 (MQ=255)
                                                                                                           aaGACCAGAAAACGTGATTTAACGCCTGAGTTGTCGTAACTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGATATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGACATTTCCGTCGCACt                                        >  1:1298934/1‑148 (MQ=255)
                                                                                                             gACCAGAACACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCGTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCg                                      <  2:3439674/148‑1 (MQ=255)
                                                                                                             gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCg                                      >  2:1515935/1‑148 (MQ=255)
                                                                                                             gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCg                                      >  1:1941679/1‑148 (MQ=255)
                                                                                                                 aGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGc                                  >  1:2197771/1‑148 (MQ=255)
                                                                                                                 aGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGc                                  <  1:3149081/148‑1 (MQ=255)
                                                                                                                   aaaaCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCtt                                <  2:873401/148‑1 (MQ=255)
                                                                                                                     aaCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAg                              >  2:257820/1‑148 (MQ=255)
                                                                                                                       cGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                            <  2:2986161/148‑1 (MQ=255)
                                                                                                                         tGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAg                          >  2:2550113/1‑148 (MQ=255)
                                                                                                                           aTTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCg                        <  1:2431965/148‑1 (MQ=255)
                                                                                                                            tttAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGa                       >  2:616259/1‑148 (MQ=255)
                                                                                                                            tttAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGa                       <  1:1254147/148‑1 (MQ=255)
                                                                                                                            tttAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGa                       <  1:793634/148‑1 (MQ=255)
                                                                                                                            tttAACGCCTGATTGGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGa                       <  2:2368557/148‑1 (MQ=255)
                                                                                                                             ttAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGAt                      >  1:1180605/1‑148 (MQ=255)
                                                                                                                               aaCGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaa                    >  1:3353343/1‑148 (MQ=255)
                                                                                                                                 cGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAac                  >  1:3224506/1‑148 (MQ=255)
                                                                                                                                   ccTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac                >  1:192151/1‑148 (MQ=255)
                                                                                                                                   ccTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac                <  2:3053886/148‑1 (MQ=255)
                                                                                                                                    cTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAAcaca               <  2:1473625/148‑1 (MQ=255)
                                                                                                                                         ttGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTa          <  2:2877165/148‑1 (MQ=255)
                                                                                                                                         ttGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTa          <  2:3342005/148‑1 (MQ=255)
                                                                                                                                            tCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAgg       >  1:2017176/1‑148 (MQ=255)
                                                                                                                                             cGTACCTGGAGTCGTCCCTTTCGCCCCCCGTCGGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGTAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGa      <  1:2887611/148‑1 (MQ=255)
                                                                                                                                             cGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGa      >  2:1787201/1‑148 (MQ=255)
                                                                                                                                             cGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGa      <  1:419235/148‑1 (MQ=255)
                                                                                                                                               tACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATa    <  2:940591/148‑1 (MQ=255)
                                                                                                                                                 ccTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAAc  >  1:1133478/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAAC  >  CP009273/4189914‑4190206

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: