Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,056,115	0	A	C	40.0%	106.2 / 15.6	71	G75G (GGT→GGG)	glnL	sensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base A (31/11);  new base C (0/28);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 6.45e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTT  >  NC_000913/4055966‑4056264                                                                                                                                                       |                                                                                                                                                       gCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAAc                                                                                                                                                       <  2:521666/151‑1 (MQ=255)        gTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑Gc                                                                                                                                                 <  2:739523/149‑1 (MQ=255)         tAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑Gcccc                                                                                                                                              >  2:1221466/1‑151 (MQ=255)          aGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑Gcccc                                                                                                                                              >  2:809337/1‑150 (MQ=255)          aGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑Gcccc                                                                                                                                              <  1:193825/150‑1 (MQ=255)           gCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGc                                                                                                                                            >  1:98130/1‑151 (MQ=255)           gCTGTTCCTGACTTAAGCGGCGCTGGATATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGCTGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑Gccc                                                                                                                                               >  1:1028940/1‑148 (MQ=255)              gTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCt                                                                                                                                          <  2:542186/150‑1 (MQ=255)                 ccTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTcc                                                                                                                                        >  2:1193204/1‑149 (MQ=255)                   tGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGACGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGAc                                                                                                                                    >  2:789474/1‑151 (MQ=255)                    gtcTTAAGCGGCGCTGTTTATCCATCGGAGCCATCTCCAGCAGGATCAGGCCGTCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCTTCACTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACt                                                                                                                                   <  2:1116620/149‑1 (MQ=255)                         aaGCGGCGCTGGTTACCCATCGGAGCCATCCCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCAAAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTACCG‑GAAAACCCTT‑GCCCCGCCTCCAGACTTTc                                                                                                                                <  1:436371/149‑1 (MQ=255)                             ggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCCCCCGTCGATGACCAGCGTCCCTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCa                                                                                                                           <  1:508398/150‑1 (MQ=255)                                   ggttacccacggggcccATCTCAAGCGGGATCGTGCCGCCCGGCATCCGCTGGGCCGCCACAAAAAGGATATGCGAGCCCCCGTCGATGCCCAGCTTCCTTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCt                                                                                                                     <  2:314956/135‑1 (MQ=255)                                     ttACCCATGGGACCCACCTCCAGCAGCACCATCCCTTCCGGCGTACGCCGGGCCGTCACAGAAAGGATAGGCGACCCCCCTTCTATGCCCACCGTCCCTTC‑TTTTCGTGTAACCCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGc                                                                                                                      <  2:411375/147‑1 (MQ=255)                                         ccATGGGAGCCACCTCCGGCGGGACCATGCCGTCCGCAATACGCTGGGCCGTAACAAAAAGGAAATGCGGGCGCCCTTCGATGACCAGCGTCATTTGGTTATCG‑TAAAACCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAtt                                                                                                              <  2:725425/151‑1 (MQ=255)                                         ccATCGGGGCCATCCCCGGCAGGATCAGGCCGCCCGGCATGCGCGGGGCCGTCACAAAAGGAATATGCGACCCCCCGTCGATGGCCAGCGTCACTTCGTTACCG‑GAAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGa                                                                                                                  <  1:211991/147‑1 (MQ=255)                                         ccAGCGGAGCCATCCCCAGCAGGTTCATGCCGTCCGGCATACGCTGGGCCGTCACAAAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGata                                                                                                                <  1:1050804/149‑1 (MQ=255)                                          cATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                                                                                             >  2:799253/1‑151 (MQ=255)                                               tggCCATCTCCAGCAGGATCTTCCCGTCCGTCCTACGTGGGGCCTTAACAGAAAGGATATGCGACCCCCCGTCGATGACCAGCGTCCCTGCTTTATCG‑TAAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTa                                                                                                            <  2:27288/145‑1 (MQ=255)                                                aGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                         >  1:560735/1‑149 (MQ=255)                                                aGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                         >  1:411967/1‑149 (MQ=255)                                                  ccATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCAGAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCCCTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGaa                                                                                                     <  1:743626/151‑1 (MQ=255)                                                  ccATCCCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCCCTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATga                                                                                                        <  1:245742/148‑1 (MQ=255)                                                   cATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATgaga                                                                                                      >  1:1039924/1‑149 (MQ=255)                                                      ccccAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCACAAAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                                                                                                 <  1:575811/149‑1 (MQ=255)                                                       tcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCATTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                                                                                                 >  2:1215361/1‑150 (MQ=255)                                                          agcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                                                                                              >  2:706057/1‑150 (MQ=255)                                                          agcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                                                                                              >  2:439138/1‑150 (MQ=255)                                                           gcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCTTCACTTCGTTATCG‑GAAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCa                                                                                             <  1:1240059/150‑1 (MQ=255)                                                              ggATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAAc                                                                                           >  1:1005151/1‑149 (MQ=255)                                                                  cATGCCGTCCGGCATACCCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGtt                                                                                       >  2:224188/1‑149 (MQ=255)                                                                   aTGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCg                                                                                    >  1:32740/1‑151 (MQ=255)                                                                     gCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATTTTTAATGAGAAGTAGCTCAACAGTTCCgg                                                                                   >  1:710123/1‑150 (MQ=255)                                                                       cGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTa                                                                                 >  2:643963/1‑150 (MQ=255)                                                                          cgggCATACGCGGGGCCGTCACAGAAAGGATATGCGCGCCCCCGTCGATGACCGTCGTCACTTCGTAATCG‑GAAAACCCTG‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGt                                                                                  <  2:427687/144‑1 (MQ=255)                                                                          ccGGCATCCGCGGGGCCGTCAAAAAAAGGATATGCGAGCGCCCGTCGATGACCGGCTTCATTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAAc                                                                               <  2:1064428/149‑1 (MQ=255)                                                                          ccGGCATACGCTGGGCCGTCACAAAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCTCTTCGTTATCG‑GAAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTaa                                                                                <  1:1084967/148‑1 (MQ=255)                                                                          ccGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCATTTCTTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTaa                                                                                <  2:58564/148‑1 (MQ=255)                                                                                  cGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCa                                                                      >  1:356513/1‑150 (MQ=255)                                                                                  cGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCa                                                                      >  1:1275519/1‑150 (MQ=255)                                                                                  cGCTGGGCCGGCAAAAAAAGGATATGGGCGCCCCCGGCGATGCCCAGCTTCTCTTCTTTATGG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCa                                                                      <  2:71873/150‑1 (MQ=255)                                                                                   gCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCa                                                                      >  1:977643/1‑149 (MQ=255)                                                                                       ggCCGCCACAAAAAGGTTATCCGCCCCCCCGTCGCCGCCCAGCTTCATTCCGTTAGCG‑GGAAATCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAAc                                                                   <  2:443279/148‑1 (MQ=255)                                                                                         ccGTCACAGAAAGGATATGCGAGCGCCGGTCGATGCCCAGCTTCATTTCGTTATCG‑GTAAACCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAAttt                                                              <  1:943217/151‑1 (MQ=255)                                                                                         ccGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAAttt                                                              >  1:659825/1‑151 (MQ=255)                                                                                            tCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCg                                                           >  2:683376/1‑151 (MQ=255)                                                                                                  aaaGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGa                                                        <  1:441688/148‑1 (MQ=255)                                                                                                  aaaGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCt                                                     >  2:77747/1‑151 (MQ=255)                                                                                                   aaGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCtt                                                    >  2:522066/1‑151 (MQ=255)                                                                                                       atagGCGAGCGCCCGTCGATGACCAGCGTCACTTCTTTATCG‑GAAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTgg                                                 <  2:1211657/150‑1 (MQ=255)                                                                                                        tatGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGc                                               >  1:1112283/1‑151 (MQ=255)                                                                                                          tGCGAGCCCCCCTGGAGGACAAGCTTCACTCCGTTGTCG‑GTAAACCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTgg                                                 <  1:659679/147‑1 (MQ=255)                                                                                                           gCGAGCGCCCGTCGATAACCACGGTCACTTCGTTATCG‑GAAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGa                                             <  1:585996/150‑1 (MQ=255)                                                                                                              aGCGCCCGTCGATGCCCAGCGTCACTTCTTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGa                                             <  1:325678/147‑1 (MQ=255)                                                                                                                  ccggtggATGCCCGGCGTCCTTTCGTTATCG‑GAAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgt                                        <  2:189315/142‑1 (MQ=255)                                                                                                                  cccGTCGATGACCAGCGTCCCTTCGTTATCG‑GTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttg                                      <  1:759348/150‑1 (MQ=255)                                                                                                                  cccGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgt                                     >  1:408951/1‑151 (MQ=255)                                                                                                                  cccGCCGATGACCGGCGTCATTTCGTTATGG‑GAAAACCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTggg                                                <  2:1078597/140‑1 (MQ=255)                                                                                                                   ccttCGATGACCCGCGTCCCTTCTTTATCG‑GAAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgtt                                    <  2:503256/148‑1 (MQ=255)                                                                                                                   ccgccgATGCCCAGCGTCTTTTCTTTATCG‑GAAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgtt                                    <  2:142924/151‑1 (MQ=255)                                                                                                                    cGTCGATGACCGGCGTCATTTCGTTATCG‑GTAAACCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGttgttg                                   <  1:1193262/151‑1 (MQ=255)                                                                                                                       cGATGACCAGCGTCACTTCGTTCTCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTgcg                                 >  2:104874/1‑150 (MQ=255)                                                                                                                             ccggcttctcttcttTATC‑‑GTAAAACCTTGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAg                           <  1:985854/141‑1 (MQ=255)                                                                                                                                 cGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTg                      >  1:1105055/1‑151 (MQ=255)                                                                                                                                      cctCGTTTTCG‑GTAAAACTTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTa                 <  2:425891/149‑1 (MQ=255)                                                                                                                                      cTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGt                  >  1:1206946/1‑150 (MQ=255)                                                                                                                                        tCGTTATCG‑GAAAACCCTT‑CCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTaa                <  2:842499/150‑1 (MQ=255)                                                                                                                                            ttatggtaaaaaCTTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATg              <  1:223598/141‑1 (MQ=255)                                                                                                                                               cG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGc        >  2:602765/1‑151 (MQ=255)                                                                                                                                                g‑gTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCg         <  1:347735/149‑1 (MQ=255)                                                                                                                                                g‑gTAAACCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCg         <  1:158352/149‑1 (MQ=255)                                                                                                                                                    aaacCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCa      <  2:67891/149‑1 (MQ=255)                                                                                                                                                      aaCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGtt  >  2:1151591/1‑151 (MQ=255)                                                                                                                                                      aaCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGtt  >  2:978904/1‑151 (MQ=255)                                                                                                                                                       aCCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAg     <  1:809366/147‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       GCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCG‑GTAAAACCTT‑GCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGTT  >  NC_000913/4055966‑4056264

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