Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA44.4% 85.6 / 52.9 81intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (36/9);  new base A (0/36);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 7.58e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAA  >  NC_000913/5012‑5308
                                                                                                                                                      |                                                                                                                                                     
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAAt                                                                                                                                                    <  2:1114512/151‑1 (MQ=255)
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAAt                                                                                                                                                    >  2:369277/1‑151 (MQ=255)
   caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaa                                                                                                                                                  <  2:1172898/150‑1 (MQ=255)
      tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                              >  1:931192/1‑151 (MQ=255)
      tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                               >  1:839441/1‑150 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                             >  1:1124824/1‑151 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                             <  1:660349/151‑1 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                             <  1:873043/151‑1 (MQ=255)
        aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                            >  1:513566/1‑151 (MQ=255)
         aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                           >  1:358155/1‑151 (MQ=255)
           tCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccgccg                                                                                                                                         >  1:395659/1‑151 (MQ=255)
              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                          <  2:749409/147‑1 (MQ=255)
              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                          <  1:1208231/147‑1 (MQ=255)
                tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                    >  1:602010/1‑151 (MQ=255)
                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                    >  2:9450/1‑150 (MQ=255)
                  aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCAt                                                                                                                                  >  1:583292/1‑151 (MQ=255)
                      tctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                                >  1:963231/1‑149 (MQ=255)
                          aaTCGGGCGGGTTTTGTTTTTATGCAGCCCGTTTTTTTTATAAAAAAATTAGGGAAAAAAATGACAGGA‑AAAAA‑GGAAAAATTCTCAATAAAGGGGTTAATTT‑AGAGATTAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCg                                                                                                                                <  1:177808/145‑1 (MQ=255)
                           aTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                         >  1:475502/1‑151 (MQ=255)
                            tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                        >  2:484902/1‑151 (MQ=255)
                               ggCGGGGTTTGTTTTTATGCACCCCGGTTTTTTTAAAAAAAAATTAGGGAAAAAAATGAAAGGA‑AAAAA‑GAAAAAATTCTAAAAAAAGGCGGTAACTT‑AGAAATTAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                      <  2:884691/150‑1 (MQ=255)
                                   gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                  >  1:1251786/1‑150 (MQ=255)
                                      tttttttttATGCAGCCCGGTTTTTTTATGAAGAAATTATGGAAAAAAATGACAGGA‑AAAAA‑GGAAAAATCCTAAAAAAATGCGGTAACTT‑AGAGTTTAGGATGGGGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                                <  1:945853/144‑1 (MQ=255)
                                             ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                        >  2:447396/1‑150 (MQ=255)
                                             ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTCCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                        <  1:1202595/150‑1 (MQ=255)
                                              tATGCAGCCCGGTTTTTTTAAAAAAAAATTAGGGAAAAAAAGGACAGGA‑AAAAA‑GGAAAAATTCTCAATAAAGGCGTAAATTT‑AGAAATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                       <  1:442460/150‑1 (MQ=255)
                                                     ccgggtttttttATAAAAAAATTATGGAAAAAAAGGACAGGA‑AAAAA‑GGAAAAATTCTAAAAAAAGGCGGTAATTT‑AGAGATTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                               <  1:445899/145‑1 (MQ=255)
                                                        gttttttttATAAAAAAATTAGGGAAAAAAATGACAGGA‑AAAAA‑GGAAAATTTCCAAAAAAATGCGTTAATTA‑AGATATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                              <  2:813006/146‑1 (MQ=255)
                                                        gtttttttaATAAAAAAATTAGGAAAAAAAAAGACAGGA‑AAAAA‑GGAAAAATTCTAAAAAAAGGCGTTAACTA‑AGAGATTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                             <  2:1258533/147‑1 (MQ=255)
                                                          tttttttATAAAAAAATTATGAAAAAAAAGGAAAGGA‑AAAAA‑GGAAAAATTCACAAAAAATGCGATAATTT‑AGAAATTGGGATTGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                            <  2:281580/148‑1 (MQ=255)
                                                           ttttttATAAAAAAATAAGGAAAAAAAAGGACAGGA‑AAAAA‑GAAAAATTCCAAAATAAAGGCGGTAACTA‑AGAAATAAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                            <  2:1127629/148‑1 (MQ=255)
                                                            tttttATAAAAAAATTATGGAAAAAAAGGACAGGA‑AAAAA‑GGAGAAATTCTAAAAAAAGGCGTTAACTA‑AAAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                         <  1:503634/150‑1 (MQ=255)
                                                             ttttATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                        >  1:287309/1‑150 (MQ=255)
                                                             ttttATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                        >  1:996502/1‑150 (MQ=255)
                                                             tttaaaaaaaaaaTAAGGGAAAAAAATGAAAGGA‑AAAAA‑GGA‑AAATTCTAAAAAAAATGGGTAAATTAAGATTTTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                        <  2:501039/140‑1 (MQ=255)
                                                               ttaaggaagaaTTTTGGAAAAAAAATGCAGGGAAAAAA‑AGAAAAATTCTTAATAATGGCGTTAATTT‑AGAGATTAGGATTGCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                       <  1:777273/141‑1 (MQ=255)
                                                                 agaaagaaATTATGGAAAAAAATGACAGGA‑AAAAA‑GGAAAAATTCTAAAAAAAGGCGGTAATTA‑AGAAATTAGGATGGGGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                    <  1:550386/147‑1 (MQ=255)
                                                                 aTGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                     >  1:1188611/1‑149 (MQ=255)
                                                                    aggaaATTATGAAAAAAAATGACAGGA‑AAAAA‑GAAAAAATCCTAAAAAAATGGGGTAATTT‑AGAGATTAGGATTGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                     <  1:1182848/144‑1 (MQ=255)
                                                                    aaaaaaTTAGGAAAAAAAATGACAGGA‑AAAAA‑GAAAAAATTCTAAAAAAAGGCGGTAATTA‑AGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                  <  2:669642/146‑1 (MQ=255)
                                                                      gaaATTATGGCGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                              >  2:462967/1‑151 (MQ=255)
                                                                        aaTTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                            >  2:1225397/1‑151 (MQ=255)
                                                                          ttATGAAAAAAAAGGAAAGGA‑AAAAA‑GAAAAAATTCTAAAAAAGGGCGTTAATTT‑AGAGATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                           <  1:1047402/150‑1 (MQ=255)
                                                                          taagagaaaaaaaaGACAGAA‑AAAAA‑AGAAAAATTCACAAAAAATGCGAAAATTA‑AGAAATGAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                            <  2:296424/141‑1 (MQ=255)
                                                                            aTGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                         >  1:26393/1‑150 (MQ=255)
                                                                                aaaaaaaTGAAAGAG‑AAAAA‑GGAAAAATTCTCAAAAATTGCGGTAACTT‑AGAGATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaaa                                                                     <  2:1059263/148‑1 (MQ=255)
                                                                                 aaaaaggaaCAGGG‑AAAAA‑GAAAAAACTGCAAATAAGTGCTTTAACTG‑ATAAATTGGGATGGAAGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                   <  2:666011/143‑1 (MQ=255)
                                                                                  aaaaaTGAAAGGA‑AAAAA‑GAAAAAATCCAAAAAAAAGGCGTTAATTA‑AGAGATTGGGTTGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                  <  1:1202237/151‑1 (MQ=255)
                                                                                        gACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                            >  1:20249/1‑151 (MQ=255)
                                                                                        gACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                            >  1:218612/1‑151 (MQ=255)
                                                                                           aggg‑aaaaagggaaaaTTTCACAAAAAGTGCGATAACTT‑AGAATTTGGGATGGCAGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                          <  2:439082/151‑1 (MQ=255)
                                                                                             ggaaaaaa‑GGAAAAATTCTCAAAAAATGGGGAAATTT‑AGAGATTGGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                         <  1:272640/147‑1 (MQ=255)
                                                                                             ggaaaaaa‑GGAAAAATTCTAAAAAAAGGGGGTAATTT‑AGAGATTAGGATTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                           <  1:861882/145‑1 (MQ=255)
                                                                                             ggaaaaaa‑GGAAAAATTCAAAATAAATGCGTTAACTG‑AGATATTAGGTTTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                          <  2:279736/146‑1 (MQ=255)
                                                                                               aaaaaa‑GGAAAAATTCTAAAAAAAGGGGGTAACTT‑AGAAATTAGGTTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                       <  1:1235385/149‑1 (MQ=255)
                                                                                                aaaaa‑GGAAAAATTCTCAAAAAGGGCTTTAACTA‑GGATATTGGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATAAATCg                                                       <  1:132352/149‑1 (MQ=255)
                                                                                                aaaaa‑GAAAAAATTCAAAAAAAATGCGTAAACTT‑AGAGATGAGGATTGAGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                        <  2:1124810/148‑1 (MQ=255)
                                                                                                  aaa‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                   >  1:6970/1‑151 (MQ=255)
                                                                                                  aaa‑GGAAAAATTCAAAATAAATGGGGAAACTA‑AGAATTGGGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                   <  1:933272/151‑1 (MQ=255)
                                                                                                        aaaaaTTCTAAAAAAATGCGGTAACTT‑AGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                              <  1:532156/149‑1 (MQ=255)
                                                                                                         gaAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             >  1:167648/1‑151 (MQ=255)
                                                                                                         gaAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             >  2:1120757/1‑151 (MQ=255)
                                                                                                          aaaTTCTAAAAAAATGCGTTAACTT‑AAAGATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             <  1:708731/150‑1 (MQ=255)
                                                                                                                 cAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                     >  1:1234778/1‑151 (MQ=255)
                                                                                                                  aaaaaagggggTAATTA‑AGAGTTTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                     <  1:899924/141‑1 (MQ=255)
                                                                                                                   aaaaaggggtaaattt‑AGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                   <  1:507967/137‑1 (MQ=255)
                                                                                                                   aTAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                   >  1:329075/1‑151 (MQ=255)
                                                                                                                    tAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                  >  1:1182700/1‑151 (MQ=255)
                                                                                                                              aaCTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                        >  1:860912/1‑151 (MQ=255)
                                                                                                                                 tt‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                     >  1:66116/1‑151 (MQ=255)
                                                                                                                                      agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                  >  1:224891/1‑150 (MQ=255)
                                                                                                                                        aTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag                >  1:386651/1‑150 (MQ=255)
                                                                                                                                         ttAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc               <  1:129822/150‑1 (MQ=255)
                                                                                                                                          tAGAATTGAGGAAAATACAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca              <  1:1031244/150‑1 (MQ=255)
                                                                                                                                            ggATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca            <  2:1062823/150‑1 (MQ=255)
                                                                                                                                             gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg          >  2:1021965/1‑151 (MQ=255)
                                                                                                                                             gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg          >  1:753655/1‑151 (MQ=255)
                                                                                                                                               tgggggAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc               <  1:786581/140‑1 (MQ=255)
                                                                                                                                                 gggaagaAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGATGc      <  2:358084/147‑1 (MQ=255)
                                                                                                                                                 gCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc      <  1:344479/151‑1 (MQ=255)
                                                                                                                                                   ggAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg    >  1:791935/1‑151 (MQ=255)
                                                                                                                                                   ggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg       <  1:1168713/148‑1 (MQ=255)
                                                                                                                                                    gagaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg    >  1:118035/1‑150 (MQ=255)
                                                                                                                                                    gagaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg     >  1:908933/1‑149 (MQ=255)
                                                                                                                                                    gagaAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc      <  1:43592/148‑1 (MQ=255)
                                                                                                                                                      gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaa  >  2:1235683/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAA‑GGAGAAATTCTCAATAAATGCGGTAACTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAA  >  NC_000913/5012‑5308

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