Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,897,407 A→C intergenic (+674/+699) ygcE → / ← queE FGGY‑family carbohydrate kinase/7‑carboxy‑7‑deazaguanine synthase QueE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,897,4070AC97.4% 276.9 / ‑2.5 77intergenic (+674/+699)ygcE/queEFGGY‑family carbohydrate kinase/7‑carboxy‑7‑deazaguanine synthase QueE
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (41/34);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAG  >  NZ_CP009273/2897274‑2897541
                                                                                                                                      |                                                                                                                                      
tGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGcc                                                                                                                                   >  1:591282/1‑139 (MQ=255)
  gagaACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCt                                                                                                                                  >  2:761816/1‑138 (MQ=255)
  gagaACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCt                                                                                                                                  >  2:642461/1‑138 (MQ=255)
  gagaACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCt                                                                                                                                  >  1:906516/1‑138 (MQ=255)
  gagaACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCt                                                                                                                                  >  2:1122799/1‑138 (MQ=255)
       cTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCagtag                                                                                                                            <  1:421971/139‑1 (MQ=255)
             tGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGt                                                                                                                      >  2:72567/1‑139 (MQ=255)
                  aTGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGc                                                                                                                 <  1:284040/139‑1 (MQ=255)
                    gAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCt                                                                                                                <  2:353075/138‑1 (MQ=255)
                    gAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCt                                                                                                                >  1:353075/1‑138 (MQ=255)
                    gAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTc                                                                                                               >  1:720902/1‑139 (MQ=255)
                     aGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTc                                                                                                               <  2:803809/138‑1 (MQ=255)
                         gtTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGtt                                                                                                          >  2:328541/1‑139 (MQ=255)
                             cGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATc                                                                                                      >  2:316159/1‑139 (MQ=255)
                               tGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATccc                                                                                                    >  1:1000131/1‑139 (MQ=255)
                                 gATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTc                                                                                                                     <  2:199983/120‑1 (MQ=255)
                                 gATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTc                                                                                                                     >  1:199983/1‑120 (MQ=255)
                                  atatGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGc                                                                                                 >  2:1157591/1‑139 (MQ=255)
                                    atGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTg                                                                                               >  1:393317/1‑139 (MQ=255)
                                       ttGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGgc                                                                                             <  1:316159/138‑1 (MQ=255)
                                         gCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGcgcg                                                                                          >  2:64513/1‑139 (MQ=255)
                                          cTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGcgcg                                                                                          <  1:1157591/138‑1 (MQ=255)
                                            tattaCAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCggg                                                                                        >  1:1100627/1‑138 (MQ=255)
                                               taCAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGaa                                                                                     >  1:739559/1‑138 (MQ=255)
                                                    gTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTg                                                                               >  2:522651/1‑139 (MQ=255)
                                                     tACTGTTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGt                                                                              <  2:591282/139‑1 (MQ=255)
                                                        tGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGtc                                                                                                                          >  2:772267/1‑92 (MQ=255)
                                                        tGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGtc                                                                                                                          <  1:772267/92‑1 (MQ=255)
                                                        tGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGtttt                                                                           >  2:1283639/1‑139 (MQ=255)
                                                         gCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGgcg                                                                                            >  2:437617/1‑121 (MQ=255)
                                                         gCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGgcg                                                                                            <  1:437617/121‑1 (MQ=255)
                                                         gCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTc                                                                          >  1:685250/1‑139 (MQ=255)
                                                          cTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCg                                                                         >  1:444580/1‑139 (MQ=255)
                                                          cTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCg                                                                         <  2:1063778/139‑1 (MQ=255)
                                                          cTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCg                                                                         >  1:261394/1‑139 (MQ=255)
                                                             atatAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCa                                                                       <  1:1122799/138‑1 (MQ=255)
                                                             atatAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTCTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCaa                                                                      <  2:1000131/139‑1 (MQ=255)
                                                              tataAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCaaa                                                                     >  2:1111579/1‑139 (MQ=255)
                                                               ataAAAACTTGAGAAAGAGATGACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAAt                                                                    >  2:1276739/1‑139 (MQ=255)
                                                               ataAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAAt                                                                    >  2:696533/1‑139 (MQ=255)
                                                                  aaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTa                                                                 <  1:1111579/139‑1 (MQ=255)
                                                                      cTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGa                                                             <  2:739559/139‑1 (MQ=255)
                                                                       ttGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTa                                                                 >  1:633838/1‑134 (MQ=255)
                                                                       ttGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTa                                                                 <  2:633838/134‑1 (MQ=255)
                                                                          agaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTAt                                                         <  1:642461/139‑1 (MQ=255)
                                                                                agaTAACGGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTa                                                    <  2:420731/139‑1 (MQ=255)
                                                                                 gaTAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTAtt                                                  <  2:444580/139‑1 (MQ=255)
                                                                                           ttATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTAt                                                         >  1:1155996/1‑123 (MQ=255)
                                                                                           ttATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTAt                                                         <  2:1155996/123‑1 (MQ=255)
                                                                                           ttATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATc                                         >  2:215752/1‑139 (MQ=255)
                                                                                            tataTGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGc                                                                                                 >  2:474664/1‑82 (MQ=255)
                                                                                            tataTGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGc                                                                                                 <  1:474664/82‑1 (MQ=255)
                                                                                                tggtggTTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCg                                     <  2:1046810/138‑1 (MQ=255)
                                                                                                  gtggtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTg                                  >  1:1261607/1‑139 (MQ=255)
                                                                                                     gtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGtc                                                                                                                          >  2:651977/1‑48 (MQ=255)
                                                                                                     gtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGtc                                                                                                                          <  1:651977/48‑1 (MQ=255)
                                                                                                     gtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCgggg                                                                                       >  2:495870/1‑83 (MQ=255)
                                                                                                     gtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCgggg                                                                                       <  1:495870/83‑1 (MQ=255)
                                                                                                     gtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTg                                                       >  2:240632/1‑115 (MQ=255)
                                                                                                     gtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTg                                                       <  1:240632/115‑1 (MQ=255)
                                                                                                     gtTTATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGgcg                               >  1:716850/1‑139 (MQ=255)
                                                                                                        tATCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCgg                            >  2:243943/1‑139 (MQ=255)
                                                                                                          tCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCa                                                                       >  1:781970/1‑94 (MQ=255)
                                                                                                          tCCCCGCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCa                                                                       <  2:781970/94‑1 (MQ=255)
                                                                                                               gCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTAtt                                                  <  1:196328/110‑1 (MQ=255)
                                                                                                               gCTGGCGCGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTAtt                                                  >  2:196328/1‑110 (MQ=255)
                                                                                                                   gcgcGGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGgc                 <  1:1127802/139‑1 (MQ=255)
                                                                                                                    cgcgGGGAACTCGACAGACCGGCCTCAGTAGTATCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGgcg                >  1:1108396/1‑139 (MQ=255)
                                                                                                                      cgGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCa              <  2:393317/139‑1 (MQ=255)
                                                                                                                      cgGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCa              <  1:696533/139‑1 (MQ=255)
                                                                                                                      cgGGGAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTGCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCa              <  2:229656/139‑1 (MQ=255)
                                                                                                                        gggAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGc                                                                                                                 >  1:318020/1‑38 (MQ=255)
                                                                                                                        gggAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGc                                                                                                                 <  2:318020/38‑1 (MQ=255)
                                                                                                                        gggAACTCGACAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACg            <  1:522651/139‑1 (MQ=255)
                                                                                                                                  cAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCa                                                                       <  2:1214869/70‑1 (MQ=255)
                                                                                                                                  cAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCa                                                                       >  1:1214869/1‑70 (MQ=255)
                                                                                                                                  cAGACCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAg  <  1:328541/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
TGGAGAACTATTTTGAACATGAGGTGTTACGTGGATATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAG  >  NZ_CP009273/2897274‑2897541

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: