Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,850,007	0	G	A	45.8%	53.6 / 39.3	59	P44S (CCG→TCG)	uhpB	sensory histidine kinase in two‑component regulatory sytem with UhpA
Reads supporting (aligned to +/- strand):  ref base G (26/6);  new base A (0/27);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 2.32e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACATGCG  >  NC_000913/3849859‑3850156                                                                                                                                                       |                                                                                                                                                      gTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGa                                                                                                                                                      >  2:1043799/1‑150 (MQ=255)       aacctcccGATCATCATAAACGGAAAATGGGTTAACCCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACGCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAAc                                                                                                                                                <  2:135472/144‑1 (MQ=255)       aaccTACCGATCATCAATAACGGAAAATGGGTTAACCCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATaa                                                                                                                                                 <  1:268109/146‑1 (MQ=255)        aaCTACCGATCATCAATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAAc                                                                                                                                                >  1:127096/1‑149 (MQ=255)         aCTACCGATCATCAATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGc                                                                                                                                             >  1:229934/1‑151 (MQ=255)            aCCGATCATCAATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCa                                                                                                                                            <  1:93294/149‑1 (MQ=255)             ccGATCATCAATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCg                                                                                                                                         <  1:356871/151‑1 (MQ=255)                 tcatcaATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCat                                                                                                                                     >  2:308720/1‑151 (MQ=255)                 tcatcaATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCat                                                                                                                                     >  1:1158382/1‑151 (MQ=255)                 tcatcaATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCat                                                                                                                                     >  2:663922/1‑151 (MQ=255)                     caATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATc                                                                                                                                  >  2:558550/1‑150 (MQ=255)                       aTAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAg                                                                                                                                >  2:832787/1‑150 (MQ=255)                            ggAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCt                                                                                                                          >  2:714581/1‑151 (MQ=255)                               aaaTGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCa                                                                                                                        >  2:238168/1‑150 (MQ=255)                               aaaGGGGTTAAACCGCCCCGCCGGCGTGAGCAAAAAAACCACCCCCCCCCCCGCGCCAAGCAAGGGGGC‑CAAGTCCCCGCGCGGGCATTGCAGCATTACCCCGAGACGCAGACAAAAGGAAAAAAAAAGCACCGCCATATCAGGgcgc                                                                                                                           <  2:317064/148‑1 (MQ=255)                                aaTGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAAc                                                                                                                      >  1:352769/1‑151 (MQ=255)                                    ggttaaccccccccgCCTGCGTTAGCAATAAAACCAGCGCCCCCCCCGCCCCCAGCAATGCGGG‑CCAGTACCCGCGGGGGCATGGCAGCATTCCCCCCAGACGCAGCCAAAACGGAAATAACGGCACCGCCATATCAGGGCGCTCAAcc                                                                                                                     <  2:1014360/137‑1 (MQ=255)                                       tAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                               >  2:4055/1‑151 (MQ=255)                                       tAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                               >  1:549142/1‑151 (MQ=255)                                       aaacccg‑ccccCCTGCGTTAGCCAGAAAACCACCACCCCCCCCGCCCCCAGAAATGCGGG‑CCAGTCCCCGCGGGGGCATGGCAGCATTACCCCCAGACGCAGACAAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGc                                                                                                               <  1:814318/147‑1 (MQ=255)                                            cg‑accgCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCt                                                                                                          >  1:154387/1‑151 (MQ=255)                                                       cGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCcac                                                                                                >  2:1168048/1‑151 (MQ=255)                                                                 aaaaCCAGCACCCCCCCCGCCCCCAGCAATACGGC‑CCAGTACCCGCGCGGGCATTGCACCATTAGCCCGAGCCGCAGACAAAAGGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAAcc                                                                                      <  2:332932/150‑1 (MQ=255)                                                                   aaTCAGCACCCCCCCCGCCCCCACGAATAGGGC‑CCAGTACCCGGGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCAt                                                                                    <  1:267808/151‑1 (MQ=255)                                                                     gccgcagcCCCCCCGCGCCCACAAGGCCGGA‑CCTGTTTCCGCGCGGGTATTGCAGCTTTCGCCCCAGCCGCAGTCCAAGAGGAAACACCGGAACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATg                                                                                   <  1:687847/147‑1 (MQ=255)                                                                     ccagccccccccccccgcCCAGCAACGCGGG‑CCATTACCCGCGGGGGCATGGCAGCATTACCCCCAGACGCAGCCAAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGc                                                                                  <  2:689698/136‑1 (MQ=255)                                                                       cgcccccccccccggcCCAAGAAAACGGGCCCGGTACCCGCGCGGCTATGACAGCATACCCCCCAGACGAAGCCAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATg                                                                                   <  1:632017/135‑1 (MQ=255)                                                                          agcCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                               >  2:514814/1‑149 (MQ=255)                                                                            cccccccccccccAGAAATACGGG‑CAATTCCCCGGGCGGGAATTGCACCATTACCCCAGGCCGAAAACAAAAGGAAAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                                                                               <  2:121209/137‑1 (MQ=255)                                                                                    cgcCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaa                                                                    >  1:194904/1‑150 (MQ=255)                                                                                        cagcacgccGGC‑CCAGTACCCGCGCGGGCATTGCAGCATTCCCCCCAGACGCAGACAAAACGAAAAAACCGGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaaa                                                                <  2:739070/142‑1 (MQ=255)                                                                                         aGCAATCCGGC‑CCAGAACCCGCGGGGGCATTGCACTATTACCCCCAGCCGCAGACAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATaaaaaa                                                                <  1:335392/149‑1 (MQ=255)                                                                                            aacaCGGG‑CCAGTACCCGCGCGGGCTTTGCAGCATTACCCCCAGACGCAGACCAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGc                                                              <  1:137687/145‑1 (MQ=255)                                                                                            aaTGCGGC‑CCAGTACCCGCGGGGCCATTGCAGCATTACCCCCAGACGAAGACAAAACGAAAATACCAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAgg                                                           <  1:528911/151‑1 (MQ=255)                                                                                                 ggg‑CAAGTCCCCGCGCGGGTATGGCAGTATTGCCCCCAGAGCCAGCCAAAGCGAAAAAACCAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATa                                                      <  1:1169618/151‑1 (MQ=255)                                                                                                     ccAGTACCCGCGCGGGAATGGCAGCATAACCCCCAGACGAAGACAAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAAc                                                    <  1:1062013/150‑1 (MQ=255)                                                                                                      cAGTACCCGCGGGGGCATTGCAGCATTACCCCCAGACGCAGACAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAAc                                                    <  1:73807/149‑1 (MQ=255)                                                                                                       aGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGt                                                 >  2:294444/1‑151 (MQ=255)                                                                                                         tctCCGGGGGGCCATCGCAGCATGACCCCCGGCCGCACACCAAACGAAAATAACAACCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTa                                                <  2:97146/148‑1 (MQ=255)                                                                                                                gcgGGCATTGCAGCATAACCCCCAGACACAGACAAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGc                                         <  1:557612/150‑1 (MQ=255)                                                                                                                  ggtcATGGCACCATAACCCCCAGGGGAAAACAAAGCGAAAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAg                                          <  2:106000/144‑1 (MQ=255)                                                                                                                  gggCATTGCAGAATCACCCCCAGACGAAGACAAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGc                                         <  1:562816/148‑1 (MQ=255)                                                                                                                  ggactttGCAGCTTTACCCCCAGCCGCAGACAAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAAt                                              <  1:1009569/138‑1 (MQ=255)                                                                                                                   ggCATGGCAGCATTACCCCCAGACGAAACCAAAACAAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGga                                     <  1:905485/151‑1 (MQ=255)                                                                                                                      aTTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgag                                    >  1:20862/1‑149 (MQ=255)                                                                                                                        tGCAGCATGACCCCCAGCCGCAGACAAAAGGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGgag                                    <  1:504122/147‑1 (MQ=255)                                                                                                                           agcaTTGCCCCCAGACGCAGACCAAACGAAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaac                              <  1:80063/150‑1 (MQ=255)                                                                                                                               tcacccccAGCCGCAAACAAAAGGAAAAAAACACCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaacaa                               <  2:202800/141‑1 (MQ=255)                                                                                                                               taacccccAGACGAAAACAAAAGGAAAAAAACAACACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTc                           <  2:140886/145‑1 (MQ=255)                                                                                                                               taacccccAGACCCAACCAAAAAGAAAATAACAGCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGaa                                  <  1:556614/138‑1 (MQ=255)                                                                                                                                tAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCtt                         >  1:664146/1‑150 (MQ=255)                                                                                                                                 aGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTtc                        >  1:793477/1‑150 (MQ=255)                                                                                                                                             aGACCAAACGAAAAAAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAGAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCAt            <  1:80305/150‑1 (MQ=255)                                                                                                                                             aGACAAAACGAAAAAAAGACCCCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCAt            <  2:983842/150‑1 (MQ=255)                                                                                                                                              gACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATc           >  2:444334/1‑150 (MQ=255)                                                                                                                                              gACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATc           >  1:441589/1‑150 (MQ=255)                                                                                                                                                 cAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAc       >  1:1135671/1‑151 (MQ=255)                                                                                                                                                  aaaCGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAAc       >  1:167170/1‑150 (MQ=255)                                                                                                                                                    aCGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACATg    >  2:1038178/1‑151 (MQ=255)                                                                                                                                                      ggAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACATgcg  >  2:636758/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCG‑ACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGG‑CCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATAACGGTAATTAAGCGGGAGAACAACGTCTTCATCACCAGCCATCAAACATGCG  >  NC_000913/3849859‑3850156

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