Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	40.0%	138.8 / 20.0	76	G566G (GGT→GGG)	gyrA	DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (41/4);  new base C (0/30);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 5.93e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACG  >  NC_000913/2337575‑2337872                                                                                                                                                     |                                                                                                                                                       cgcgcgcgCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAcc                                                                                                                                                      >  1:81211/1‑151 (MQ=255)  gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCt                                                                                                                                                     >  1:1107708/1‑151 (MQ=255)  gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCt                                                                                                                                                     <  1:634503/151‑1 (MQ=255)        gcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCg                                                                                                                                              <  1:156195/151‑1 (MQ=255)              aGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCacg                                                                                                                                         >  2:610974/1‑150 (MQ=255)                  gCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCt                                                                                                                                    >  2:826335/1‑151 (MQ=255)                         gCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTc                                                                                                                             >  2:586468/1‑151 (MQ=255)                         gCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTc                                                                                                                             >  1:457158/1‑151 (MQ=255)                           aaCTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGt                                                                                                                           >  1:305726/1‑151 (MQ=255)                           aaCTGATAAACTTTCATCGAATAGACGCGACCACGGCGGGAGAAGCAAAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGGGCGGCAGATTTCCCTT‑TCCCGCCACGACGCTGCGCtt                                                                                                                              <  2:501274/148‑1 (MQ=255)                            aCTGAAAAACTTTCAGCGAATAGACGCAACCCCGGCTGGAGAAGCACAGAATATGGTGGGGAGTGTTCCCCCCCAGCAGTCGGTCAAAAAAGTTTTCTTCTTTAATCCGGGCGAAAGATTTCCCTT‑CCCCGCCACGACGCTGCGCTTc                                                                                                                             <  1:780416/148‑1 (MQ=255)                                ataatttttCATCGAATAGACGCGACCCCGGCTGGAGAAGCAAAAAATATGGTCGTGAGTTTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                       <  2:362189/144‑1 (MQ=255)                                aTAATCTTTCATGGAATGGCCGCACCCACGGTGGGAGAACCACAAAATATGGTCGTGCGTGTCCCCCCCCAGCAGTCGGTCGATAAATTTTTCTTCTTTAACACGTGCGGGAAATTTTCTTT‑TCCCGCCACGACGCTGCGCTTc                                                                                                                             <  2:543679/144‑1 (MQ=255)                                  aaccTTCCATCGAATAGCCGCGACCACGGTTGGAGAAGCAAAGAATATGTTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTTTTTTTTAATACGTGGGGCAGATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAg                                                                                                                     <  2:517814/147‑1 (MQ=255)                                  aaaCTTTCATCAAATAGCCGCCACCACGCGGGGAGAACCACAAAATTGGGTGGGGGGTTTCCCCCCCCAGCAGGCGGTCAATAAGGTCTTCTTCTTAATTACGTGGGGAAGTTTTCCTTT‑TCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                       <  1:122528/148‑1 (MQ=255)                                   aaCTTTCATCGAATAGACGCGCCCACGGCTGGAAAAGCACAAAATATGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGaa                                                                                                                   <  2:143663/151‑1 (MQ=255)                                   aaCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGa                                                                                                                    >  1:612846/1‑150 (MQ=255)                                     cTTTCATCGAATAGCCGCGCCCGCGGCTGGAGAAGCACAAAATATGGTCGTGATTGTTCGCCACCAGCAGTCGGCCTATAAAGTCTTCTTCTTTAATGGGTGGGGCAGTTTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGaa                                                                                                                   <  1:241979/149‑1 (MQ=255)                                          aTCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                            >  2:896417/1‑151 (MQ=255)                                           tCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGGGAGTGCTCGCCACCAGCAGTCGGTCGATAAAGTATTCTTCTTTAATACGGGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAg                                                                                                                     >  1:97730/1‑141 (MQ=255)                                              aaaagccgcgCCCACGGCTGGAGAAGCAAAGAATATGGTGGTGAGTTTTCCCCCCCAGCAGTGGGTCGAAAAAGTCTTTTTTTTTAACACGTGCGGAAGATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTg                                                                                                           <  1:793828/142‑1 (MQ=255)                                                      cgACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                                >  1:293086/1‑151 (MQ=255)                                                      cgACCACGGCTGGAGAAGCACAAAATATGGTCGTGAGTTTTCGCCACCAGCAGTGGGTCGATAAAGTTTTTTTTTTTAATACGTGCGGCAGATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACg                                                                                                 <  2:365477/150‑1 (MQ=255)                                                       acccACGGCGGGAGAAGCACAGAATAGGGCCGTGAGTTTCCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACg                                                                                                 <  1:266282/147‑1 (MQ=255)                                                             ggCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTg                                                                                         >  2:691627/1‑151 (MQ=255)                                                                tGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                      >  2:907429/1‑151 (MQ=255)                                                                 ggAGAAGCACAAAATATGGTTGTAAGTGTTCCCCCCGAGTTGTCGTCCAATAATTTTTTCTTTTTAAATCCG‑GCGGCAAGATTACCTTCTTCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTga                                                                                     <  2:722965/151‑1 (MQ=255)                                                                 ggAAAAGCACATAATTTGGTGGGGGTTGCCCCCCCCCAGCTGTCTGTCGATAATTTTTTTTTCTTTATCGCGGCGGGCACTTTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                      <  1:816201/150‑1 (MQ=255)                                                                   agaAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                      >  1:1028253/1‑148 (MQ=255)                                                                         acaGAATATGGTCGTGAGTTTTCCCCCCCAGCAGTCGGTGGATAAAGTCTTTTTCTTTAATCCGTGCGGAAGATTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCg                                                                              <  2:502354/150‑1 (MQ=255)                                                                           aaattatGGTCGTGAGTGTCCGCCACCAGCAGTCGTTCGATAAAGTTTTCTTTTTTAATGCGTGGGGCATATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTc                                                                            <  1:22074/146‑1 (MQ=255)                                                                           aaaaTATGGTCGGGAGTGTCCGCCACCAGCAGTCGGTCGAAAAAGTCTTCTTCTTTAATACGTGCGGCAGTTTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                                           <  2:1071501/149‑1 (MQ=255)                                                                             aaattGGTCGGGGGTGTTCGCCCCCCGAATTCGGTCAATAAAGTTTTTTTTTTTAATCCGTGCGGAAGTTTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACg                                                                         <  1:834768/147‑1 (MQ=255)                                                                                atGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGAcc                                                                      >  2:703586/1‑151 (MQ=255)                                                                                atGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGAcc                                                                      >  1:502293/1‑151 (MQ=255)                                                                                  ggtggtGAGTTTTCCCCACCAGCAGTCGGTCGAAAAAGTCTTCTTCTTTAATACGTGCGGCAGTTTTCCTTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACca                                                                     <  1:849822/150‑1 (MQ=255)                                                                                        gAGTGTTCCCCCCCAGCAGTCGGTCGATAAATTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCtt                                                               <  2:310659/150‑1 (MQ=255)                                                                                        gAGTGTTCACCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTc                                                              >  2:181118/1‑151 (MQ=255)                                                                                             ttCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                                                         >  2:775282/1‑151 (MQ=255)                                                                                               cGCCACCAGCAGTGGGTCGATAAATTCTTCTTTTTTAATACGTGCGGAAGTTTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTTACGACCACATCTTCCTGGGTg                                                       <  1:226812/151‑1 (MQ=255)                                                                                                 cccccaGCAGTCGGTGGATAAATTCTTCTTCTTTAATCCGTGCGGCAGATTCCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGt                                                        <  2:744232/145‑1 (MQ=255)                                                                                                  caccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                    >  1:808150/1‑151 (MQ=255)                                                                                                        cagTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCt                                              >  2:16447/1‑151 (MQ=255)                                                                                                        cagTCGGTCAATAAAGTTTTTTTTTTTAATACGTGCGGCAGTTTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAg                                                  <  1:409555/147‑1 (MQ=255)                                                                                                         agTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCtt                                             >  2:260561/1‑151 (MQ=255)                                                                                                             ggtcaAAAAAGTCTTTTTTTTTAATGCGTGCGGCAAATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAg                                         <  2:964067/151‑1 (MQ=255)                                                                                                              gtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                        >  2:440209/1‑151 (MQ=255)                                                                                                              gtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAg                                         >  2:170048/1‑150 (MQ=255)                                                                                                                cgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt                                      >  2:1093244/1‑151 (MQ=255)                                                                                                                    aaaGTCTTCTTTTTTAATCCGTGCGGAAGTTTTCCTTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                                  <  2:175718/151‑1 (MQ=255)                                                                                                                    aaaGTCTTCTTTTTTAATACGTGCGGCAGATTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAt                                   <  2:612947/150‑1 (MQ=255)                                                                                                                    aaaGTCTTCTTTTTAAATACGTGCGGCAGATTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt                                      <  2:388412/147‑1 (MQ=255)                                                                                                                    aaaGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                                  >  2:736626/1‑151 (MQ=255)                                                                                                                     aaGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGt                                 >  2:286450/1‑151 (MQ=255)                                                                                                                      aGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTc                                >  1:311382/1‑151 (MQ=255)                                                                                                                       gTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCt                               >  1:328837/1‑151 (MQ=255)                                                                                                                       gTCTTCTTCTTAAATACGGGCGGAATTTTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt                                      <  1:836493/144‑1 (MQ=255)                                                                                                                            cttcttTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCt                          >  2:1098672/1‑151 (MQ=255)                                                                                                                              tcttTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt                        >  1:468340/1‑151 (MQ=255)                                                                                                                                tttAATCCGTGCGGAATATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt                        <  1:877151/149‑1 (MQ=255)                                                                                                                                tttAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                      >  2:47407/1‑151 (MQ=255)                                                                                                                                  tAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGc                    >  2:970797/1‑151 (MQ=255)                                                                                                                                  tAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGc                    >  1:1105272/1‑151 (MQ=255)                                                                                                                                   aatgcggggggCAGTTTTCCTTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCg                   <  1:92572/142‑1 (MQ=255)                                                                                                                                   aaTGCGGGCGGCAATTTTCCCTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCg                   <  1:1019634/151‑1 (MQ=255)                                                                                                                                            ggCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTc           >  2:167760/1‑150 (MQ=255)                                                                                                                                            ggCAAATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTc           <  2:266214/150‑1 (MQ=255)                                                                                                                                               agttttccctt‑tcccGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTc           <  1:420898/144‑1 (MQ=255)                                                                                                                                               aGATTTCCTTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAg         <  2:533626/149‑1 (MQ=255)                                                                                                                                               aGATTTCCCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGt        <  2:192058/150‑1 (MQ=255)                                                                                                                                                gATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGCGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGGTTTCAGTa       >  1:65303/1‑150 (MQ=255)                                                                                                                                                gATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTa       >  2:800100/1‑150 (MQ=255)                                                                                                                                                 aTTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacg     >  1:717149/1‑151 (MQ=255)                                                                                                                                                 aTTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacg     >  1:960029/1‑151 (MQ=255)                                                                                                                                                  tttccctt‑tcccGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTac      <  1:714724/149‑1 (MQ=255)                                                                                                                                                  tttACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacg     >  1:954376/1‑150 (MQ=255)                                                                                                                                                    tACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacgacg  >  2:148403/1‑151 (MQ=255)                                                                                                                                                    tACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGGGATTTCAGTacgacg  >  1:517169/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTT‑TCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACG  >  NC_000913/2337575‑2337872

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