Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,017,350 G→A 5.6% W18* (TGG→TGA) ‡ U712_09990 → Prophage‑derived protein uncharacterized protein yozM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,017,3500GA 5.6% 224.9 / 3.0 72W18* (TGG→TGA) ‡U712_09990Prophage‑derived protein uncharacterized protein yozM
Reads supporting (aligned to +/- strand):  ref base G (40/27);  new base A (2/2);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

ATAAAATATAAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAATAGGCATGGATCAAACCAA  >  CP006881/2017216‑2017484
                                                                                                                                      |                                                                                                                                      
ataaAATATAAACGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTTTGATTGTGAATGAAGAAGAGGTTAATTTGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGta                                                                                                                                    <  5:18889/139‑1 (MQ=255)
 taaAATATAAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGtat                                                                                                                                   <  4:253464/139‑1 (MQ=255)
 taaAATATAAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGtat                                                                                                                                   <  8:197404/139‑1 (MQ=255)
 taaAATATAAAGGTACATATTGATATCCGTAGTATTAATTAAGGAGTGATGCTTTTCGAGGTGGGAGTGTGTGATGGTTAATGAAGAAGAGGCTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGta                                                                                                                                    <  3:365739/138‑1 (MQ=255)
  aaaaTATAAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGtat                                                                                                                                   >  3:289409/1‑138 (MQ=255)
    aaTATAAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGtat                                                                                                                                   >  7:374678/1‑136 (MQ=255)
    aaTATAAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGtat                                                                                                                                   <  8:374678/136‑1 (MQ=255)
       ataAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAAGtt                                                                                                                              >  7:222952/1‑138 (MQ=255)
                aataatGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAAt                                                                                                                     >  2:154702/1‑138 (MQ=255)
                aataatGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAAt                                                                                                                     >  5:205456/1‑138 (MQ=255)
                aataatGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTGTAATCGAAt                                                                                                                     >  5:335513/1‑138 (MQ=255)
                 ataatGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATca                                                                                                                   <  3:15391/139‑1 (MQ=255)
                    atGATATCCGTAGTATTAATAAAGGAGCGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATcaaac                                                                                                                <  3:116059/139‑1 (MQ=255)
                     tGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATcaaac                                                                                                                >  2:320652/1‑138 (MQ=255)
                       atatCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATcaaacaa                                                                                                              >  3:399333/1‑138 (MQ=255)
                                gtaTTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGta                                                                                                                                    >  6:225557/1‑107 (MQ=255)
                                gtaTTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGta                                                                                                                                    <  5:225557/107‑1 (MQ=255)
                                   taaaaaaaGGAGACTTTCTTTTCGACGCGGTCTAGTGTGATGGGCATTGAAGAAGAGGTTAATTGGATTTTTGGTCTTTGTTCCTGCTTTGATTTTGTGGGGTATAACTTTAATCGAATCAAACAAAACGACCCCTGTa                                                                                                 <  2:123661/134‑1 (MQ=255)
                                     aataaAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAg                                                                                                <  1:154702/138‑1 (MQ=255)
                                      ataaAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCATACAAAAAGACCCCTGTAg                                                                                                >  4:180509/1‑137 (MQ=255)
                                      ataaAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAg                                                                                                >  1:324522/1‑137 (MQ=255)
                                           ggAGGGATTCTTTTCCATGAAGGAGAGTGTAATGGTTAATGAAGAAGAGGCTAATCGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTCAGGta                                                                                                                                    <  3:439780/96‑1 (MQ=255)
                                           ggAGGGATTCTTTTCCATGAAGGAGAGTGTAATGGTTAATGAAGAAGAGGCTAATCGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTCAGGta                                                                                                                                    >  4:439780/1‑96 (MQ=255)
                                           ggAGGGATTCTTTTCCATGAAGGAGAGTGTAATGGTTAATGAAGAAGAGGCTAATCGCATTTTTGGTCTTAGTTCCTGCTTTGATTATGTCAGGTATTACTTTAATCGAAGCAAATAAAAAGTCCCCTTTGGAAGtaat                                                                                         >  8:348606/1‑136 (MQ=255)
                                           ggAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGttt                                                                                          >  3:258280/1‑138 (MQ=255)
                                           ggAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGttt                                                                                          >  6:331091/1‑138 (MQ=255)
                                           ggAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGttt                                                                                          >  6:382496/1‑138 (MQ=255)
                                           ggAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACACCTGAAGAAGtttt                                                                                         >  5:27753/1‑139 (MQ=255)
                                              gagaTTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGa                                                                                      <  5:277379/139‑1 (MQ=255)
                                               agaTTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGa                                                                                      >  8:17651/1‑138 (MQ=255)
                                                 aTTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAg                                                                                   >  7:99413/1‑139 (MQ=255)
                                                 aTTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAg                                                                                   >  2:361362/1‑139 (MQ=255)
                                                      tttCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCtt                                                                              >  8:164081/1‑139 (MQ=255)
                                                      tttCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCtt                                                                              >  3:177927/1‑139 (MQ=255)
                                                      tttCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCtt                                                                              >  7:423533/1‑139 (MQ=255)
                                                      ttgagccgaGGGAGAGTGTTCTTGTGAATGAAGAATAGGTTAATTGGATTTTTGGTCTTTGTTCCTGCTTTTTTTATGTGGGGTATAACTTTAATCTAATAAAAAAAAACGACCCCTGTAGAAGTTTTAGAAAGTGCt                                                                               <  4:51864/132‑1 (MQ=255)
                                                          gacgaGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAACGACCCCTGTAGAAGTTTTAGAAAGTGCTTggg                                                                           <  7:403477/138‑1 (MQ=255)
                                                          gacgaGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTggg                                                                           <  5:382496/138‑1 (MQ=255)
                                                            cgaGGGATAGTGTGATGGTGAATGAAGAAGTGGTTAATTGTATTTTTTGTCTTAGTTCCTGCTTTGATTATGTGGGTTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGa                                                                          <  6:205456/137‑1 (MQ=255)
                                                                ggattgtgtGATTGTGAATTAAGAAGCGTTTTTTTGTATTTTTGGTCTTAGTTGTTTCTTTGTTTATGTGGGGTATAATTTTAATCGAATAAAAAAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAAt                                                                     <  6:27753/133‑1 (MQ=255)
                                                                 gagaGTTTGATGGTGAATGAATACTCGGTTAATTGGATTTTTGGTCTTAGTGCCTGCTTTTATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAAttt                                                                   <  5:82427/139‑1 (MQ=255)
                                                                          aTGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTAttt                                                          >  5:70361/1‑139 (MQ=255)
                                                                          aTGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCAGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTAtt                                                           >  2:116813/1‑138 (MQ=255)
                                                                              tgaatgaaGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGtt                                                      >  5:211414/1‑139 (MQ=255)
                                                                               gaatgaaGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGttt                                                     >  1:378036/1‑139 (MQ=255)
                                                                                aatgaaGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGtttt                                                    >  5:280518/1‑139 (MQ=255)
                                                                                aatgaaGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGtttt                                                    >  6:25438/1‑139 (MQ=255)
                                                                                  tgaaGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGTTGAATTTGGGTTATTTAGTTTTg                                                   >  1:328275/1‑138 (MQ=255)
                                                                                    aagaagAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGaaa                                                >  5:131795/1‑139 (MQ=255)
                                                                                    aagaagAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGaa                                                 <  2:378036/138‑1 (MQ=255)
                                                                                     agaagaGGTTAATTGGTTTTTTGGTCTTAGTTCCTGCTTTTATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAAt                                               <  6:209888/139‑1 (MQ=255)
                                                                                            gTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAAt                                         >  6:251849/1‑138 (MQ=255)
                                                                                            gTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAAt                                         >  8:232970/1‑138 (MQ=255)
                                                                                            gTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAAt                                         >  5:3141/1‑138 (MQ=255)
                                                                                            gTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAAt                                         >  8:409411/1‑138 (MQ=255)
                                                                                              tAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATaa                                       >  7:387236/1‑138 (MQ=255)
                                                                                                      tttttGGTCTTAGTTCATTCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAAAAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCg                               <  3:362142/138‑1 (MQ=255)
                                                                                                       ttttGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCgg                              >  5:139869/1‑138 (MQ=255)
                                                                                                          tGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAAt                          >  4:334905/1‑139 (MQ=255)
                                                                                                           ggTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAAt                          <  6:139869/138‑1 (MQ=255)
                                                                                                           ggTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAAt                          <  7:164081/138‑1 (MQ=255)
                                                                                                           ggTCTTAGTTCCTGCTTTGATTATGTCAGGTATTACTTTAATCGAAGCAAATAAAAAGTCCCCTTTGGAAGTAATTGATAATATTCGGGATGAATTTGGTATATTCAGTGTGCAAATTGGTCAAACAGATCCGGCAAt                          <  7:348606/138‑1 (MQ=255)
                                                                                                                tAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACaa                     <  7:17651/138‑1 (MQ=255)
                                                                                                                 aGTTCCTTCTTTTATTAGTTGGGGTATAAATTTATTCTAATAAAACAAAACGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACaa                     <  6:269135/137‑1 (MQ=255)
                                                                                                                 aGTTCCTGCTTTGATTATGTGTTTAATAAATTTAATCGAATAAAACAAAACGACCCCTGTATAAGTTTTAGAAAGTTCTTGGGATGAATTTTGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAAt                    <  1:116813/138‑1 (MQ=255)
                                                                                                                 aGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACaa                     >  3:208950/1‑137 (MQ=255)
                                                                                                                 aGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACaa                     <  8:99413/137‑1 (MQ=255)
                                                                                                                       tGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAATAGGCATg             >  1:370630/1‑139 (MQ=255)
                                                                                                                            tGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAATAGGCATGGAt          <  8:122569/137‑1 (MQ=255)
                                                                                                                            tGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAATAGGCATGGAt          <  3:180509/137‑1 (MQ=255)
                                                                                                                              aTTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAATAGGCATGGATCaaa      >  6:199508/1‑139 (MQ=255)
                                                                                                                                   gtgGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAATAGGCATGGATCAAACCaa  >  3:149671/1‑138 (MQ=255)
                                                                                                                                      |                                                                                                                                      
ATAAAATATAAAGGTCAATAATGATATCCGTAGTATTAATAAAGGAGAGATTCTTTTCGACGAGGGAGAGTGTGATGGTGAATGAAGAAGAGGTTAATTGGATTTTTGGTCTTAGTTCCTGCTTTGATTATGTGGGGTATAACTTTAATCGAATCAAACAAAAAGACCCCTGTAGAAGTTTTAGAAAGTGCTTGGGATGAATTTGGGTTATTTAGTTTTGAAATTGGAATAACAGATCCGGCAATTACAATAGGCATGGATCAAACCAA  >  CP006881/2017216‑2017484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑