Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	48.3%	58.9 / 30.4	61	G566G (GGT→GGG)	gyrA	DNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (24/7);  new base C (0/29);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 7.62e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGT  >  NC_000913/2337573‑2337873                                                                                                                                                       |                                                                                                                                                       acCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTa                                                                                                                                                        <  2:703279/151‑1 (MQ=255)    gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCt                                                                                                                                                     >  1:399684/1‑151 (MQ=255)     cgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCtt                                                                                                                                                    >  1:109637/1‑151 (MQ=255)     cgcgcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCtt                                                                                                                                                    <  2:30629/151‑1 (MQ=255)        gcgcCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTcc                                                                                                                                                 >  1:875776/1‑151 (MQ=255)           ccACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGc                                                                                                                                              >  2:225982/1‑151 (MQ=255)            cACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGcc                                                                                                                                             <  2:512679/151‑1 (MQ=255)            cACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGcc                                                                                                                                             <  2:347346/151‑1 (MQ=255)                        ccGGCAACTGAAAATCTTTCATCAAAAAGCCGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCCCCACCAGCAGTCGGTCGATAAAGTCTTTTTTTTTAATCCGTGGGGAAGATTTCCCTTTCCCGCCACGACGCTgc                                                                                                                                   <  1:491105/149‑1 (MQ=255)                          ggaaaCTGAAAACCTTTCATCAAATAGACGCGACCGCGGCGGGAAAAGCAAAGAATATGGTGGTGATTGTTCCCCCCCAGCGGTGGGTCGAAAAAGTTTTTTTTTTAAATCCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCt                                                                                                                                <  1:824634/147‑1 (MQ=255)                             aaCTGATAACCTTTCATCAAATGGCCGCGCCCCCGGCTGGAGAACCCCAGAATATGTTCGTGTGTGTTCCCCCCCAGCAGTCGTTCAATAAAGTCTTTTTCTTTAATCCGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCtt                                                                                                                               <  2:245701/148‑1 (MQ=255)                             aaCTGATAAACTTTCATCGAATAGCCGCGCCCCGGGCGGGAGAAGCAAAAAATATGGTCGTGAGTGTTCGCCACCAGCAGTGGGTCGATAAAGTCTTCTTTTTTAATGCGTGCGGCAGTTTTCCTTTCCCCGCCACGACGCTGCGCTTCg                                                                                                                             <  1:456962/150‑1 (MQ=255)                               cTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTAt                                                                                                                          >  1:274127/1‑151 (MQ=255)                                    aaccTTTCATCAAATAGACGCGCCCACGGCTGGAGAAGCACAAAATATGGTCGGGAGTTTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATACGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAg                                                                                                                      <  2:170939/147‑1 (MQ=255)                                    aaaCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTc                                                                                                                        >  1:300076/1‑148 (MQ=255)                                      aCTTTCATCGAATAGACGCGCCCACGGCGGGAGAACCACAAAATAGGGTGGTGAGTGTTCCCCCCCAGCAGTGGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaaa                                                                                                                   <  2:423525/151‑1 (MQ=255)                                        tttCATCAAATAGAGGGGCCCACGGTTGGAGAAGAACAAAATATGGTGGGGAGTGTCCGCCCCCAGCAGTCGTTCGAAAATGTCTTCTTTTTTAATACGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaaa                                                                                                                   <  2:343332/149‑1 (MQ=255)                                        tttCACCGAAAGGACGCGCCCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTCCCCCCCCAGCAGTCGTTCGAAAAAGTCTTCTTTTTTAATACGTGCGGAAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaaa                                                                                                                   <  2:980025/149‑1 (MQ=255)                                          tCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCgg                                                                                                               >  1:755093/1‑151 (MQ=255)                                             tCGAAAAGCCGCGCCCCCGGCGGGAGAAGCAAAGATTTTGGTGGGGAGTGTCCCCCCCAAGCAGTGGTTGAAAAAAGTCTCCTTTTTAAATACGTGCGGAATATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                             <  1:116933/150‑1 (MQ=255)                                                    gccgcgACCCCGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTTTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa                                                                                                     <  2:56472/149‑1 (MQ=255)                                                        cgcccACGGCTGGAAAAGCACAGAATAGGGTCGTTAGTGTTCGCCACCAGCAGTGGGTCGATAAAGTCTTCTTCTTTAATCCGGGGGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACtt                                                                                                      <  2:816663/143‑1 (MQ=255)                                                        cgACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                                                 >  1:389866/1‑151 (MQ=255)                                                           ccACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTCCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAg                                                                                               <  2:972045/150‑1 (MQ=255)                                                                 aaagaaacgCACAGAATTGGCTGGGGAGTGTCCGCCAACACAAGGGGGTCGATAACGTTTTTTTATTTAAAGTATGCGAAAGATTGCCGTTCGCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                        <  1:30964/143‑1 (MQ=255)                                                                   ggAGAAGCACAGAATATGGTGGGGGGTGTTCCCCCCCAGCATGCGGTCGAAAAATTCTTCTTTTTTAATCCGTGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                        <  2:389583/149‑1 (MQ=255)                                                                   ggAGAAGCACAGAATATGGTCGTGAGTGTCCGCCACCAGCAGTCGGTCGATAAAGTTTTCTTTTTTAATACGTGCGGAAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGt                                                                                        <  1:575799/149‑1 (MQ=255)                                                                        aGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagag                                                                                   >  1:1071169/1‑149 (MQ=255)                                                                           acaGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGt                                                                              >  1:266631/1‑151 (MQ=255)                                                                            caGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTc                                                                             >  1:616514/1‑151 (MQ=255)                                                                             aGAATATGGTCGTGAGTGTTCCCCACCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATGCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGt                                                                              <  1:560300/149‑1 (MQ=255)                                                                               aaTATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACg                                                                          <  2:193657/151‑1 (MQ=255)                                                                               aaTATGGTCGTGAGTGTCCGCCCCCAGCAGTCGGTCAATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTc                                                                             <  2:789573/148‑1 (MQ=255)                                                                                        gTGAGTGTTCGCCCCCAGCTGGCGTTCTATAACGCTTTTTTTTTAAATACGTGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACAt                                                                   <  1:279553/149‑1 (MQ=255)                                                                                           aGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATc                                                                  >  1:207829/1‑147 (MQ=255)                                                                                             tgtTCGCCCCCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATCCGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                                             <  2:51659/150‑1 (MQ=255)                                                                                                  gCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTg                                                        >  1:803089/1‑150 (MQ=255)                                                                                                   ccaccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGAt                                                      >  1:139203/1‑151 (MQ=255)                                                                                                    caccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                     >  2:387708/1‑151 (MQ=255)                                                                                                       cagcagTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTTACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                     >  2:211972/1‑148 (MQ=255)                                                                                                          cagTCGGTCAAAAAAGTCTTTTTTTTTAATACGGGGGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGAt                                                 <  1:299364/149‑1 (MQ=255)                                                                                                                gtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                         >  2:1103099/1‑151 (MQ=255)                                                                                                                 tcgATAAATTTTTCTTTTTTAATCCGTGCGGCAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                         <  1:837454/150‑1 (MQ=255)                                                                                                                 tcaaTAAATTTTTTTTCTTTAATACGTGCGGCAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTcc                                            <  1:970410/144‑1 (MQ=255)                                                                                                                  cgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt                                       >  2:601493/1‑151 (MQ=255)                                                                                                                    aTAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGa                                     <  1:861281/151‑1 (MQ=255)                                                                                                                      aaaGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg                                   >  1:667506/1‑151 (MQ=255)                                                                                                                         tttttcttcttAAATAGGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTc                                 <  1:248265/147‑1 (MQ=255)                                                                                                                            ttcttcttAAATCCGGGCGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTc                                 <  1:33140/147‑1 (MQ=255)                                                                                                                               ttcttAAATACGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg                             <  2:762738/148‑1 (MQ=255)                                                                                                                                tttttAATACGGGCGGAATTTTTTCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTg                          <  1:415274/148‑1 (MQ=255)                                                                                                                                tcttTAATTCGTGCGGAAGTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTg                          <  1:802553/150‑1 (MQ=255)                                                                                                                                  tttAATACGGGCGGCAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGtt                        <  1:340164/150‑1 (MQ=255)                                                                                                                                   ttAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTgg                      >  1:586818/1‑151 (MQ=255)                                                                                                                                   ttAATACGTGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt                         <  1:823417/148‑1 (MQ=255)                                                                                                                                   ttAATACGGGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGt                         <  2:635984/148‑1 (MQ=255)                                                                                                                                     aaTACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCg                    >  1:451600/1‑151 (MQ=255)                                                                                                                                       tACGTGCGGCATTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCgg                   <  2:511293/150‑1 (MQ=255)                                                                                                                                          gtggggaataTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGa                <  1:968851/141‑1 (MQ=255)                                                                                                                                             cGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTc            >  2:762772/1‑151 (MQ=255)                                                                                                                                                 aGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGt         >  2:154043/1‑150 (MQ=255)                                                                                                                                                  gATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTac       >  1:909303/1‑151 (MQ=255)                                                                                                                                                       aCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGt  >  1:362680/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       ACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGT  >  NC_000913/2337573‑2337873

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