Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA38.3% 75.2 / 25.3 60intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (29/8);  new base A (0/23);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 3.58e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5013‑5305
                                                                                                                                                      |                                                                                                                                                   
atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATa                                                                                                                                                 <  1:418360/151‑1 (MQ=255)
   aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATaa‑ca                                                                                                                                             <  2:1051707/151‑1 (MQ=255)
      aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATaa‑caac                                                                                                                                           <  1:486009/150‑1 (MQ=255)
      aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAAcc                                                                                                                                          >  2:478870/1‑151 (MQ=255)
                      ctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGt                                                                                                                          >  2:469360/1‑151 (MQ=255)
                          aTCAGGCGGGTTTTGTTTTTAGGCACCCGGGTTTTTTTATG‑AAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAT‑TTCTAAAAAAAGGCGTTAACTTAG‑AGATT‑AGGATTGGGGAGAAAAA‑CAACCGCCGTTCTCATCGAGTAAtc                                                                                                                      <  2:832030/151‑1 (MQ=255)
                              ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                  >  2:901554/1‑151 (MQ=255)
                                cgggttttGTTTTTAGGCAGCCCGGTTTTTTAATG‑AAAAAATTATGGAGAAAAATGAAGGGGAAAAAGGAAAAT‑TTCTAAAAAAAGGGGTTAACTTAA‑AATTG‑GGGATTGCGAAAAAAAA‑AAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                 <  1:295785/145‑1 (MQ=255)
                                   agTTTGCTTTTGTGCAGGCCGGCTTTTTTATG‑AAGAAATTATGGTGAAAAATGACAGGTAAAAATGAGTAA‑TTCTCAATAAATGCGGTAACTTAG‑AGAGT‑AGGAGTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                             >  1:481999/2‑151 (MQ=255)
                                     ttttCTTTAATGCAG‑CCGGTTTTTTTAAGAAAAAAATTAGGGAAAAAAAGAAGGGGAAAAAAGGAAAAA‑TTCACAAAAAATGCGTTCTCATAG‑AGTTG‑AGGTGGGCGAAAAAAAA‑AAACCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                              <  1:442570/148‑1 (MQ=255)
                                     tttGTTTTTAAGCAGCCCGTTTTTTTTATA‑AAAAAATTATGGAAAAAAAAAACAGGAAAAAAGGAGAAA‑TTCTAAATATATGGGTTAACTTAG‑AAATA‑AGGATTGCGGAAAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                            <  2:580127/150‑1 (MQ=255)
                                         cTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                       >  1:169440/1‑151 (MQ=255)
                                           tttaggcgtcccGTCTTTTTTAAG‑AAAAAATTTTGAAAAAAATGGAAAGGAAAAAAAGAGAAA‑TCATCAAAAACGGGGTTATTTTAG‑GGATT‑ATGGTGGCGAAAAACAA‑AAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                      <  2:532310/141‑1 (MQ=255)
                                             tagcaaGCCGGGTTTTTTTAAA‑AAAAAATTATGGAAAAAAAAGACAGGAAAAAAGGAAAAA‑TTCTAAATAAATGCGTTAACTAAG‑AGATT‑AGGATGGGGGAAAAAAA‑AAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                    <  2:10977/145‑1 (MQ=255)
                                             tATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                    >  2:795496/1‑150 (MQ=255)
                                              aTGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                   >  1:950117/1‑150 (MQ=255)
                                               gacaccctgttttttttAGA‑AAAAAATTA‑GGAAAAAAAAGAAAGGAAAAAAAGAAAAACTTCTCAATAATTGCGTTAAATTAA‑AATTT‑AGATTGGCGAAAAAAAA‑CACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                 <  1:493377/140‑1 (MQ=255)
                                                  aGCCCGTTTTGTTTATA‑AAAAAATTATGGAAAAAAAGGAAAGGAAAAAAGGAAAAA‑TTATCAATAAGTGCGGTAAATAAG‑AGATT‑AGGAGGGCGAAGAATAC‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                  <  1:940768/147‑1 (MQ=255)
                                                  aGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                              >  2:684153/1‑151 (MQ=255)
                                                    ccgggtttttttAAA‑AAGAAATTAGGAAAAAAAATGACGGGAAAAAAGAAAAAA‑TTCTAAATAAACGCGGTAACTAAG‑AGAGA‑AGGATGGGGGAAAAAAA‑CAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                 <  1:302050/140‑1 (MQ=255)
                                                         tttttttAAA‑AAAAAATTAGGAAGAAAAAGGACGGAAAAAAAGGAAAAA‑TTCTAAAAAAATGGGTTAACTAGG‑AGATT‑AGGATTGCGAAAAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                         <  1:1090574/148‑1 (MQ=255)
                                                            ttttATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                     >  2:886386/1‑150 (MQ=255)
                                                            ttttAGA‑AAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAT‑TTCAAAAAAAATGCGTTAATTTAG‑AGATT‑AGGATTGGGAAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                     <  2:693075/150‑1 (MQ=255)
                                                                 tG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                               >  1:40159/1‑151 (MQ=255)
                                                                    aaaaaaTTAGGAAAAAAAATGAAAGGAAAAAAGGAAAAA‑TTCTAAAAAAAGGCGGTAACTTAG‑AGATT‑AGGATGGGGGAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                               <  2:602446/146‑1 (MQ=255)
                                                                       aaaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                          >  1:93680/1‑151 (MQ=255)
                                                                        aaTGATGGAAAAAAAGGACGGAAAAAAGGGAGAAA‑TTCCAAATAAGTGCGTTAA‑TTAA‑AAATTGAGTAGGGGTGAGAAAAA‑AAACCGCCGTTCTCAGCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                         <  1:301904/151‑1 (MQ=255)
                                                                          ttATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                       >  2:1079367/1‑151 (MQ=255)
                                                                          ttATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                       >  1:449245/1‑151 (MQ=255)
                                                                            aTGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                      >  2:625068/1‑150 (MQ=255)
                                                                             ggaagaAAAAGCAAAGGAAAAAAAGAAAAC‑TTCTAAAAAAATGCGGAATTTTAG‑ATATG‑A‑GATTGGGGAAAAAAACCACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGTCAATGATAAAAGGAGTAACCtggg                                                                       <  1:230151/145‑3 (MQ=255)
                                                                                 aaaaaaGGAGAGGAAAAAAGGAAAAA‑TCCTCAAAAAAGGCGTAACCTAAG‑AAATT‑AGAATTGAGGAAAAAAA‑AACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                   <  1:451991/147‑1 (MQ=255)
                                                                                  aaaaaTGAGGGAAAAAAGGGAAAAA‑TTCTAAAAAAAGGGGGTAATTAAG‑AAATT‑AGGTTGGGGAAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                <  2:697503/150‑1 (MQ=255)
                                                                                  aaaaaTGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                   >  2:827380/1‑147 (MQ=255)
                                                                                  aaaaaTGACAGGAAAAAAGGATAAA‑TTCACAATAAAGGCGGTACCTTAG‑AGATT‑AGGATTGCGGAGAAAAA‑CACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa                                                                    <  2:405676/146‑1 (MQ=255)
                                                                                     aaagaaggggAAAAAGGAGAAA‑TTTTAAATAAGGGC‑GTTATTGAGAAGATT‑GGAATGGCGAAAAACAA‑CAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                            <  2:1127577/144‑1 (MQ=255)
                                                                                     aaTGAAAGGAAAAAAGGAAAAA‑TTCTAAAAAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                            <  2:800321/151‑1 (MQ=255)
                                                                                       tGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                          >  1:988559/1‑151 (MQ=255)
                                                                                             ggAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCTAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                      >  2:722908/1‑149 (MQ=255)
                                                                                               aaaaaGGAAAAT‑TTCAAAAAAACTCCGCTAATTTGA‑AGATT‑AGGATGGTGAAAAATAA‑GAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                       <  1:769566/146‑1 (MQ=255)
                                                                                                aaaaGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                 >  2:161404/1‑151 (MQ=255)
                                                                                                aaaaGGAAAAA‑TTCACAATAAAGGCGTTATCTAAG‑AGTTT‑AGGTTGGCGAAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                 <  1:280967/151‑1 (MQ=255)
                                                                                                       gaAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                          >  1:451993/1‑151 (MQ=255)
                                                                                                                 aaaaaaTGCGGTAATTTAG‑AGATT‑AGGATTGCGGAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                  <  2:987450/147‑1 (MQ=255)
                                                                                                                 aaTAAAGGCGTTAATTGAG‑ATATT‑GGGTTTGCGGAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                 <  2:1030871/151‑1 (MQ=255)
                                                                                                                 aaTAAAGGCGGAAACTTAG‑AGATT‑AGGATTGCGGAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                   <  2:16211/149‑1 (MQ=255)
                                                                                                                    aaaTGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                              >  2:384198/1‑151 (MQ=255)
                                                                                                                       tGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATAGAGAAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                            >  1:998845/1‑150 (MQ=255)
                                                                                                                             aaCTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                     >  2:812321/1‑151 (MQ=255)
                                                                                                                              aCTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                    >  1:142384/1‑151 (MQ=255)
                                                                                                                                 ttaaagttttgggATTGCGGAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                         <  2:1108855/134‑1 (MQ=255)
                                                                                                                                 tAG‑AGATT‑AGGATTGCGGAGAATAA‑GAACCGCCGTAATCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATACAAGGAGTAACCAGTGAAAAAGATGCAATCTATCGTACTCCCACTTTCCCGGGATCTGGTGGATCCCa                   >  1:49866/1‑149 (MQ=255)
                                                                                                                                 tAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                 >  1:561123/1‑151 (MQ=255)
                                                                                                                                   g‑agaTT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc               >  1:643703/1‑151 (MQ=255)
                                                                                                                                     agaTT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc               >  1:659468/1‑150 (MQ=255)
                                                                                                                                      agtgattgtatGCGAAAAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAta                 <  1:41897/138‑2 (MQ=255)
                                                                                                                                           aggtgggaggAAAAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                  <  1:309834/134‑1 (MQ=255)
                                                                                                                                            ggATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac          >  1:120651/1‑149 (MQ=255)
                                                                                                                                               tgggaaagaAAAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                            <  1:459955/122‑1 (MQ=255)
                                                                                                                                               tggCGGAAAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg       <  2:216224/147‑1 (MQ=255)
                                                                                                                                                  cGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg  >  1:731091/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATG‑AAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAACTTAG‑AGATT‑AGGATTGCGGAGAATAA‑CAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5013‑5305

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