Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,190,061 T→A 68.4% 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,0610TA68.4% 439.8 / 286.3 313intergenic (+39/‑148)yjaG/hupADUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (39/60);  new base A (110/104);  total (149/164)
Fisher's exact test for biased strand distribution p-value = 5.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

TTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTAT  >  CP009273/4189912‑4190210
                                                                                                                                                     |                                                                                                                                                      
    gACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAg                                                                                                                                                    >  1:3886065/1‑150 (MQ=255)
     aCTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGt                                                                                                                                                   <  1:224486/150‑1 (MQ=255)
      cTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTc                                                                                                                                                  >  2:5159460/1‑150 (MQ=255)
        tttAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCtt                                                                                                                                                <  2:2324066/150‑1 (MQ=255)
          tAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTcc                                                                                                                                              <  2:2782569/150‑1 (MQ=255)
           aGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTccc                                                                                                                                             >  1:1893537/1‑150 (MQ=255)
           aGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTccc                                                                                                                                             <  1:4437272/150‑1 (MQ=255)
             cTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCtt                                                                                                                                           >  2:2648997/1‑150 (MQ=255)
               gagaGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTc                                                                                                                                         <  2:4214861/150‑1 (MQ=255)
                     gaagaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGccccc                                                                                                                                  >  2:948677/1‑150 (MQ=255)
                          aCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCt                                                                                                                             <  2:3601791/150‑1 (MQ=255)
                           cgcgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTg                                                                                                                            <  1:3791814/150‑1 (MQ=255)
                            gcgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTgg                                                                                                                           <  1:2167939/150‑1 (MQ=255)
                            gcgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTgg                                                                                                                           >  2:342529/1‑150 (MQ=255)
                             cgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGt                                                                                                                          >  2:3546720/1‑150 (MQ=255)
                              gATATCGAGCTGCTAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTc                                                                                                                         <  2:1838913/150‑1 (MQ=255)
                               atatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCt                                                                                                                        <  2:4428652/150‑1 (MQ=255)
                               atatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCt                                                                                                                        <  2:3133421/150‑1 (MQ=255)
                               atatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCt                                                                                                                        >  2:2597774/1‑150 (MQ=255)
                                tatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTa                                                                                                                       >  1:249628/1‑150 (MQ=255)
                                  tCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACa                                                                                                                     <  1:3663292/150‑1 (MQ=255)
                                   cGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAt                                                                                                                    <  1:5031187/150‑1 (MQ=255)
                                     aGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAttt                                                                                                                  <  1:4396138/150‑1 (MQ=255)
                                     aGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAttt                                                                                                                  <  1:4320130/150‑1 (MQ=255)
                                         gATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgggg                                                                                                              >  1:3518372/1‑150 (MQ=255)
                                          aTAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAAAAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTggggg                                                                                                             <  1:1426120/150‑1 (MQ=255)
                                           tAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgggggg                                                                                                            >  2:2380385/1‑150 (MQ=255)
                                            aaaaGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGc                                                                                                           <  2:5144141/150‑1 (MQ=255)
                                             aaaGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCg                                                                                                          >  1:2324307/1‑150 (MQ=255)
                                               aGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                                                                                        >  2:4994997/1‑150 (MQ=255)
                                               aGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaa                                                                                                        >  1:3400062/1‑150 (MQ=255)
                                                   cTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaaa                                                                                                    >  2:1185510/1‑150 (MQ=255)
                                                   cTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaaa                                                                                                    >  2:661832/1‑150 (MQ=255)
                                                      aGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGt                                                                                                 <  2:316520/150‑1 (MQ=255)
                                                      aGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGt                                                                                                 >  2:4840590/1‑150 (MQ=255)
                                                          cAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCt                                                                                             <  1:3193529/150‑1 (MQ=255)
                                                          cAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCt                                                                                             <  2:1377326/150‑1 (MQ=255)
                                                                 gCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTg                                                                                      <  1:1331762/150‑1 (MQ=255)
                                                                    tGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgt                                                                                   >  2:4299448/1‑150 (MQ=255)
                                                                     gAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtg                                                                                  <  2:2889984/150‑1 (MQ=255)
                                                                     gAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtg                                                                                  >  2:3699452/1‑150 (MQ=255)
                                                                      aGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt                                                                                 <  1:1945079/150‑1 (MQ=255)
                                                                      aGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt                                                                                 <  2:1328048/150‑1 (MQ=255)
                                                                         cGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATg                                                                              <  1:1806040/150‑1 (MQ=255)
                                                                         cGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATg                                                                              >  2:2110773/1‑150 (MQ=255)
                                                                           ggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa                                                                            >  2:3225036/1‑150 (MQ=255)
                                                                           ggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa                                                                            >  1:4833597/1‑150 (MQ=255)
                                                                             tGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                          >  2:2039635/1‑150 (MQ=255)
                                                                              gaggGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGga                                                                         >  1:1173579/1‑150 (MQ=255)
                                                                              gagaGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGga                                                                         >  1:4610626/1‑150 (MQ=255)
                                                                                 agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt                                                                     >  1:1074251/1‑150 (MQ=255)
                                                                                 agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                      >  1:2953565/1‑150 (MQ=255)
                                                                                 agCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                      >  1:3248343/1‑150 (MQ=255)
                                                                                  gCAATATTTGTATAATTTTTCAGCAATAAGACCAGAAAACTTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt                                                                     <  1:2700345/150‑1 (MQ=255)
                                                                                  gCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt                                                                     <  1:2076898/150‑1 (MQ=255)
                                                                                  gCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt                                                                     <  1:3732492/150‑1 (MQ=255)
                                                                                   cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTg                                                                    >  1:3304221/1‑150 (MQ=255)
                                                                                   cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTg                                                                    >  1:3152994/1‑150 (MQ=255)
                                                                                       atTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTAGCGGTGCGTGTATGCAGGAGAGTGCTa                                                                 >  2:1006536/1‑150 (MQ=255)
                                                                                        ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAtt                                                               >  1:3609341/1‑150 (MQ=255)
                                                                                         tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTc                                                              <  1:596799/150‑1 (MQ=255)
                                                                                         tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTc                                                              >  1:1961535/1‑150 (MQ=255)
                                                                                         tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTc                                                              >  2:9364/1‑150 (MQ=255)
                                                                                          ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCt                                                             <  2:692004/150‑1 (MQ=255)
                                                                                          ggTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCt                                                             >  1:1651172/1‑150 (MQ=255)
                                                                                            tataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTgg                                                           >  2:3885149/1‑150 (MQ=255)
                                                                                              taATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCa                                                         <  1:4746890/150‑1 (MQ=255)
                                                                                               aaTTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                                        >  1:457116/1‑150 (MQ=255)
                                                                                               aaTTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                                        >  1:2299762/1‑150 (MQ=255)
                                                                                                          aaTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc                                             >  2:1150609/1‑150 (MQ=255)
                                                                                                           aTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                            >  1:3223401/1‑150 (MQ=255)
                                                                                                           aTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                            >  2:5052461/1‑150 (MQ=255)
                                                                                                            tAAGACCAGAAAACGTGAATTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTc                                           <  1:1525877/150‑1 (MQ=255)
                                                                                                               gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt                                        <  2:3155386/150‑1 (MQ=255)
                                                                                                                aCCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                       <  2:1900039/150‑1 (MQ=255)
                                                                                                                  cAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCt                                     >  1:4577866/1‑150 (MQ=255)
                                                                                                                     aaaaCGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAg                                  >  2:3999888/1‑150 (MQ=255)
                                                                                                                     aaaaCGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAg                                  >  2:4181127/1‑150 (MQ=255)
                                                                                                                      aaaCGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGc                                 <  1:3253302/150‑1 (MQ=255)
                                                                                                                       aaCGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                                <  1:288044/150‑1 (MQ=255)
                                                                                                                         cGTGATTTAACGCCTTATTTGTCGTACCTGGATTCTTCCCTTTCGCCCCCCGTATGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAg                              <  1:690240/150‑1 (MQ=255)
                                                                                                                         cGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAg                              <  1:5159460/150‑1 (MQ=255)
                                                                                                                         cGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAg                              <  1:875510/150‑1 (MQ=255)
                                                                                                                         cGTGATTTAACGCCTGATTTGTCGGACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAg                              <  2:3244350/150‑1 (MQ=255)
                                                                                                                          gTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGc                             >  1:2622077/1‑150 (MQ=255)
                                                                                                                           tGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCg                            >  1:3358302/1‑150 (MQ=255)
                                                                                                                               ttAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaa                        <  2:1451789/150‑1 (MQ=255)
                                                                                                                                  aCGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAaca                     <  2:5082529/150‑1 (MQ=255)
                                                                                                                                   cGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac                    <  2:2179099/150‑1 (MQ=255)
                                                                                                                                     ccTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAt                  <  1:2787427/150‑1 (MQ=255)
                                                                                                                                     ccTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAt                  >  2:2544649/1‑150 (MQ=255)
                                                                                                                                     ccTGATTTGTCGGACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATGTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAGTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAt                  <  2:4608381/150‑1 (MQ=255)
                                                                                                                                      cTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAtt                 >  2:57107/1‑150 (MQ=255)
                                                                                                                                       tGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAATAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTg                <  2:294247/150‑1 (MQ=255)
                                                                                                                                       tGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTg                <  1:136865/150‑1 (MQ=255)
                                                                                                                                         aTTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTa              >  1:2438415/1‑150 (MQ=255)
                                                                                                                                           ttGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTTCATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAg            <  2:82098/150‑1 (MQ=255)
                                                                                                                                               cGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATa        <  1:3984888/150‑1 (MQ=255)
                                                                                                                                                gTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATaa       >  2:1563092/1‑150 (MQ=255)
                                                                                                                                                 tACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAAc      <  1:1013051/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTCCGACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTAT  >  CP009273/4189912‑4190210

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: