Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,511,772	0	A	C	40.0%	113.5 / 23.5	72	V312G (GTG→GGG)	mntH	manganese/divalent cation transporter
Reads supporting (aligned to +/- strand):  ref base A (38/4);  new base C (0/28);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.29e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGT  >  NC_000913/2511624‑2511922                                                                                                                                                       |                                                                                                                                                       agCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGGCGGATCTAATCCCACCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGTCCTGCATcaccac                                                                                                                                                       <  2:205763/150‑1 (MQ=255)     aTACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGcc                                                                                                                                                  >  1:1414349/1‑151 (MQ=255)      tACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGccc                                                                                                                                                 >  1:1017141/1‑151 (MQ=255)           aaaCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGcc                                                                                                                                             >  2:580634/1‑150 (MQ=255)           aaaCTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAAGGACAAAAAAGGACGGCAACATGGTGACTGTCCGACGCACCCACAGCGGAATATA‑AAACCG‑AATGAAGCCCGGCATCACCACCGGCCCCGCCa                                                                                                                                            <  2:229970/151‑1 (MQ=255)              cTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAgtg                                                                                                                                         >  1:1196831/1‑151 (MQ=255)                 aaCAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTcc                                                                                                                                      >  1:486556/1‑151 (MQ=255)                                aTAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGa                                                                                                                         >  1:1041794/1‑149 (MQ=255)                                aTAACAAGAACCCGTGGCGGATCAACCCCCACCAGAAACACAAAAAATGACGGAAACATGGTGCCTGTACGACCCCCCCCCGGGGGAATAGA‑AAGGCG‑AATGAGCCCCTGCTCCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAg                                                                                                                       <  2:12947/151‑1 (MQ=255)                                  acccAAAATCCGTGTCGGATCAAACCCCATCAGAATAAAAAAAAAGGACGGAAACAGGGTGACTGTACGACGCCCCCACAGGGGGATATA‑AAAGCA‑AATGAAGCCCTCAATCACCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGa                                                                                                                      <  1:1383400/148‑1 (MQ=255)                                  aacaaaaacccGTGCAGGATCTAATCCCAACAAAATCACAAAAAATGGCGGAAATTTGGGAACTCACCAACCCGCCCAAAGGGGAATGGG‑AAAGGA‑AATGAGCCCCGGCACAACCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAg                                                                                                                       <  1:192531/140‑1 (MQ=255)                                   aCCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGAc                                                                                                                     >  1:1027206/1‑150 (MQ=255)                                    ccAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGAc                                                                                                                     >  1:148671/1‑149 (MQ=255)                                    ccAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGAc                                                                                                                     >  2:759631/1‑149 (MQ=255)                                      aGATTCCGTGTGGGACCTAACCCCATCAGAAGAAAAAAAAAGGGCGACAACAGGGTGACGGTCCGACCCCCCCACCGCGGGATATG‑AAGGCG‑AATGACGCCCTCCAACCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACa                                                                                                                    <  2:343538/148‑1 (MQ=255)                                             gtgtCGGATCTAATCCCATCAGAATGACAATACATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATAACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCgcagc                                                                                                          >  1:468505/1‑151 (MQ=255)                                               gtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCgcagca                                                                                                         >  1:982590/1‑150 (MQ=255)                                               gtCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCgcagca                                                                                                         >  1:519821/1‑150 (MQ=255)                                                 cGGTTCTAATCCCAACAAAATGAAAAAAAATGACGGCAACAGGGTGACTCTCCGCCGCCCCCACGGCGGTAAATA‑AAGGCG‑AAGGAAGCCCTCAACCCCCACCGGCACCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCg                                                                                                              <  1:399107/143‑1 (MQ=255)                                                  ggttcaaagccaagcAGAATGACAAAAAATGAGGGCAACAATGTGACTGCCCGACGCCCCCACAGCGGGATATG‑AAAGTA‑AAGGAAGCCCGGCACCGCCCCCCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAAcc                                                                                                      <  1:205165/136‑1 (MQ=255)                                                  ggATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAAcc                                                                                                      >  2:905710/1‑150 (MQ=255)                                                  ggATCTAACCCCACCAGAATAACAAAAAATGCCGGCAACATGGTGACTGTACGCCGCCCCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCa                                                                                                     <  1:1302454/151‑1 (MQ=255)                                                    aTCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAgg                                                                                                   >  2:477896/1‑151 (MQ=255)                                                     tCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGc                                                                                                  >  1:277390/1‑151 (MQ=255)                                                      cTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCt                                                                                                 >  1:215084/1‑151 (MQ=255)                                                       tAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCtt                                                                                                >  1:1273256/1‑151 (MQ=255)                                                           cccaaaaGAATAAAAAAAAATGCCGCAAACAGGGTGACTGTCCGCCCCCCCCACGGCGGAATATG‑AAAGCA‑AAGGAACCCCTCCACCCCCCCCGGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATc                                                                                            <  2:191699/145‑1 (MQ=255)                                                              aTCAAAATGACAAAAAAGGACGGCAACATGGTGACTGTACGACGCCCCCACAGGGGGATATG‑AAAGCA‑AATGAAGCCCTGCATCCCCCCCGGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCa                                                                                          <  1:109798/150‑1 (MQ=255)                                                                cAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCaa                                                                                         >  1:1093669/1‑149 (MQ=255)                                                                 agaaagaaAAAAAATGACACCAACTTGTTGCCTGTACGACCCCCCCCCAGCGGGA‑ATGAAAAGCA‑AAGGAACCCCTGCACCACCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCa                                                                                          <  2:369392/147‑1 (MQ=255)                                                                 aGAATGAAAAAAAATGACGGCAAAGTGGTGACTCGACGACGCCCCCACACCGGAATA‑G‑AAAGGGAAATGAACCCCTGCACCCCCGCCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCaaa                                                                                        <  1:274316/149‑1 (MQ=255)                                                                       aCAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTg                                                                                >  2:1414334/1‑151 (MQ=255)                                                                       aCAAAAAATGGCGGCAACAGGGTGACTGTGCGCCTCCCCCACAGCGGGATAGG‑AAAGCG‑AATAAGGCCCTGCATCTCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGtt                                                                                 <  1:354409/150‑1 (MQ=255)                                                                         aaaaaaGGAGGGAAACATGGTGACTGTCCGACCCCCCCACGGCGGGATATA‑AAAGCG‑AATGAAGCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGc                                                                               <  1:813583/147‑1 (MQ=255)                                                                         aaTAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGcc                                                                              >  2:1145506/1‑151 (MQ=255)                                                                         aaTAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGc                                                                               >  2:1346225/1‑150 (MQ=255)                                                                          aaaaaTGACGAAAACGTGGTGACTGTACGCCGCACCCCCAGCGGAATAAA‑AAAGCG‑AATGAAGCCCTGCACCCCCCCCCCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCg                                                                             <  2:102196/149‑1 (MQ=255)                                                                            aaatcacgaaatcagggTGATTGTACGACGCACCCACAGGGGTAAACA‑AAGGCG‑AATAAAGCCCGGCACCCCCCCCCGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCa                                                                           <  1:391562/136‑1 (MQ=255)                                                                               tGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGc                                                                         >  2:143868/1‑150 (MQ=255)                                                                                         aTGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAAc                                                               >  2:1036420/1‑150 (MQ=255)                                                                                           gtggACTGACCGACGCACCCACAGCGGGATATG‑AAAGCG‑AAGGAACCCCTGCCCCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAg                                                             <  2:208097/147‑1 (MQ=255)                                                                                           ggTGACTGTACGACGCCCCCACAGCGGAATATA‑AAAGCG‑AATGAACCCCTGCACCACCACCGCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAg                                                             <  1:237848/150‑1 (MQ=255)                                                                                               aCTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGtt                                                        >  2:726630/1‑151 (MQ=255)                                                                                                  gtccgccccccccACAGCGGAATATA‑AAGGCA‑AAGGACCCCCTCAACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCgg                                                          <  1:407136/136‑1 (MQ=255)                                                                                                  gTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCACCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTCGc                                                      >  2:30549/1‑150 (MQ=255)                                                                                                        cgccccccaagcgcaatatA‑AAAGCG‑AATAATGCCGTGAACCCCCCCCGGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCa                                               <  1:494970/136‑1 (MQ=255)                                                                                                        ccccccccCAGGGGGAAAAG‑AAAGCA‑AATGAGGCCCTGCACCACCCCCGGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAg                                                    <  2:533744/142‑1 (MQ=255)                                                                                                         gCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTc                                                >  1:798350/1‑149 (MQ=255)                                                                                                         gCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCC‑CCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGCCCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTc                                                >  2:500792/1‑148 (MQ=255)                                                                                                            cccACGGCGGCATAAA‑AAACCG‑AAGGACCCCCACCATCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGAt                                            <  2:601807/150‑1 (MQ=255)                                                                                                            cccACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGa                                             >  2:1067697/1‑149 (MQ=255)                                                                                                            cccACAGCGGGATAAG‑AAGGCG‑AATAAGCCCCGGCATCCCCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCa                                               <  1:1207148/147‑1 (MQ=255)                                                                                                                    gggATATG‑AAAGCA‑AATGACGCCCTGCATCACCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtca                                    <  1:607659/150‑1 (MQ=255)                                                                                                                    ggaatatG‑AAAGCG‑AATAAACCCCTGCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtc                                     <  1:117182/146‑1 (MQ=255)                                                                                                                    gagATAAG‑AAAGCG‑AATGAACCCCGCCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtca                                    <  2:486549/148‑1 (MQ=255)                                                                                                                     ggATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtcatc                                  >  1:441713/1‑151 (MQ=255)                                                                                                                         atG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtcatca                                 >  2:544534/1‑148 (MQ=255)                                                                                                                         atG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtcatca                                 >  2:499680/1‑148 (MQ=255)                                                                                                                             aaaGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCtcatca                                 >  1:24827/1‑145 (MQ=255)                                                                                                                              aacgccaaccaaGCCCGGCAGCACCGCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCaa                                <  2:496093/136‑1 (MQ=255)                                                                                                                                 gcaaatgacccctggCATCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCgg                          <  1:167636/135‑1 (MQ=255)                                                                                                                                  g‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAac                      >  2:1142469/1‑151 (MQ=255)                                                                                                                                    aaggAACCCCTGCACCCCCCCCGGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCa                        <  1:494199/145‑1 (MQ=255)                                                                                                                                    aaTGAAGCCCTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAac                      <  1:175208/150‑1 (MQ=255)                                                                                                                                     agaagccccTCCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAacac                    <  1:1414411/145‑1 (MQ=255)                                                                                                                                        aaGCCCTCCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCa                  <  2:1224689/150‑1 (MQ=255)                                                                                                                                          gCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGt                >  1:1449474/1‑150 (MQ=255)                                                                                                                                           cccTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTAt              >  2:1075641/1‑151 (MQ=255)                                                                                                                                           cccTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTAt              >  1:363765/1‑151 (MQ=255)                                                                                                                                           cccTGCACCCCCCCCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTAt              <  1:726646/151‑1 (MQ=255)                                                                                                                                            ccTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATg             >  1:17934/1‑151 (MQ=255)                                                                                                                                            ccTGCACCCCCCCCTCCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAacac                    <  2:548495/144‑1 (MQ=255)                                                                                                                                               gCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCATGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGAc           >  1:561888/1‑150 (MQ=255)                                                                                                                                                       accTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGt  >  1:970194/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AGCGATACCAAAACTTAACAGTACCTGACTCATAACCAGAATCCGTGTCGGATCTAATCCCATCAGAATGACAATAAATGACGGCAACATGGTGACTGTACGACGCACCCACAGCGGGATATG‑AAAGCG‑AATGAAGCCCTGCATCACCACCTGCCCCGCCAGTGTCCCCACCACCGTTGAAGACAGCCCCGCAGCAACCAGGCTTAATCCAAAGACCGTTGCCGCAGCGTGGCTTAACAGCGGTTGCAGCGTCAGATAAGCCTCATCAAGATCGGCAACACCAGTATGACCGGAAAAGT  >  NC_000913/2511624‑2511922

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: ‑