Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,190,061 T→A 73.7% 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,0610TA73.7% 383.3 / 162.5 209intergenic (+39/‑148)yjaG/hupADUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (23/32);  new base A (70/84);  total (93/116)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 209 total aligned reads displayed.

GACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAAC  >  CP009273/4189916‑4190206
                                                                                                                                                 |                                                                                                                                                 
gACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTgg                                                                                                                                                 >  2:3000200/1‑148 (MQ=255)
 tcTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATATAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAAGATTGGTATTATTTTTCAGCAATAAGACCAGAAATCGTGATTTATCGCCTGATTTTTGGTACCAGGa                                                                                                                                                <  2:376816/147‑1 (MQ=255)
 aCTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGa                                                                                                                                                <  2:3173916/148‑1 (MQ=255)
  cTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAg                                                                                                                                               <  1:2783548/148‑1 (MQ=255)
  cTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAg                                                                                                                                               >  2:3259400/1‑148 (MQ=255)
  cTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAg                                                                                                                                               >  1:852950/1‑148 (MQ=255)
      tAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCtt                                                                                                                                           <  1:2530338/148‑1 (MQ=255)
       aGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTc                                                                                                                                          >  2:1036608/1‑148 (MQ=255)
        gCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTcc                                                                                                                                         <  1:1328971/148‑1 (MQ=255)
           gAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCtt                                                                                                                                      >  2:1234265/1‑148 (MQ=255)
              tGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCg                                                                                                                                   >  1:2658840/1‑148 (MQ=255)
              tGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCg                                                                                                                                   >  2:2539422/1‑148 (MQ=255)
              tGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCg                                                                                                                                   >  2:1133347/1‑148 (MQ=255)
                  aagaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGcccc                                                                                                                               <  1:1434822/148‑1 (MQ=255)
                  aagaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGcccc                                                                                                                               <  1:1433141/148‑1 (MQ=255)
                    gaaCGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGcccccc                                                                                                                             >  2:1504923/1‑148 (MQ=255)
                       cgcgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTc                                                                                                                          <  1:1470278/148‑1 (MQ=255)
                         cgATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTg                                                                                                                        >  2:1463348/1‑148 (MQ=255)
                          gATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTgg                                                                                                                       <  1:3362399/148‑1 (MQ=255)
                            tatCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTc                                                                                                                     >  2:2790616/1‑148 (MQ=255)
                                  gCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAtt                                                                                                               <  2:1399489/148‑1 (MQ=255)
                                  gCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAtt                                                                                                               <  2:873171/148‑1 (MQ=255)
                                   cTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGGCTGGTCTACAttt                                                                                                              >  1:2400005/1‑148 (MQ=255)
                                    tGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTg                                                                                                             <  2:1759210/148‑1 (MQ=255)
                                    tGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTg                                                                                                             <  2:3521131/148‑1 (MQ=255)
                                     gATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgg                                                                                                            <  2:2381558/148‑1 (MQ=255)
                                      aTAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTggg                                                                                                           <  2:231029/148‑1 (MQ=255)
                                        aaaaGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTggggg                                                                                                         >  2:1619372/1‑148 (MQ=255)
                                         aaaGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgggggg                                                                                                        >  2:485830/1‑148 (MQ=255)
                                         aaaGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTgggggg                                                                                                        <  2:3113726/148‑1 (MQ=255)
                                           aGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCg                                                                                                      >  2:406507/1‑148 (MQ=255)
                                           aGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCg                                                                                                      >  1:764165/1‑148 (MQ=255)
                                              ccTTCGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaa                                                                                                   <  2:3144712/148‑1 (MQ=255)
                                                  aGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaaaa                                                                                               <  1:3046899/148‑1 (MQ=255)
                                                   gggCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAg                                                                                              >  1:1679495/1‑148 (MQ=255)
                                                    ggCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGt                                                                                             <  1:3150140/148‑1 (MQ=255)
                                                      cAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTgg                                                                                           <  1:2506295/148‑1 (MQ=255)
                                                       aGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGc                                                                                          <  1:1406990/148‑1 (MQ=255)
                                                           cTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATc                                                                                      <  1:1341386/148‑1 (MQ=255)
                                                               gTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTg                                                                                  <  2:1962679/148‑1 (MQ=255)
                                                                    gCGGGTGAGAGCAATATTGGTATAATTTTTCAGAAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt                                                                             <  1:3255900/148‑1 (MQ=255)
                                                                     cGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTa                                                                            <  2:3364421/148‑1 (MQ=255)
                                                                         tGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCa                                                                        >  2:2586312/1‑148 (MQ=255)
                                                                          gagaGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAg                                                                       >  2:764358/1‑148 (MQ=255)
                                                                            gagCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGga                                                                     >  1:1122463/1‑148 (MQ=255)
                                                                                aaTATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGt                                                                 <  1:3291690/148‑1 (MQ=255)
                                                                                 atatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTg                                                                <  1:2310577/148‑1 (MQ=255)
                                                                                  tatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGc                                                               >  1:3208799/1‑148 (MQ=255)
                                                                                  tatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGc                                                               >  1:2438094/1‑148 (MQ=255)
                                                                                  tatTGGTATAATTTTTCAGAAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGc                                                               >  2:1541555/1‑148 (MQ=255)
                                                                                    ttGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTa                                                             >  1:595666/1‑148 (MQ=255)
                                                                                     tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAt                                                            >  1:356811/1‑148 (MQ=255)
                                                                                        tataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCt                                                         >  1:1813307/1‑148 (MQ=255)
                                                                                        tataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCt                                                         >  1:3159678/1‑148 (MQ=255)
                                                                                         ataATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGCTTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCGGGGCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAATCTg                                                        >  1:1878854/1‑148 (MQ=255)
                                                                                            aTTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCa                                                     <  2:1799759/148‑1 (MQ=255)
                                                                                             tttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                                    <  2:542704/148‑1 (MQ=255)
                                                                                             tttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                                    >  1:3460488/1‑148 (MQ=255)
                                                                                             tttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAt                                                    >  1:2937960/1‑148 (MQ=255)
                                                                                              ttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAtt                                                   >  1:3090760/1‑148 (MQ=255)
                                                                                               tttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAttt                                                  <  2:1692318/148‑1 (MQ=255)
                                                                                                ttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTc                                                 <  1:2649998/148‑1 (MQ=255)
                                                                                                   aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGt                                              <  2:2147902/148‑1 (MQ=255)
                                                                                                       aTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTCTTCTGGCATTTCCGTCGCa                                          >  2:2069610/1‑148 (MQ=255)
                                                                                                        tAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc                                         >  1:1635493/1‑148 (MQ=255)
                                                                                                        tAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc                                         <  2:788561/148‑1 (MQ=255)
                                                                                                        tAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCAc                                         <  1:255574/148‑1 (MQ=255)
                                                                                                         aaGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                        <  1:1959091/148‑1 (MQ=255)
                                                                                                         aaGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                        <  1:3448919/148‑1 (MQ=255)
                                                                                                         aaGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                        >  2:2507968/1‑148 (MQ=255)
                                                                                                          aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTc                                       <  2:2412608/148‑1 (MQ=255)
                                                                                                           gACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCg                                      >  1:1676490/1‑148 (MQ=255)
                                                                                                            aCCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGa                                     >  2:2818444/1‑148 (MQ=255)
                                                                                                             ccAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCCTTTCCGTCGCACTCGAt                                    >  2:1199782/1‑148 (MQ=255)
                                                                                                             ccAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt                                    <  1:22848/148‑1 (MQ=255)
                                                                                                             ccAGAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt                                    <  1:2719184/148‑1 (MQ=255)
                                                                                                               aGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGc                                  <  1:2346333/148‑1 (MQ=255)
                                                                                                               aGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGc                                  <  2:1122463/148‑1 (MQ=255)
                                                                                                                gAAAACGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCt                                 >  2:307948/1‑148 (MQ=255)
                                                                                                                 aaaaCGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCGTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCtt                                <  2:2355484/148‑1 (MQ=255)
                                                                                                                     cGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                            <  2:1772567/148‑1 (MQ=255)
                                                                                                                     cGTGATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCa                            <  2:706652/148‑1 (MQ=255)
                                                                                                                       tGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAg                          <  1:477564/148‑1 (MQ=255)
                                                                                                                        gATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGc                         >  2:672081/1‑148 (MQ=255)
                                                                                                                        gATTTAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGAGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGc                         >  1:1528734/1‑148 (MQ=255)
                                                                                                                          tttAACGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGa                       >  1:1863924/1‑148 (MQ=255)
                                                                                                                           ttAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGAt                      >  1:2207483/1‑148 (MQ=255)
                                                                                                                             aaCGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaa                    >  2:2278555/1‑148 (MQ=255)
                                                                                                                               cGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAac                  >  2:719997/1‑148 (MQ=255)
                                                                                                                               cGCCTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAac                  <  2:236167/148‑1 (MQ=255)
                                                                                                                                 ccTGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac                >  1:2060464/1‑148 (MQ=255)
                                                                                                                                   tGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCACGCGATAAACACAt              >  2:1921914/1‑148 (MQ=255)
                                                                                                                                   tGATTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACAGTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGATTAGCAAGCGATAAACACAt              >  1:1040764/1‑148 (MQ=255)
                                                                                                                                     aTTTGTCGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTg            <  2:798255/148‑1 (MQ=255)
                                                                                                                                           cGTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGa      <  1:2880542/148‑1 (MQ=255)
                                                                                                                                            gTACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt     <  2:2937960/148‑1 (MQ=255)
                                                                                                                                             tACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATa    <  1:1056000/148‑1 (MQ=255)
                                                                                                                                             tACCAGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATa    <  1:609278/148‑1 (MQ=255)
                                                                                                                                              aCCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCgtgt                                                                             >  1:1488722/1‑74 (MQ=255)
                                                                                                                                               ccTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAAc  >  1:526001/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GACTTTTAGCTGAGTGCGAAGAACGCGATATCGAGCTGATAAAAGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAAC  >  CP009273/4189916‑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: