Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,333,757	0	T	G	46.4%	84.4 / 17.3	57	G16G (GGT→GGG)	ispB	octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (4/26);  new base G (26/0);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.24e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCA  >  NC_000913/3333609‑3333906                                                                                                                                                     |                                                                                                                                                       gtgCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGgtgt                                                                                                                                                      <  2:123791/151‑1 (MQ=255)        tttttCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGCg                                                                                                                                               >  1:190497/1‑151 (MQ=255)          tttCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCgg                                                                                                                                              <  1:427887/150‑1 (MQ=255)               gTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATcc                                                                                                                                        <  1:754925/151‑1 (MQ=255)                   aaTCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGa                                                                                                                                    <  1:225890/151‑1 (MQ=255)                     tCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGCGGGAAACCttta                                                                                                                                    >  2:623459/1‑147 (MQ=255)                      cACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGagca                                                                                                                                 <  2:785447/151‑1 (MQ=255)                      cACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGagca                                                                                                                                 <  1:354601/151‑1 (MQ=255)                       aCCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGCGTTTATTGGGGAAACCtttggt                                                                                                                                  >  1:18608/1‑145 (MQ=255)                          cTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGCGGGAATTATTTAGGA‑GCtt                                                                                                                            >  1:748264/1‑151 (MQ=255)                            atatTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGCGGCAAACCTTTGGCA‑GCtt                                                                                                                            >  1:949742/1‑149 (MQ=255)                              atTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAAtt                                                                                                                        <  2:621282/151‑1 (MQ=255)                                 ccTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTcc                                                                                                                      <  1:391914/150‑1 (MQ=255)                                  cTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGAGGGAAACCTTTAGGA‑GCttta                                                                                                                          >  1:585903/1‑143 (MQ=255)                                    aaCCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACg                                                                                                                  <  1:205387/151‑1 (MQ=255)                                    aaCCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACg                                                                                                                  <  2:72866/151‑1 (MQ=255)                                          aaaCCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATTAGGCAAACCCTTAGCA‑GGTTAATTTCCACGTcc                                                                                                               >  1:216924/1‑148 (MQ=255)                                           aaCCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATGCGGCAAACCTTTTGCA‑GCCTAATTTTCACGTCCCAc                                                                                                            >  2:254021/1‑150 (MQ=255)                                            aCCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGTCGCAATTCTTTTGCA‑GCTTTATTTCGACGcccca                                                                                                             >  1:338356/1‑146 (MQ=255)                                             cccTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGa                                                                                                         <  2:885050/151‑1 (MQ=255)                                                 aaGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGatc                                                                                                       <  2:781054/149‑1 (MQ=255)                                                   gTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATAAGGCAAACCTTTAGCAGGCTTATTTCC‑ACGGCCAAATGatcaa                                                                                                     >  2:835687/1‑149 (MQ=255)                                                      gCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAg                                                                                                 <  1:244832/150‑1 (MQ=255)                                                         tttGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGTGGGAAACCTTTAGGA‑GGCTAATTTCGACGGCCCAATGgtcaaataat                                                                                                >  2:383743/1‑143 (MQ=255)                                                               cacaGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTAtt                                                                                       <  1:279278/151‑1 (MQ=255)                                                               cacaGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTTTTTGGCAATTCTTTAGCA‑GGCTAATTTCGACCTGCCACTTCTCCATCAATTTAGCTAtt                                                                                       >  1:299299/1‑151 (MQ=255)                                                               cacaGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGCGGGAAACCTTTTGCA‑GCCTTATTTCCACGGCCAACTGATCAAACAGTTATGCTAtt                                                                                       >  2:810172/1‑151 (MQ=255)                                                                 caGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGCGGGCATCCTTTGGCA‑GCTTTATTTCGACGGCCCAATGgtcaaacaagttgggtat                                                                                        >  1:186018/1‑138 (MQ=255)                                                                      aGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGt                                                                                <  2:238944/151‑1 (MQ=255)                                                                       ggTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTc                                                                               <  2:53217/151‑1 (MQ=255)                                                                         tAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTAATTCGGCGATCCTTTAGCA‑GCTTAATTCCGCCGTCCAACTGATCAATCAGTTAGGCTATTAAATCGTc                                                                               >  1:849029/1‑149 (MQ=255)                                                                                ggCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgcggcggc                                                                      <  2:626705/151‑1 (MQ=255)                                                                                 gCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgcggcggcg                                                                     <  1:432368/151‑1 (MQ=255)                                                                                   gAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGCGGCAAACCTTTAGCA‑GCTTAATTTCGACGGCCAACTGATCAAAAAATTAGGCTATTAAATCGTCAGcgggggcgg                                                                    >  2:909053/1‑150 (MQ=255)                                                                                   gAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAAAGCGGCAAACCTTTAGCA‑GCCCAATTTCCAAGTCCAACTGATCAATAAGATAAGGTATTTAATCGTCAGcggcggcgg                                                                    >  2:562112/1‑150 (MQ=255)                                                                                   gAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGAAAATGCGGCGATCCCTTGGCA‑ACGGCATTTCTACGTGCAAAAGGTCCATATATTAGGCTATTATCTCGCCAgcgggcg                                                                       >  1:698836/1‑144 (MQ=255)                                                                                         cccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGt                                                             <  2:120079/151‑1 (MQ=255)                                                                                         cccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGt                                                             <  2:172792/151‑1 (MQ=255)                                                                                              cTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGGTAATAAGGCAATCCTTTAGGA‑GCTTAATTTTCACGTCCAACTGATCAATCAGTTAGGGTATTTAATCTTCAGCGGCGGGGGGAAAccgatt                                                          >  2:268552/1‑144 (MQ=255)                                                                                               ttttGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGAGGCAAAACCTTAGCAGGCTTAATTCC‑ACGTCCCACCGATCAAACAGTTAGGCTATTACATCGTCAGCGGCGGGGGGGAACGGATTTGt                                                       >  2:963258/1‑151 (MQ=255)                                                                                                    cGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCgat                                                  <  1:872345/151‑1 (MQ=255)                                                                                                      aTGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATTAGGGAATACTTGAGGA‑GGTTTATTCCCACGGCCCAAAGATCAAACAAATAAGGGATTATATCCACAGcggggggggggaaaaaatttgt                                                       >  1:277397/1‑131 (MQ=255)                                                                                                        gAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGAtt                                              <  2:937030/151‑1 (MQ=255)                                                                                                          aTTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATGCGGCAATCCTTTAGCA‑GCCTAATTTCCACGTCCAAATGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGGAAAAGTATTCGTCCGCTGAttt                                             >  1:593141/1‑149 (MQ=255)                                                                                                            ttAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAAAGCGGGCATCCTTTAGCA‑GCTTTATTTCGACGGCCCACTGATCAAACAATTAGGCTATTACATCGTCAACGGCGGCGGGAAAAGGATTTGTCCGCTGAtt                                              >  1:827612/1‑147 (MQ=255)                                                                                                                  aaaaTCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATAATGCAAACCTTTAGCA‑GCTTAATTTCCACGGCCCACATATCAATCAATTAAGCTATTAAATCGTTAGCGGCGGCGGGGAAAGGAATTGTTCGCTGATTTCTGt                                         >  2:424720/1‑146 (MQ=255)                                                                                                                  aaaaTCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATTCGGCCATCCTTTAGCA‑GCCTAATTTCGACGTCCAACTGATCAAACAGTTAGGCGATTAAATCGTCAGCGGCGGCGGGAAAAGAATTCGTCCGCTGATTTCTGTACt                                      >  1:305977/1‑149 (MQ=255)                                                                                                                   aaaTCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTAAAAGGAAAAACTTTAGCA‑GCTTAAATTCGACGGCCAAAAGAAAAATCAATTAGGGTATTACATCGTCAGCGGGGGGGGGGAAAATATTCGTCCGGAGGTTGTTGTATTgg                                    >  2:316413/1‑150 (MQ=255)                                                                                                                     aTCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATTCGGCAATCCCTGGGCA‑GCTTAATTTCCACGTCCAACTGATCAATCAGTTAGGGTATTAAATCGTCAGCGGCGGCGGGAAAAGAATTCGTCCGCTGATTTCTGTACTgg                                    >  1:282902/1‑148 (MQ=255)                                                                                                                        aaTGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCAc                              <  1:623522/151‑1 (MQ=255)                                                                                                                         aTGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACg                             <  2:772203/151‑1 (MQ=255)                                                                                                                                 aCCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTg                     <  1:845256/151‑1 (MQ=255)                                                                                                                                       cAAGATATGGCGGGGGTTAATGCGGGGAACCCTGAGCA‑GCTCAATTTCGACGTCCAACTGATCAATCAGTTTTGGTATTTTATTGTCAGCGGGGGCGGGGGAAGTATTTGTGCGAAGATTGCTGTTCTGGGTGCGCGAGCTGTTggggt                 >  2:265866/1‑146 (MQ=255)                                                                                                                                        aaGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGa              <  1:810253/151‑1 (MQ=255)                                                                                                                                         aGATATGGCGGGGGTTTATGCGGGAATCCTTTAGCA‑GCTTAATTTCGACGTCCAACTGATCAAACAGGTAGGCTATTACATCGTCAGCGGCGGCGGGAAAAGTATTCGTCCGCTGATTGTTGTACTGGCTGCACCAGCTGTTgg                    >  2:890559/1‑144 (MQ=255)                                                                                                                                                   ggTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATgcgc   <  2:614700/151‑1 (MQ=255)                                                                                                                                                   ggTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATgcgc   >  1:65242/1‑151 (MQ=255)                                                                                                                                                   ggTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATgcgc   >  2:126486/1‑151 (MQ=255)                                                                                                                                                    gtgtTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCa  <  1:953595/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GTGCCATTTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCA‑GCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCA  >  NC_000913/3333609‑3333906

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