Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,748	0	C	A	39.4%	60.2 / 44.3	66	E55* (GAG→TAG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (28/12);  new base A (0/26);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 2.15e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTG  >  NC_000913/3177603‑3177896                                                                                                                                                  |                                                                                                                                                     tCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaaa                                                                                                                                                  <  2:1068312/150‑1 (MQ=255)       tctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTccc                                                                                                                                           >  1:1221114/1‑151 (MQ=255)         ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                         >  2:600299/1‑151 (MQ=255)          tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTccc                                                                                                                                           <  1:174631/148‑1 (MQ=255)           cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg                                                                                                                                       >  2:873033/1‑151 (MQ=255)           cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg                                                                                                                                       >  2:1041917/1‑151 (MQ=255)           cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                        <  2:265211/150‑1 (MQ=255)                 ccccGGCACCCACCTCCAGTAGCCAAGGGGTTTCGCTGGTGCCGTGCGCGGCAACCCGATTCTCCACAATCAGCACAACTCCACCCCGCCCTGGGTCAAAGGGTAGAAAGACGGGGGGGGGCCCGCGCTCAAAATTTTCCCGCCGCACCTc                                                                                                                                 <  1:930970/150‑1 (MQ=255)                     ggCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATg                                                                                                                              >  2:221198/1‑150 (MQ=255)                         aCCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTc                                                                                                                          >  1:856030/1‑150 (MQ=255)                          ccATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCAt                                                                                                                        >  1:1153527/1‑151 (MQ=255)                          ccATCCCCAGTAGCCAGGGGTTTTCGCTGGTGTCGTCCGCGGCAACCCGAATCTGCCCAATCAGCACAACTTCACCCCGCACGGGGTAAAGGGGTAGCAAGACGGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCAt                                                                                                                        <  1:744371/151‑1 (MQ=255)                          ccATCCCCAGTAGCAAGGGGTTTTCGCTGGTGTCGGACGCGGCACTCCGAATCTGCTCAACCACCACACCTTCACCCCGCACTGGCTAAAAGGGTAGCAAGCCTGCGGGGGCCCCCCCCTAAAAATTTTCCCGCCGCACCTCATGACTCa                                                                                                                         <  1:669199/150‑1 (MQ=255)                              ctcCAGTAGCAAGGGGTTTTCGCTGGTGTCGTCCGCGGCAATCCGAATCTCCTCAAACACCACAATTTCACCCCGCACTGGGTAAAAGGGTACCAAGACTGGGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTg                                                                                                                     <  1:1111208/150‑1 (MQ=255)                                    tAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                              >  2:867083/1‑151 (MQ=255)                                        cAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCCCAATCAGCACAACTTCATCACGCACTGGGTCAAAAGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                           >  1:1190641/1‑150 (MQ=255)                                         aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                           >  2:827431/1‑149 (MQ=255)                                                   gCTGGTGTCGTACGCGGCAACCCGAATCTGCACAACCACCACAACTTCACCACGCACGGGGTAAAGGGGTAGCAAGACTGGGGCGTCCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGAc                                                                                                 <  1:665873/149‑1 (MQ=255)                                                      ggTGTCGTACGCGGCAACCCGAATCTGCTCAACCAGCACAACTTCACCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt                                                                                            <  1:371883/151‑1 (MQ=255)                                                      ggTGTCGTACGCGGCAACCCAAATCTCCTCAACCACCACAACTTCACCACGCACGGGTTCAAGGGGTAGAAAGCCTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                <  1:753282/147‑1 (MQ=255)                                                         gtgggaCGCGGCAATCCAAATCTGCTCAATCAGCACAACTTCACCACGCACGGGGCAAAGGGGTAGCAAGCCTGCGGCGTGCCCGCGCTCAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                          <  1:854660/145‑1 (MQ=255)                                                           cgcgcgcgGCAACCGGATTGTGCTCCATTACCACATCTTCATCACTCCCGGGGTAAAAGGATAGAATGTCGCCGGTATGCCCGCGCCAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCta                                                                                         <  1:1082326/145‑1 (MQ=255)                                                            gTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                       >  1:361843/1‑150 (MQ=255)                                                                gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa                                                                                      >  2:148993/1‑147 (MQ=255)                                                                gcgGCAATCCGAATCTGCTCAACCAGCACAACTTCACCACGCGCTGGTGCAAAGGGTAGCAAGACTGCGGCGTGGCCTCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                                   <  1:36955/150‑1 (MQ=255)                                                                 cgGCAATCCGAATCTGCTCAACCAGCACAACTTCATCACGCGCTGGGTAAAAGGGTAGCAAGCCGGCGGCGTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGAt                                                                                  <  2:461477/150‑1 (MQ=255)                                                                  ggCAATCCGAATCTCCTAAACCACCACAACTCCACCACCCCCTGGGTAAAAGGGTAGAAAGCCTGGGGGGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                                 <  1:115991/150‑1 (MQ=255)                                                                  gaaaaccaaaacttGCTCAACCAGCACATCTTCACCACGAATGGGGTCAAAGGGAAGCAAGCCTGGGGCGTCCCCGCGCACAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                                 <  1:647138/137‑1 (MQ=255)                                                                          gAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaa                                                                         <  2:1030515/150‑1 (MQ=255)                                                                           aaTCTGCTCAATCAGCACACCTCCACCAGGCACTGGGTCAAGGGGTGACAAGACGGGGGCGTAACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                       <  2:643906/151‑1 (MQ=255)                                                                           aaTCTGCTCAATCAGCACAACTTCATCACGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                       <  2:1047167/151‑1 (MQ=255)                                                                           aaTCTGCCCAACCAGCCAAACTACACCACGCACTGGGTAAAAGTTTAGCAAGCCGGGGGCGTGCCCGCTCTCAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                                           <  2:1182723/147‑1 (MQ=255)                                                                            aTCTACTCAACCACCACAACTTCCCCCCGCACGGGTTAAAAGGGTAGAAAAACTGCGGCGTGGCCGCGCAAAAAAATTTGCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                       <  1:967711/150‑1 (MQ=255)                                                                             tCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                      >  2:885978/1‑150 (MQ=255)                                                                                gCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCAACTCATGACTCATTTGCCCGTTGAATAGACAATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc                                                                  >  1:585943/1‑151 (MQ=255)                                                                                 ctaactCAGCACAACTCCACCACCCACTGGTCCAAGGGGTAGCAAGACTCCGGCTTGCCCGCTCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc                                                                  <  1:174824/145‑1 (MQ=255)                                                                                   cAACCACCACAACTTCATCACGCCCTGGGTAAAGGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                               <  2:805384/151‑1 (MQ=255)                                                                                       cAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATa                                                            >  2:870017/1‑150 (MQ=255)                                                                                          cacaACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgt                                                        >  2:997215/1‑151 (MQ=255)                                                                                            caACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                      >  2:184544/1‑151 (MQ=255)                                                                                              aCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                    >  1:1044485/1‑151 (MQ=255)                                                                                               cTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtc                                                   >  2:461928/1‑151 (MQ=255)                                                                                                  catccCGCACGGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                 <  2:1092810/150‑1 (MQ=255)                                                                                                  catcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCg                                                >  1:396665/1‑151 (MQ=255)                                                                                                  catcaCGCACTGGGTCAAAGGGTAGCAAACCTGCGGCGTGACCGCGCTCAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                 <  2:935095/150‑1 (MQ=255)                                                                                                  caccaCCCATTGGGTCAAAGGTTAGCACGACGGCGGCGGGACCGCCAAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                 <  1:558930/147‑1 (MQ=255)                                                                                                   atcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGt                                               >  2:945988/1‑151 (MQ=255)                                                                                                   atcaCGCACGGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtct                                                  <  2:1005640/148‑1 (MQ=255)                                                                                                      aCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                             >  1:277815/1‑150 (MQ=255)                                                                                                          aCTGGGTCAAAGGGAAGCAAGCCGGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaata                                         <  2:1104551/150‑1 (MQ=255)                                                                                                             gggTCAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaata                                         <  1:438282/147‑1 (MQ=255)                                                                                                             gggTAAAAGGGAAGAAAGCCGGGGGGGTGCCCGCCCAAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataat                                       <  1:1173355/149‑1 (MQ=255)                                                                                                                tCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCt                                   >  2:819027/1‑150 (MQ=255)                                                                                                                   aaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACAt                               >  2:1207725/1‑151 (MQ=255)                                                                                                                    aGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                              >  1:814369/1‑151 (MQ=255)                                                                                                                     gggTAGAAAGACGGCGGCGTCACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                              <  1:620087/150‑1 (MQ=255)                                                                                                                            aaaacggcgctgTGACGGCGCAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtttt                         <  1:1036289/137‑1 (MQ=255)                                                                                                                            aaGACTGCGGGGTGACCGCGCACAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGc                      <  2:612401/151‑1 (MQ=255)                                                                                                                                   cggcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaaa                  >  1:1029478/1‑148 (MQ=255)                                                                                                                                    ggcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGt               >  2:253882/1‑150 (MQ=255)                                                                                                                                     gcgTGCCCGCGCAAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTaa             <  1:918956/151‑1 (MQ=255)                                                                                                                                       gTGCCCGCCCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACg           <  1:671400/151‑1 (MQ=255)                                                                                                                                        tgcccGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACg           <  1:927725/147‑1 (MQ=255)                                                                                                                                           ccGCACACAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACggg         <  1:618064/149‑1 (MQ=255)                                                                                                                                              cgcTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAgg     <  2:1186355/150‑1 (MQ=255)                                                                                                                                                cTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTg  >  2:691926/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTG  >  NC_000913/3177603‑3177896

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