Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	44.6%	118.6 / 14.0	65	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/2);  new base C (0/29);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.85e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGA  >  NC_000913/2984292‑2984588                                                                                                                                                      |                                                                                                                                                    gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                    >  1:925057/1‑151 (MQ=255) gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                    <  2:269304/151‑1 (MQ=255)  aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAATTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                   >  2:1563479/1‑151 (MQ=255)  aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                   >  2:1477009/1‑151 (MQ=255)   atatGTCATGAGAATGGAAAAATGAGGGAGGCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTTAATGCCacacac                                                                                                                                                  <  2:508792/151‑3 (MQ=255)   atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                  <  2:600755/151‑1 (MQ=255)              gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCTccgcc                                                                                                                                       <  1:1125031/151‑1 (MQ=255)                     aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACa                                                                                                                                 >  2:517145/1‑150 (MQ=255)                     aaaaGGGGGGAGCCAAAAAAATTTATGTCGATTGGCCAGCAATCATTATTCAACAGTTAAATGACGGCAGATATTTATCGGTAGGGAGAAAATGCAGAGCTTTTCTTCACGTTAAATGGAAAATGCCACCCCCGGGCCCCCGCCGATACaa                                                                                                                                <  2:1413518/151‑1 (MQ=255)                        ctGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGc                                                                                                                             >  2:1247781/2‑151 (MQ=255)                          gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                                                                           >  1:706871/1‑151 (MQ=255)                                  ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                   >  1:1563672/1‑151 (MQ=255)                                  ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                   >  1:81924/1‑151 (MQ=255)                                      aaattttttgtCGACTGCCCACCCACAATAATACACCAAATAAATGCCGGTAGTTATTTAGGGGAAGGGGGAAAAGGCAGACCTTTTCGTCACGTTCATGGAACAGTGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                               <  2:298789/143‑1 (MQ=255)                                      aaaaTGTATGCCGCTTGGCCACCTATATTTTTTCAACAGATTATTGCGGGAAGATATTTTGCGGTAGGGAGGAAAGGCAGAGCTATTCTTCCGTTTTACTGGTCAGGCCCCCCCCCTATCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                                    <  1:1069698/146‑1 (MQ=255)                                       aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                              >  1:919697/1‑151 (MQ=255)                                       aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                              >  2:160382/1‑151 (MQ=255)                                          tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                           >  2:551069/1‑151 (MQ=255)                                              tgtCAACTTGCCGGCCATTTTTTTTCAACAGTTAAAGGGCGCAAAATATTTTTGGGAAAGAGGGAAAGGCAAAGTTTTTCCCCACTTTAAAGGGTCAGTCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                                        <  1:1446345/150‑1 (MQ=255)                                                 cgccTTGCCAGCCATCATTATTAAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTCAATGGTAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                     <  1:316592/147‑1 (MQ=255)                                                 cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                    >  2:574260/1‑151 (MQ=255)                                                  gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                   >  1:1389146/1‑151 (MQ=255)                                                  gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                   >  2:850821/1‑151 (MQ=255)                                                   aCTTGCCGCCCATCATTATTAAAAAAATTAATGCCGGCAATTATTTAGCGGTAAGGAGAAAGTGAGGGTTTTTTCTTCAGGTTCAATGGTCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                      <  2:745186/147‑1 (MQ=255)                                                      tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTc                                                                                               >  2:1515265/1‑151 (MQ=255)                                                        ccaccctTCATTATTCAAAAAATTAATGCCGGCAGATTTTTATCGGTAGGGAGAAAATGCAGAGTTTTTGGCCCCTTTAAGTGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                 <  2:1343988/147‑1 (MQ=255)                                                        caaccaatcatTATTAAACAGATTAATGACGGCAGATTTTTATCGGTAAGGAGAAAATGAAGAGCTTTTCGTCCCTTTAAATGGACAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                <  1:95767/142‑1 (MQ=255)                                                            ccaTCATTATTCAACAGATTAATGAGGGCAGATATTTAGGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                            <  1:78339/148‑1 (MQ=255)                                                                catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                     >  1:1389927/1‑151 (MQ=255)                                                                  ttattCAACAGTTAAAGGCCGGCAGATATTTATCGGTAAGGGGGAAATGCAGAGCTTTTCGTCACTTTCAAGGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                      <  2:523691/148‑1 (MQ=255)                                                                  ttataaaaCAGATAAATGACGGCAAATTTTTACGGGAAAGGAGAAAATGCAGAGCTTTTCTTCACGTTAAATGGTAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                       <  1:600491/141‑1 (MQ=255)                                                                    attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                 >  1:106631/1‑151 (MQ=255)                                                                    attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                 >  1:653265/1‑151 (MQ=255)                                                                     ttCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                >  2:1059381/1‑151 (MQ=255)                                                                        aagaGATTAATCACGGCGGATTTTTTGCGTTAAGGAGGAAAGGCAGAGTTTTTATTCCCGTTCAAGGGATAATGCCTCCCCCTTTCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                        <  1:563531/137‑1 (MQ=255)                                                                        aaCAGATTAATGCCGGCAAATTTTTATGGGAAAGGAGAAAAGGAAGAGTTTTTCTTCACTTTAAATGTTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                 <  1:86829/147‑1 (MQ=255)                                                                         aCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGCCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                              >  2:979336/1‑149 (MQ=255)                                                                            gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                          >  1:998242/1‑150 (MQ=255)                                                                                 aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                    >  2:1143215/1‑151 (MQ=255)                                                                                 aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                    >  1:849102/1‑151 (MQ=255)                                                                                  tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                    >  1:5176/1‑150 (MQ=255)                                                                                      ggaaGATTTTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTAAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                  <  1:394863/145‑1 (MQ=255)                                                                                      ggCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCATTTTAAATGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                   <  2:446111/147‑1 (MQ=255)                                                                                        aaGATATTTATGGGTAAGGAGAAAAGGCAGAGCTTTTCGTCACGTTCAATGGTAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                               <  1:271007/148‑1 (MQ=255)                                                                                          gttttttAGCGGTAAGGAGAAAATGCAGAGTTATTCGTCACGTTAAATGGTCAACGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                             <  2:1300384/145‑1 (MQ=255)                                                                                            tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                         >  2:1189784/1‑151 (MQ=255)                                                                                            tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                         >  1:496963/1‑151 (MQ=255)                                                                                                tATCGTTAAGGAGGAAATGCAGAGTTTTTCGTCACGTTAAATGGTCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                      <  2:904765/150‑1 (MQ=255)                                                                                                tATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                      >  1:701588/1‑150 (MQ=255)                                                                                                    ggTAAGGAGGAAATGCAGAGTTTTTCGTCACGTTCAATGGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                   <  1:341165/149‑1 (MQ=255)                                                                                                     gTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                                >  1:612177/1‑151 (MQ=255)                                                                                                       agggagaaAATGAAGACTTTTTCGCCACTTTAAATGTTAAGTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                               <  2:1282520/148‑1 (MQ=255)                                                                                                         ggaggaAATGCAGGGCTTTTCGTCCTTTTCAATGTTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                            <  2:1399970/151‑1 (MQ=255)                                                                                                         ggagaaAAGGAAGAGTTATTCGTCACGTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                            <  2:759446/151‑1 (MQ=255)                                                                                                           agaaaaTGAAGAGCTTTTCGTCACTTTCAATGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                           <  1:296535/147‑1 (MQ=255)                                                                                                                aTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                     >  2:431331/1‑151 (MQ=255)                                                                                                                 tGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                    >  1:988514/1‑151 (MQ=255)                                                                                                                  gCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                   >  1:1053651/1‑151 (MQ=255)                                                                                                                           ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                          >  1:853210/1‑151 (MQ=255)                                                                                                                            tCGTCACGTTCAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                          <  1:1438566/150‑1 (MQ=255)                                                                                                                               tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                      >  1:663180/1‑151 (MQ=255)                                                                                                                               tCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                      >  1:737197/1‑151 (MQ=255)                                                                                                                                cACGTTAAATGGTAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                     <  1:746784/151‑1 (MQ=255)                                                                                                                                    tacAATGATCAATCCCCCCCCCTCCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccac                   <  2:50099/147‑1 (MQ=255)                                                                                                                                       aaTGGTCAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt               <  2:549361/150‑1 (MQ=255)                                                                                                                                           gTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat          >  1:1220234/1‑151 (MQ=255)                                                                                                                                              aaTGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattga       <  2:191053/151‑1 (MQ=255)                                                                                                                                              aaTCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat          <  2:398246/148‑1 (MQ=255)                                                                                                                                                   cacaCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGAGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg   >  2:4490/1‑150 (MQ=255)                                                                                                                                                    acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa  >  2:1208411/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGA  >  NC_000913/2984292‑2984588

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