Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,332,791	0	A	C	41.4%	102.7 / 14.0	58	V161G (GTG→GGG)	basS	sensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/24);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.79e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332643‑4332941                                                                                                                                                     |                                                                                                                                                       aTGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcc                                                                                                                                                      >  1:865551/1‑151 (MQ=255)   gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcccc                                                                                                                                                    <  2:308404/151‑1 (MQ=255)   gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAccc                                                                                                                                                     <  2:818621/150‑1 (MQ=255)     tGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCg                                                                                                                                                  >  1:941929/1‑151 (MQ=255)        tAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGcc                                                                                                                                                >  2:351865/1‑150 (MQ=255)           ttACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCAAGCAGTTCCAGATGCAAACGCACCCCCGCCa                                                                                                                                               <  1:422748/148‑1 (MQ=255)            tACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTgg                                                                                                                                           >  2:124270/1‑151 (MQ=255)                ggAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGtt                                                                                                                                       >  2:172077/1‑151 (MQ=255)                    gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                    >  1:480162/1‑150 (MQ=255)                       aaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGtt                                                                                                                                >  2:341632/1‑151 (MQ=255)                          gACTGCCGGGCCCGCGCCATTGGCAGCGGTTGGGAGCCCCTCTCCATCATCTGATAAAGCCGTGCCACTAACGGTGCTACATCAATGTGTTGCTTTTTCGCCAGCAGTTCAAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                              <  1:474518/150‑1 (MQ=255)                                ccGGCCCGCGCCAGTTGCAGCAGCGGGGAGCCCCTCCCCATCATTTAACCACCCCGTCCCCCAACCGGTGCTACATCAATGTGGTGCTTTTCCGCCAGCATTCCCGGAGCCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                         <  2:591109/149‑1 (MQ=255)                                   gCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACg                                                                                                                     >  1:433004/1‑150 (MQ=255)                                       gcgcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGc                                                                                                                >  1:453184/1‑151 (MQ=255)                                        cgcCAGTTGCAGCAGCTGGGAGACGCTCCCCATCATCTGAACAAGCCGTGCCACTAGCGGTCCTACATCAATGTGATGCTTTTCCGCCAGCGTTCCCAGTTGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCa                                                                                                                  <  1:477854/148‑1 (MQ=255)                                          caggtagcagggggggggAGCCCCTCCCCCTCTTCTAACCAACCCGTGCCACTAGCGGTACTACATCATGGTGGGGCTTTTCCCCCAGCAGTTCCAGATGAAACCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                             <  2:623545/136‑1 (MQ=255)                                            aGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCgg                                                                                                              >  1:260898/1‑148 (MQ=255)                                             gtggcagcaGCGGGGAGCCCCCCTCCATCTTCTGATCAAGCGGTCCCACTAACGGTGCTACATCATGGTGATGCTTTTCCGCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaa                                                                                                           <  2:176215/147‑1 (MQ=255)                                              tcgcggcagcTGGGAGACGCTCCCAATCATCTGACCAAGCCGTCCCACAACCGGTACTACATCAATGTGATGCTTTTCCGCCGGCAGTCCCAGACGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTa                                                                                                            <  2:407600/143‑1 (MQ=255)                                                 cagcagCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaacaac                                                                                                      >  2:229848/1‑151 (MQ=255)                                                   ggaggtGGGACCCGCTCCCCGCCTTCTGAGCACGCCGGCCCAATAAGGTTTCTACACAAATGTGAGGTTTTTTCCCGAGCTCTTCCAAATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCt                                                                                                    <  2:217439/146‑1 (MQ=255)                                                         gggggcccctctCCATCTTCTGACAAACCCGTCCCACTAACGGTGCTACATAAATGTGTTGCTTTTCCGCCGGCAGTTCCAGAGGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCtt                                                                                                   <  1:1053596/138‑1 (MQ=255)                                                         gggAGACGCTCTCCTCTTTCTGATCAACCCGTGCCACAACCGGTTCAACACCAATGTGATGGTTTTCCGCCGGCAGTCCGGGAGGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                        <  2:142951/141‑1 (MQ=255)                                                          ggggaCGCTCCCCACCATTTAATCACGCGGTGCCATAAAGGGTGCAACACCAATGTGATTTGTTTCCGCCAGCGTTCCCGGAGGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                                  <  1:634138/143‑1 (MQ=255)                                                               ccctctCCATCATCTAACAAACCCGTGCCACAAACGGTGCTACATAAATGTGATGCTTTTCCCCCAGCAGTTCCAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                        <  1:983501/149‑1 (MQ=255)                                                               cGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATcc                                                                                          >  1:855316/1‑149 (MQ=255)                                                                 ctccCCATCATCTGATCAACCCGTCCCACTAACGGTGCTACATAATTGTGATGCTTTTTCGCCAGCAGTCCCAGATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                        <  2:456344/149‑1 (MQ=255)                                                                     ccATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCt                                                                                  >  1:818171/1‑151 (MQ=255)                                                                       atcatcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTgg                                                                                >  1:302826/1‑151 (MQ=255)                                                                              gacaaggCCGTCCCACTACGGGTGCTACATCAATGTGATGCTTTTTCGCCAGCATTTCCAGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCa                                                                             <  1:131413/141‑1 (MQ=255)                                                                               aTCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGCGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                         >  1:366508/1‑150 (MQ=255)                                                                                     ccGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                  >  2:836797/1‑151 (MQ=255)                                                                                          ccACTACCGGTGCTACAACAATGTGATGCGTTTTCCCCAGCAGTTCAAGATGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCTACTGACCAGAt                                                              <  1:17472/150‑1 (MQ=255)                                                                                             cTAACGGTGCTACACCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGtt                                                          >  1:189734/1‑151 (MQ=255)                                                                                               aagggTTCTCCATAAATGTGTTGGTTTTCCCCCAGCAGTCCCAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                             <  1:75887/143‑1 (MQ=255)                                                                                               aaCGGTTCTACATCAATGTGTTGCGTTTTCGCCAGCAGTCCCAAATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAg                                                                <  1:1023890/143‑1 (MQ=255)                                                                                                  gtggttacagaattgtgATGCTTTTCCCCCGCCAGTCCAAGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTa                                                        <  1:495600/135‑1 (MQ=255)                                                                                                  ggTGCTACATAAGTGTGTTGCTTTTTCGCCGGCAGTTCCAGATGAAACCCCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                             <  1:42603/143‑1 (MQ=255)                                                                                                    tGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                                   >  2:280152/1‑151 (MQ=255)                                                                                                    tGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                                   >  2:925493/1‑151 (MQ=255)                                                                                                     gCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                                  >  1:407082/1‑151 (MQ=255)                                                                                                      cTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                                  >  1:112806/1‑150 (MQ=255)                                                                                                         cataattgtgATGGTTTTCCGCCGGCCTTTCGAGATACAACCACCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                                    <  2:652111/139‑1 (MQ=255)                                                                                                             attggagtgtgTTTTCGCCATCAGTCCCATATCAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                                 <  2:367755/134‑1 (MQ=255)                                                                                                             aaTGTGATGCTTTTTCGCCAGCAGTTCCAGATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                          <  1:977762/151‑1 (MQ=255)                                                                                                             aaTGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                          >  2:438313/1‑151 (MQ=255)                                                                                                             aaTGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                          >  1:988896/1‑151 (MQ=255)                                                                                                               tgtgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGc                                        >  2:152702/1‑151 (MQ=255)                                                                                                                       tttttCGCCAGCATTCCCAGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAt                                  <  1:651052/148‑1 (MQ=255)                                                                                                                              ccAGCAGTTCAAGATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                         <  1:901476/151‑1 (MQ=255)                                                                                                                                    gtccctaaTGCAACCGCCCCACCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAATCCACCGCTTCGATTTCGAGGGTGgc                      <  1:902832/141‑1 (MQ=255)                                                                                                                                      tCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAT‑CAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                >  1:241575/1‑151 (MQ=255)                                                                                                                                       ccggagccaaccccACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg                  <  2:343552/136‑1 (MQ=255)                                                                                                                                       ccAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                >  2:843395/1‑151 (MQ=255)                                                                                                                                            tGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg            <  1:776447/150‑1 (MQ=255)                                                                                                                                             gCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATgg           >  2:91643/1‑150 (MQ=255)                                                                                                                                               aaaCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGaa              <  1:336748/145‑1 (MQ=255)                                                                                                                                                    cACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGc   >  1:770774/1‑151 (MQ=255)                                                                                                                                                     aCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCg  >  1:357566/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       ATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332643‑4332941

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