Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	44.1%	85.3 / 25.4	59	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (30/3);  new base A (0/26);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 9.23e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436369‑436663                                                                                                                                                    |                                                                                                                                                    tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCCCTGGCGGTGAAGCTTACGCGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGCGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCa                                                                                                                                                    >  2:590469/1‑149 (MQ=255)           ttGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                         >  2:778539/1‑150 (MQ=255)               aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATc                                                                                                                                    >  2:162061/1‑151 (MQ=255)                ccGGACCAGGCGCTGCGCTGGGCGCTCTCGGGCGGTGAAGATTACGAGTTGTGTTTCCCTGTGCCGGAACTAACCCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCGGTTTACCTGTAGCGGGAAAATGACCGCCGATATc                                                                                                                                    <  2:193629/150‑1 (MQ=255)                  ggAACGGGCGCGGCGCGGGGCGCTCTCGGGGGGGGAAGATTCCGAGTTGTTTTTCACTGTGCCGGAACTGACCGGTGGCGCGCTGGAGGTGGCTCTCGCACCCCTGGGCGTCCCTTTTACCTGTATCGGGAAAATGACCGCCGATATCGaa                                                                                                                                 <  2:160610/151‑1 (MQ=255)                           gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                        >  2:882878/1‑151 (MQ=255)                             tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                      >  2:788637/1‑151 (MQ=255)                              gcgctgGGCGCTCTTGGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGGGGCGCGCGGGATGTGGCTCTCGGACACCTGGGCGTCCCGTTTCCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                       <  1:131500/149‑1 (MQ=255)                               cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                      >  2:739771/1‑149 (MQ=255)                                gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                   >  2:225908/1‑151 (MQ=255)                                 ctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                  >  1:12403/1‑151 (MQ=255)                                   gggCGCTCTCTGGGGGTAAAGATTACGAGTTGTTTTTCCTTGGGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTACCGTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                  <  1:256593/149‑1 (MQ=255)                                          ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                         >  2:222455/1‑151 (MQ=255)                                             tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                       >  2:57384/1‑150 (MQ=255)                                              gggggTGAAGTTTAGGATTTTTGTTTCCCTGTGCGGGAACTGACCGGGGGCGGGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCTTTTTACTTGTGGCGGGAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                        <  2:538332/145‑1 (MQ=255)                                                 ggagAAGATTCGGGGTTGTGTTTCAGGGTGCGGAAACTGACCCGTGGCGGGCTGGTTGGGGCTCTGGGCCACCTGGGCTTCCCTTTTCCCTGTAGCGGACAAAGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                  <  2:304454/148‑1 (MQ=255)                                                 ggTGAAGTTTACGATTTGTGTTTCCCTGTGCCGGAACTGAACGGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTCCCTTTTCCCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                   <  1:394945/150‑1 (MQ=255)                                                 ggTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                  >  1:599293/1‑151 (MQ=255)                                                 ggTAAAGATTACGAGTTGTGTTTCATTGTGCCGGAACTGACCCGGGGCGCGTTGGATGTGGCTCTCGGACACCTGGGCGTCCCGTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                   <  1:954066/150‑1 (MQ=255)                                                   tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                >  1:1067314/1‑151 (MQ=255)                                                     agatttACGGGTTGTTTTTCCCTGGGCGGAAACGGACCCGGGGCGCGTTGGTGGGGGCTCGCGGCCCCCGGGGCGTTCTGTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                               <  2:413660/146‑1 (MQ=255)                                                     aaGATTACGAGTTGTGTTTCACTGTGCCGGAGCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                              >  2:236118/1‑151 (MQ=255)                                                             gAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                      >  1:670552/1‑151 (MQ=255)                                                             gAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                      >  2:880873/1‑151 (MQ=255)                                                              aGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                     >  1:902150/1‑151 (MQ=255)                                                                 tgtgTTTCCCTGTGCCGGAACGGACCCGTGGCGCGTTGGAGGTGGCTCTCGGACCCCTGGGCGTACCGTTTCCCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                   <  1:322700/150‑1 (MQ=255)                                                                   tgtTTCCCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCCTTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                                 <  1:295275/150‑1 (MQ=255)                                                                    gtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                               >  1:600668/1‑151 (MQ=255)                                                                     tttCACGGTGCGGGAACGGACCCGGGGCGCGTGGGTTGTGGCTCTCGGACCCCGGGGCGTCCCGTTTCCCTGTACCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                   <  1:191260/146‑1 (MQ=255)                                                                         aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                          >  2:519238/1‑151 (MQ=255)                                                                               ccGGACCTGACCCGTGGCGCGTTGGTTGTGGCTCTCGGCCCCCGGGGGGTCCCTTTTACCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                    <  1:1096157/151‑1 (MQ=255)                                                                               ccGGAACTGACCCGTGGCGCTTTGGATGTGCCTCTCGGACCCTTGGGCGTCCCTTTTACTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                    <  2:161974/151‑1 (MQ=255)                                                                                      tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                             >  1:757212/1‑151 (MQ=255)                                                                                            gaggcgctttggcggtgGCTCTGGGCCACTGGGGCGTCCGTTTTACTTGTTCCGGGAAAAGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                            <  2:914689/132‑1 (MQ=255)                                                                                            gTGGCGCGTTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                             <  1:575704/145‑1 (MQ=255)                                                                                              ggCGCGTGGGGTGGGTCTCTCGGACCCCTGGGCGTCCCTTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                      <  2:815530/150‑1 (MQ=255)                                                                                              ggCGCGCTGGTTGTGTCTCTCGCCCTCCTGGGCGCCCCTTTTACCTGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                         <  1:500586/147‑1 (MQ=255)                                                                                              ggCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                      >  2:410973/1‑150 (MQ=255)                                                                                               gcgcgcGGGATGTGCCTCTCGGCCACCGGGGCGTACCTTTTCCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAg                                                    <  1:383226/151‑1 (MQ=255)                                                                                                   gttGGATGGGGCTCCGGGCCCCCGGGGGGTCCCTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                 <  1:711242/148‑1 (MQ=255)                                                                                                     tGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                              >  1:106965/1‑151 (MQ=255)                                                                                                          ggggCTCCGGGACCCCGGGGCGTCCCTTTTCCCTGTTCGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGt                                                <  1:909596/142‑1 (MQ=255)                                                                                                            ggCTCTGGGTCACCGGGGGGTTCCTTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAg                                       <  1:364367/151‑1 (MQ=255)                                                                                                                 tcGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc                                  >  2:1025904/1‑151 (MQ=255)                                                                                                                  cGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                                 >  2:595014/1‑151 (MQ=255)                                                                                                                      ccgcTGGGCGTCCCTTTTACCTTTAGCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGCCTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                   <  2:152522/142‑1 (MQ=255)                                                                                                                       acCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                               >  1:180502/1‑148 (MQ=255)                                                                                                                          ggggCGTACTTTTTACCTGTAGGGGAAAAAGGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGa                         <  2:289550/150‑1 (MQ=255)                                                                                                                               gtccgttttgcttGTATGGGAAAAAAGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt                    <  1:241459/139‑1 (MQ=255)                                                                                                                                tACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                   >  2:186563/1‑151 (MQ=255)                                                                                                                                  ccGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                 >  1:114804/1‑151 (MQ=255)                                                                                                                                     tttccttttttggggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGt                <  1:528320/137‑1 (MQ=255)                                                                                                                                     tttccttgattCGGGAAAAGTACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg               <  2:422863/140‑1 (MQ=255)                                                                                                                                     tttACCTGTAGCGGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg               <  2:982610/150‑1 (MQ=255)                                                                                                                                        aCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCAt           >  1:1108042/1‑151 (MQ=255)                                                                                                                                           tGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt         >  1:98486/1‑150 (MQ=255)                                                                                                                                               tCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTg    >  2:123611/1‑151 (MQ=255)                                                                                                                                                 gggCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt  >  2:11246/1‑151 (MQ=255)                                                                                                                                                 gggAAAAAGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt     <  2:393277/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436369‑436663

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