Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,494	0	T	G	41.8%	104.7 / 12.9	55	intergenic (‑512/+91)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (2/30);  new base G (23/0);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.61e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987346‑987644                                                                                                                                                     |                                                                                                                                                         tcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGtt                                                                                                                                                         <  1:683729/150‑1 (MQ=255) tcttTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGtt                                                                                                                                                         <  1:955220/150‑1 (MQ=255)    ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTg                                                                                                                                                        <  2:964575/148‑1 (MQ=255)    ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑a                                                                                                                                                     <  1:935130/150‑1 (MQ=255)             ccAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGG‑aaa                                                                                                                                                   >  2:488030/1‑143 (MQ=255)               aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCtt                                                                                                                                         <  1:816718/151‑1 (MQ=255)                  aTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAc                                                                                                                                      <  2:904784/151‑1 (MQ=255)                   tGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑CAAAAAATTGCtttgct                                                                                                                                     >  1:965061/1‑147 (MQ=255)                    gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGG‑aaaaaaaattatttactt                                                                                                                                    >  2:231050/1‑136 (MQ=255)                       aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGT‑GAATAAAATGTTTgtgttt                                                                                                                                   >  2:323826/1‑144 (MQ=255)                       aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑AAATAAATTTATTGAt‑tt                                                                                                                                   >  2:396590/1‑148 (MQ=255)                       aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑AAATAAATTGTTTGAA‑tt                                                                                                                                   >  2:1043283/1‑148 (MQ=255)                        aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTaaca                                                                                                                               <  2:920199/151‑1 (MQ=255)                               gTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGAt                                                                                                                        <  2:212293/151‑1 (MQ=255)                                 aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATa                                                                                                                       <  1:259332/150‑1 (MQ=255)                                 aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGG‑AAATAGAATGTTTTTC‑TTtacaaaaaaaaa                                                                                                                       >  1:487157/1‑144 (MQ=255)                                    aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGTTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGCGTACC‑GAACACCCTACATGGC‑TTAACAAACAGaccgaa                                                                                                                    >  1:169623/1‑145 (MQ=255)                                       aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGC‑GAATAGACAGGTTGAT‑TTaaaagaaaaataaatt                                                                                                                   >  2:628744/1‑134 (MQ=255)                                            gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGaaa                                                                                                           <  2:334076/151‑1 (MQ=255)                                            gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGaaa                                                                                                           <  2:743282/151‑1 (MQ=255)                                                ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑AGAAAAATAGGTTGAC‑TTTACAAAAAAAAAAATTTACACaaaa                                                                                                          >  1:850056/1‑148 (MQ=255)                                                ttAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑AAAAAAATTTTTTTAC‑TTTACAAAAAAAAAAATTTACAGaaaaaa                                                                                                        >  1:935134/1‑148 (MQ=255)                                                      aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATaaaaaa                                                                                                 <  1:576804/151‑1 (MQ=255)                                                      aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATaaaa                                                                                                   <  1:975758/149‑1 (MQ=255)                                                       aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑AAAAAGAATTCTTTAA‑TTTAAAAAAAAAAAAATTGACAGAAAAAATaaaaaa                                                                                                 >  2:524614/1‑150 (MQ=255)                                                         aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                                                                                              <  2:856391/151‑1 (MQ=255)                                                          aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag                                                                                             <  2:364736/151‑1 (MQ=255)                                                           aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                            <  1:653910/151‑1 (MQ=255)                                                           aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                            <  1:560832/151‑1 (MQ=255)                                                            aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                           <  1:226896/151‑1 (MQ=255)                                                             aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                          <  1:1019586/151‑1 (MQ=255)                                                             aaCTAATCCGCAATAAGAATGAGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGT‑AAAAAAATTTTTTTAT‑TTAAAAAAAAAAAAAAATTACAGaaaaaaaaaaaaaaaa                                                                                              >  2:302421/1‑145 (MQ=255)                                                                  aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGT‑GAATAAAATTCTTTAA‑TTAAAAAAAAAAAAACATTACAGAAAAAATaaaaaaaaaaacaacggg                                                                                     >  2:499196/1‑140 (MQ=255)                                                                    ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                                   <  2:1051577/151‑1 (MQ=255)                                                                    ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                    <  2:919277/150‑1 (MQ=255)                                                                      gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGG‑AAAAAAATTGCTTGAC‑TTTACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCggg                                                                                     >  1:1040858/1‑147 (MQ=255)                                                                       cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑GAAAAAATTGGTTGGCTTTAAAAAACAAAAAAATTATCCGAAAAGAAAAAAAAAAAG‑CAGCggggggt                                                                                 >  2:689113/1‑146 (MQ=255)                                                                       cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑AAAAAAAATTCTTGTC‑TTTACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCggg                                                                                     >  1:345885/1‑146 (MQ=255)                                                                        aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg                                                                               <  2:260222/151‑1 (MQ=255)                                                                          tAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTGTTTACT‑TTAAAAAAAAAAAAAATTTAAAGaaaaaaaaaaaaaaa                                                                                               >  2:543790/1‑133 (MQ=255)                                                                           aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGG‑AAATAAAATGTTTTTC‑TTAAAAAAAAGATAAAATTAAAGCAAAAATAAAAAAAAAGACAACGGgggcgg                                                                                >  2:738552/1‑145 (MQ=255)                                                                            aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTG‑AAAAAAATTGCTTTTT‑TTTAAACAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAGACAGAGGGCGCggg                                                                               >  1:586019/1‑147 (MQ=255)                                                                                   ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                    <  1:161448/151‑1 (MQ=255)                                                                                    gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACg                                                                   <  2:1064281/151‑1 (MQ=255)                                                                                             ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGG‑AAATAAATTGATTGTC‑TTTACAAACAAATAAATTGACAGAAAAAATAAAAAAAAAAACAACGGGCGCCGGGTCTTGGAGAAGCTcaa                                                              >  2:975419/1‑147 (MQ=255)                                                                                               gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                         <  2:684164/150‑1 (MQ=255)                                                                                                     aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTgg                                                   <  1:168608/150‑1 (MQ=255)                                                                                                     aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                  <  1:557848/151‑1 (MQ=255)                                                                                                         gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                               <  1:341559/150‑1 (MQ=255)                                                                                                                 gTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGC‑AAAAAAAATTTTTTAC‑TTAAAAAAAAAATAAATTGACAGAAAAGAAAAAAAAAAAACAAGCGGGGGCCGGGTCTTTTAGAAGCTAAATTAAAAAAAGGGGGTACGCgg                                         >  2:882341/1‑148 (MQ=255)                                                                                                                   aaTATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTAc                                    <  1:489647/151‑1 (MQ=255)                                                                                                                           gTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                            <  1:927246/151‑1 (MQ=255)                                                                                                                              aCAGGATGATCGTTATGGCGGGGGTGG‑AAAAAAATTTGTTGAT‑TTAACAAAAAAATAAAATGACAGGAAAGAAAAAAAAAAAAACAACGGGGGGCGGGTCTTTTAGAAGCTAAAATATAAAAAAGCGTTAAGCGGAGTAAGgggggaggg                          >  2:227516/1‑145 (MQ=255)                                                                                                                              aCAGGATGATCGTTATGGCGGGGGTGC‑AAAAAAATTGTTTTTC‑TTAACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAGCCCGCGGGGGCCGGGTCTTGAAGACGCCAAAATAAAAACTGGCGGTTTGCGGAGTACGTggg                               >  1:949889/1‑145 (MQ=255)                                                                                                                                     gatCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                  <  1:937120/151‑1 (MQ=255)                                                                                                                                                  gggTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg     <  1:281866/151‑1 (MQ=255)                                                                                                                                                    gTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc    >  1:1019705/1‑150 (MQ=255)                                                                                                                                                     ttGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc  >  1:819427/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                         TCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGC‑AAATAGATTGCTTGAC‑TTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAC  >  NC_000913/987346‑987644

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