Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,190,061 T→A 36.2% 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,0610TA36.2% 207.8 / 220.1 193intergenic (+39/‑148)yjaG/hupADUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (57/66);  new base A (27/43);  total (84/109)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

CCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACT  >  CP009273/4189914‑4190207
                                                                                                                                                   |                                                                                                                                                  
 cGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTg                                                                                                                                                   <  2:1790135/148‑1 (MQ=255)
   aCTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGa                                                                                                                                                 >  2:618326/1‑148 (MQ=255)
   aCTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGAGTTGTCGTACCAGGa                                                                                                                                                 >  1:597232/1‑148 (MQ=255)
    cTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAg                                                                                                                                                >  2:2633253/1‑148 (MQ=255)
    cTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAg                                                                                                                                                >  1:1841611/1‑148 (MQ=255)
         aGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTc                                                                                                                                           <  1:2385776/148‑1 (MQ=255)
          gCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTcc                                                                                                                                          >  2:985024/1‑148 (MQ=255)
          gCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTcc                                                                                                                                          >  2:1843484/1‑148 (MQ=255)
           cTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTccc                                                                                                                                         <  1:2217824/148‑1 (MQ=255)
            tGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGACTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCt                                                                                                                                        <  2:584231/148‑1 (MQ=255)
             gAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCtt                                                                                                                                       >  1:702339/1‑148 (MQ=255)
              aGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCttt                                                                                                                                      >  2:1225994/1‑148 (MQ=255)
                 gCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGc                                                                                                                                   >  1:1085648/1‑148 (MQ=255)
                    aagaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcccc                                                                                                                                <  2:2318540/148‑1 (MQ=255)
                     agaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGccccc                                                                                                                               >  2:1285357/1‑148 (MQ=255)
                      gaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGcccccc                                                                                                                              >  2:515469/1‑148 (MQ=255)
                       aaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCg                                                                                                                             >  2:45505/1‑148 (MQ=255)
                        aCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGt                                                                                                                            <  2:2264272/148‑1 (MQ=255)
                        aCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGt                                                                                                                            >  1:1706525/1‑148 (MQ=255)
                           cgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTg                                                                                                                         >  1:442499/1‑148 (MQ=255)
                            gATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTgg                                                                                                                        <  2:316596/148‑1 (MQ=255)
                             atatCGAGCTGATAAAGGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGt                                                                                                                       >  2:1835321/1‑148 (MQ=255)
                                tCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTa                                                                                                                    >  1:1457041/1‑148 (MQ=255)
                                 cGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTAc                                                                                                                   >  2:2545625/1‑148 (MQ=255)
                                   aGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAt                                                                                                                 <  1:2330360/148‑1 (MQ=255)
                                   aGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAt                                                                                                                 >  1:3207498/1‑148 (MQ=255)
                                    gCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAtt                                                                                                                <  2:1576785/148‑1 (MQ=255)
                                     cTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAttt                                                                                                               <  2:1711416/148‑1 (MQ=255)
                                         tAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgggg                                                                                                           >  1:360751/1‑148 (MQ=255)
                                         tAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgggg                                                                                                           >  1:2980590/1‑148 (MQ=255)
                                             aGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCg                                                                                                       >  1:1304744/1‑148 (MQ=255)
                                             aGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCg                                                                                                       >  2:1426787/1‑148 (MQ=255)
                                               gCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                                                                                     >  2:468999/1‑148 (MQ=255)
                                               gCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACAAGGAGTCTTCCATTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                                                                                     >  2:1008375/1‑148 (MQ=255)
                                                ccTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaa                                                                                                    >  1:1724849/1‑148 (MQ=255)
                                                ccTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATAGTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaa                                                                                                    <  2:3197270/148‑1 (MQ=255)
                                                  ttAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaa                                                                                                  <  1:826714/148‑1 (MQ=255)
                                                      ggCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGt                                                                                              <  2:2447110/148‑1 (MQ=255)
                                                        cAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTCTAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTgg                                                                                            <  2:730823/148‑1 (MQ=255)
                                                        cAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTgg                                                                                            <  2:791224/148‑1 (MQ=255)
                                                        cAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCGTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTgg                                                                                            <  2:177143/148‑1 (MQ=255)
                                                        cAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTgg                                                                                            <  1:774678/148‑1 (MQ=255)
                                                         aGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGc                                                                                           <  2:685905/148‑1 (MQ=255)
                                                           aCCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTa                                                                                         >  1:3089464/1‑148 (MQ=255)
                                                            ccTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAt                                                                                        >  2:855396/1‑148 (MQ=255)
                                                               gCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCgg                                                                                     >  2:2072924/1‑148 (MQ=255)
                                                               gCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCgg                                                                                     >  1:1816075/1‑148 (MQ=255)
                                                                   gAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCg                                                                                 >  1:2794629/1‑148 (MQ=255)
                                                                      gCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt                                                                              <  2:3004537/148‑1 (MQ=255)
                                                                       cGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTa                                                                             <  1:2545625/148‑1 (MQ=255)
                                                                         ggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATg                                                                           <  1:1166857/148‑1 (MQ=255)
                                                                           tGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa                                                                         <  2:2817186/148‑1 (MQ=255)
                                                                           tGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa                                                                         <  1:515469/148‑1 (MQ=255)
                                                                             agagCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                       <  2:56518/148‑1 (MQ=255)
                                                                               agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGTAGgag                                                                     <  2:2642152/148‑1 (MQ=255)
                                                                               agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                                     <  2:786976/148‑1 (MQ=255)
                                                                               agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGgag                                                                     <  2:1103850/148‑1 (MQ=255)
                                                                                 cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                   <  2:758052/148‑1 (MQ=255)
                                                                                 cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                   <  2:810673/148‑1 (MQ=255)
                                                                                   atatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGACTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGGCTGGTCTACATTTGGGGGGCGTAAAAAAGTGGCTATCGGTGCGTGTATGAAGGAGAGTg                                                                 >  2:1589791/1‑148 (MQ=255)
                                                                                     atTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCt                                                               >  1:305871/1‑148 (MQ=255)
                                                                                      ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTa                                                              <  1:2867379/148‑1 (MQ=255)
                                                                                      ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTa                                                              >  2:3224394/1‑148 (MQ=255)
                                                                                        ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTTACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAtt                                                            <  1:3112096/148‑1 (MQ=255)
                                                                                        ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAtt                                                            <  1:1486823/148‑1 (MQ=255)
                                                                                            taATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTgg                                                        >  1:3265486/1‑148 (MQ=255)
                                                                                            taATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTgg                                                        >  2:368894/1‑148 (MQ=255)
                                                                                            taATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGCGGGGCGAAAAAAAGTGGCTCTCGGTGCGAGTATGCAGGAGAGTGCTATTCTgg                                                        >  2:2003951/1‑148 (MQ=255)
                                                                                                ttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCGATTCTGGCAtt                                                    >  2:2690913/1‑148 (MQ=255)
                                                                                                 tttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAttt                                                   >  1:1243568/1‑148 (MQ=255)
                                                                                                  ttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTc                                                  <  2:943721/148‑1 (MQ=255)
                                                                                                    cAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCg                                                <  1:1619107/148‑1 (MQ=255)
                                                                                                    cAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCg                                                >  2:1524575/1‑148 (MQ=255)
                                                                                                     aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGt                                               <  2:328719/148‑1 (MQ=255)
                                                                                                     aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGt                                               <  2:2931341/148‑1 (MQ=255)
                                                                                                         aTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCa                                           >  1:2523866/1‑148 (MQ=255)
                                                                                                            aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACGTTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTc                                        >  1:1950571/1‑148 (MQ=255)
                                                                                                            aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTc                                        >  1:1823132/1‑148 (MQ=255)
                                                                                                             gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCg                                       >  1:187248/1‑148 (MQ=255)
                                                                                                                 aGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGc                                   <  1:545615/148‑1 (MQ=255)
                                                                                                                       cGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                             <  1:618326/148‑1 (MQ=255)
                                                                                                                       cGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                             <  2:2899517/148‑1 (MQ=255)
                                                                                                                         tGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAg                           <  1:2433027/148‑1 (MQ=255)
                                                                                                                          gATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGc                          >  2:2361012/1‑148 (MQ=255)
                                                                                                                           aTTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCg                         <  2:2776269/148‑1 (MQ=255)
                                                                                                                             ttAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTAGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGAt                       >  1:2705456/1‑148 (MQ=255)
                                                                                                                                    cTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAAcaca                <  1:985024/148‑1 (MQ=255)
                                                                                                                                    cTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAAcaca                <  2:1003362/148‑1 (MQ=255)
                                                                                                                                    cTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAAcaca                <  2:3244078/148‑1 (MQ=255)
                                                                                                                                     tGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAt               <  2:3079709/148‑1 (MQ=255)
                                                                                                                                      gATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAtt              <  2:354517/148‑1 (MQ=255)
                                                                                                                                      gATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAtt              >  2:2208418/1‑148 (MQ=255)
                                                                                                                                       aTTTGTCGTACCTGGAGTCGTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTg             <  2:2413334/148‑1 (MQ=255)
                                                                                                                                         ttGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTa           >  2:2763038/1‑148 (MQ=255)
                                                                                                                                          tGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTaa          <  2:258993/148‑1 (MQ=255)
                                                                                                                                            tCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAgg        <  2:632055/148‑1 (MQ=255)
                                                                                                                                            tCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAgg        <  1:503750/148‑1 (MQ=255)
                                                                                                                                             cGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGa       <  1:32757/148‑1 (MQ=255)
                                                                                                                                              gTACCATGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTTGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt      <  1:1008375/148‑1 (MQ=255)
                                                                                                                                                aCCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATaa    >  1:2012573/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACT  >  CP009273/4189914‑4190207

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: