Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,014,0230TG43.1% 85.5 / 14.3 58G255G (GGT→GGG) pqiBparaquat‑inducible, SoxRS‑regulated MCE domain protein
Reads supporting (aligned to +/- strand):  ref base T (10/23);  new base G (25/0);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 1.13e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGG  >  NC_000913/1013884‑1014163
                                                                                                                                            |                                                                                                                                             
aCCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGG‑GCTGGATTTAgg                                                                                                                                    >  2:409526/1‑151 (MQ=255)
   gACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGATTTTggg                                                                                                                                   >  2:922901/1‑149 (MQ=255)
   gACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGGCAGCTTTTTTGGGCCGGAAGG‑GCTGGATTTTggg                                                                                                                                   >  1:85442/1‑149 (MQ=255)
       ggTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCa                                                                                                                             <  1:596041/151‑1 (MQ=255)
          gACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGt                                                                                                                           <  1:1031903/150‑1 (MQ=255)
          gACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCCTTTATGTGGCGGGAGG‑GCTGGGTTTAGGGCAGCCAGTg                                                                                                                          >  2:839951/1‑151 (MQ=255)
           aCCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGTAGCTTTGATGTGGCGGGAGG‑GCTGGTTTTTGGGCAGCCAg                                                                                                                            >  1:293730/1‑148 (MQ=255)
             caacaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGTGAGCTTTTATGTGGCGGGAGG‑GCTGGGGTTTAGGCAGCCAGTgg                                                                                                                         >  1:941063/1‑149 (MQ=255)
              aacaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTGAACTTTTATGTGCCGGGAAG‑TGTGGGTTTTTGGCAGCCCGTggggc                                                                                                                      >  2:594165/1‑148 (MQ=255)
               acaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCCGCTTTTATGTGGCGGGAGG‑GCTGGGTTTTGGGCAGCCCGTGGGAcc                                                                                                                     >  2:824348/1‑151 (MQ=255)
                caaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTTGGGGAGCCAGTGGCAcc                                                                                                                     >  2:44688/1‑150 (MQ=255)
                caaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCg                                                                                                                    <  2:258034/151‑1 (MQ=255)
                 aaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGG‑GCTGGGTTTTGGGGAGCCCGTGGCAcc                                                                                                                     >  2:498789/1‑149 (MQ=255)
                    gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGaaaa                                                                                                                <  2:186015/151‑1 (MQ=255)
                     ttCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTGATGGGCCGGGAGG‑GCTGGGTTTAGGGCAGCCAGGGGGACCGaaaa                                                                                                                >  2:841804/1‑150 (MQ=255)
                       cGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACa                                                                                                              <  1:524277/150‑1 (MQ=255)
                       cGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAg                                                                                                             <  2:976445/151‑1 (MQ=255)
                           tCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCCTTTGTGTGCCGGGAGGTTTTGGGTTTTGGGGGGCC‑GGGGCCCCCAAAACAACttt                                                                                                         >  2:517056/1‑151 (MQ=255)
                              ggAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGGAGG‑GCTGGGTTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCg                                                                                                       >  2:727552/1‑150 (MQ=255)
                                   gATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGGCGGGAGG‑GGTGGGGTTAGGGGGGCCAGTGGCACCGGAAAAAGCTTTCGtttt                                                                                                   >  2:607937/1‑149 (MQ=255)
                                    aTAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGGCGGAAGG‑GGTGGATTTTGGGGAGCCAGGGGGACCGAAAACAACTTTTGtttt                                                                                                   >  1:152347/1‑148 (MQ=255)
                                     tAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTTATGTGGCGGGAGG‑GCTGGGTTTTGGGGAGGCAGGGGGACCGGAAAAAGCTTTTGtttttt                                                                                                 >  2:588969/1‑147 (MQ=255)
                                        tGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTGCGCTTTTTTGGGCCGGAAGG‑GCTGGGTTTTGGGGAGCCAGGGGCACCGCAAAAAACTTTCGTTTTTTat                                                                                               >  1:10346/1‑148 (MQ=255)
                                          gCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGGCTGCTTTTATGTGCCGGTAGG‑TCTGGCATTAGGGGAGGCAGCGGGCCCGAAAAAAGCATTTGTtttgtttgtt                                                                                            >  2:670779/1‑147 (MQ=255)
                                           cATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATgatgat                                                                                         <  1:301021/151‑1 (MQ=255)
                                           cATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCCGAAGG‑TCTGGATTTAGGGCAGCCCGGGGCACCGAAAACAGCCTTCGtttttt                                                                                                 >  2:337647/1‑141 (MQ=255)
                                                cGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTTGGGGAGCCCGGGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCag                                                                                      >  2:904472/1‑149 (MQ=255)
                                                    tGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCCGCCTTTATGTGCCGGAAGG‑GCTGGGTTTAAGGCAGCCCGGGGGCCCGAAAAAAGCTTTTGTTTTTTATGATGATTagaa                                                                                    >  2:753651/1‑147 (MQ=255)
                                                          gACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCa                                                                              <  2:830370/146‑1 (MQ=255)
                                                            cGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTTAGCTTTTTTGTGCCGGGAGG‑GCTGGGTTTAGGGGAGCCCGGGGGACCGGAAACAGCTTTCGTTTTTTTTGATTATCAGAAGA‑GCATTCAgg                                                                        >  2:700365/1‑150 (MQ=255)
                                                             gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGTCCGCGTTTATGTGCCGGAAGG‑GCTGGATTTAGGGGAGCCCGAGGCGCCGAAAACAGCTTTCGTTTTGTATGATGTTAAGCAGA‑GCATTCAgg                                                                        >  1:780962/1‑149 (MQ=255)
                                                                 gcagGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑GCTGGATTTAGGGCAGCCCGTGGGACCGGAAACAACTTTTGTTTTTTATGATGATCAGAAGA‑GCATTCAGGATTCg                                                                   >  2:162187/1‑150 (MQ=255)
                                                                  cagGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGtt                                                                 <  2:394633/151‑1 (MQ=255)
                                                                   agGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTg                                                                <  1:1029751/151‑1 (MQ=255)
                                                                       aTGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGTGGCTTTTATGGGCCGGGAGG‑TCTGGGTTTTTGGGCGCCCGGGGGACCGAAAAAAGCCTTCGTTTTTTATGATGATTAAAAAA‑GGACTCAGGGATCGtt                                                                 >  1:159167/1‑146 (MQ=255)
                                                                         gCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACAcc                                                          <  1:953918/151‑1 (MQ=255)
                                                                          cGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTTGCGGTGTCAGCTCTGATGTGCCGGAAGG‑TGTGGATTTAGAGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTAGGATGATCAGAAGA‑GCATTCAGGATTCGTTGTAcccc                                                          <  1:449183/150‑1 (MQ=255)
                                                                             gtgGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGAt                                                       <  2:143040/150‑1 (MQ=255)
                                                                               ggAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATcac                                                    <  1:387753/151‑1 (MQ=255)
                                                                               ggAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTTGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTTTATGATGATCAGAAGA‑GCATTCCGGGTTCGTTGTACACCCATca                                                     >  2:271144/1‑150 (MQ=255)
                                                                                     gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGAt                                              <  2:921850/151‑1 (MQ=255)
                                                                                      ggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATTTGGCCGGAGG‑GCTGGGTTTTGGGCAGCCAGGGGCACCGAAAAAAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGAtt                                             >  2:571939/1‑151 (MQ=255)
                                                                                      ggCTCATTGACAACGCTGCTGAGTGGCGGGGTCAGCTTTTATGTGCCGGGAGG‑GCTGGGTTTAGGGGAGGCAGGGGGACCGGAAAAAGCTTTTGTTTTGTATGATGATCAGAAGAGGCTTCCAGGA‑TCGGTGTACACCGactcccattttt                                             >  1:623750/1‑140 (MQ=255)
                                                                                         tCATTGACAACGCTGCTGAGTTGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATc                                          >  2:65445/1‑151 (MQ=255)
                                                                                         tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGAAGG‑GCTGGGTTTTGGGGAGCCCGTGGCACCGAAAAACACTTTCGGTTTTTTTGATGATCAGAAGA‑GCATTCAAGGTTCGTTGTTCAACGGTCACATTGTTTAt                                           >  2:509980/1‑150 (MQ=255)
                                                                                              gACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGG‑GCTGGATTTAGGGCAGCCAGGGGGACCGAAAACAGCTTTCGTTTTTTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCCATCACATTTATTATCTGAt                                      >  2:540328/1‑150 (MQ=255)
                                                                                                    gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt                               <  2:189045/151‑1 (MQ=255)
                                                                                                               gCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGt                    <  1:32136/151‑1 (MQ=255)
                                                                                                                cGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGt                    >  1:772396/1‑150 (MQ=255)
                                                                                                                cGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTa                   <  2:335369/151‑1 (MQ=255)
                                                                                                                 ggTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAc                  >  2:272859/1‑151 (MQ=255)
                                                                                                                        gCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCt           <  1:867535/151‑1 (MQ=255)
                                                                                                                            tGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCa        >  2:152755/1‑150 (MQ=255)
                                                                                                                             gATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAc      >  1:986055/1‑151 (MQ=255)
                                                                                                                             gATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAc      <  1:922541/151‑1 (MQ=255)
                                                                                                                              aTGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAcc     <  2:816272/151‑1 (MQ=255)
                                                                                                                                  gCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCggg  <  1:36789/150‑1 (MQ=255)
                                                                                                                                  gCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCggg  <  1:286512/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
ACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGG‑TCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGA‑GCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGG  >  NC_000913/1013884‑1014163

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