Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,835,219	0	A	C	42.2%	107.1 / 22.5	64	V736G (GTT→GGT)	hypF	carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base A (29/8);  new base C (0/27);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 4.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGAT  >  NC_000913/2835083‑2835361                                                                                                                                         |                                                                                                                                               gcgcGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGcccc                                                                                                                                     >  1:91930/1‑148 (MQ=255)    cgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGAc                                                                                                                               >  2:1101645/1‑151 (MQ=255)     gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACa                                                                                                                              <  1:170951/151‑1 (MQ=255)          aTCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAAcc                                                                                                                          <  1:248903/150‑1 (MQ=255)            ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGc                                                                                                                        <  2:265818/150‑1 (MQ=255)            ccGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGc                                                                                                                        <  1:44329/150‑1 (MQ=255)             cGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGcc                                                                                                                       >  1:177187/1‑150 (MQ=255)              gACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTCCACCCGCTAACCAACGCGCCGCAGCAATACCCCCCTGCCCCAGAGACAAACCGCCa                                                                                                                      <  2:694157/150‑1 (MQ=255)                ccTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATc                                                                                                                    >  2:1146385/1‑150 (MQ=255)                     ggggACGCATGTGTAGACCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCGCCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCAccc                                                                                                                <  1:812181/149‑1 (MQ=255)                       ggACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgccg                                                                                                            >  1:702865/1‑151 (MQ=255)                        ggcGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCgg                                                                                                           <  1:761198/149‑1 (MQ=255)                        gACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCgg                                                                                                           <  1:240642/151‑1 (MQ=255)                                    ggCGGGATAAGGCTTTTGCGCGGCTCCGGGCAATGGTGCCAAAATCCACCCCGTTTTATCCTTTCTGGACTTCCCCCGCTACCCACCGCGCCGCAGAAATACCCCCCGGCCCAAGAGAAAAACCGCCACCCCCCGCCGGTAAACTCTGTgg                                                                                               <  2:395487/151‑1 (MQ=255)                                     gCCGGATAAGGCGTTTCCCCCGCATCCGGCAAGGGTGCCAAAATGCAACACTTTTTACCCTTTCTGCTTTTCCCCCGTTAACCACCGCGCCGCACAAAACCCCCCCGGCCCAAGAGACAAACCCCCATCCCCCGCCGGTAAACTCTGTGGa                                                                                              <  1:1192485/151‑1 (MQ=255)                                     gCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTgg                                                                                               >  1:281497/1‑150 (MQ=255)                                        ggATAGGGCTTTTACGCCGCATCGGGCAATGTTGTCCAAATGCAACCCTTTTTTTCCTTTCTGGCCTCCCCCCGCTACCCAACGCGCCGCAGCAATACCCCCCTGCCCAAGAGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGaa                                                                                             <  1:657771/149‑1 (MQ=255)                                         gATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGaa                                                                                             >  2:309904/1‑148 (MQ=255)                                          aTAGGGCTTTTACGCCGCTTCCGGCATTGGTGTCAAAATGCAACCCGTTTTATCCGTTCTGGACTTCCCCCGCTAACCACCGCGCCGCAGCAATACCCCCCTGCCCCAGAGACAACCCCCCATCCCCCGCCGGTAAACTCTGTGGAAag                                                                                           <  1:783244/149‑1 (MQ=255)                                            aggggtttccccccgcACCCGGAATTGGTGCCCAAATCAACCACGTTTTACCCTTTCTGGCCTCCCCCCCCAAACCACCGCGCCGCACAAATACCCCCCTCCCCAAGAGCAAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGc                                                                                        <  2:294639/138‑1 (MQ=255)                                              ggCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAAt                                                                                     >  2:289484/1‑151 (MQ=255)                                                cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGCGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAAtg                                                                                    >  2:194943/1‑150 (MQ=255)                                                cGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAAtgt                                                                                   >  1:1102351/1‑151 (MQ=255)                                                 gTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAAtgtg                                                                                  >  1:1060948/1‑151 (MQ=255)                                                   ttACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGa                                                                                 >  2:763616/1‑150 (MQ=255)                                                    tACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGaa                                                                                >  2:673630/1‑150 (MQ=255)                                                      cgccccACCCGGAATGGGTTTCAAAATCAAACCCTTTTTACCCGTTCGGGCCTCCCCCCGCAACCCACCGCCCCGAAGAAACACCCCCCCCCCCAAAAAAAAACCCCCCACCCCCCGCCGATAAACTCTGTGGAAAGAGCAATGTGAAATc                                                                             <  1:868268/151‑1 (MQ=255)                                                       gccgcATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATc                                                                             >  1:473908/1‑150 (MQ=255)                                                       gccgcATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTCCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATc                                                                             >  2:605462/1‑150 (MQ=255)                                                        cggcaccctccAATGGTGCCAAAATGCACCCTTTTTTACCGGTTCTGCACTTCCCCCGCAACCCAACGCCCCGCAGAAATACCCCCCTGCCCAAGGTAAAAACCCCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCa                                                                            <  1:651178/140‑1 (MQ=255)                                                        ccgcACCCGGAAATGGTGTCAAAAGGAAACCCTTTTTATCCGTTCTGGCCTTCCCCCGAAAACAACCGCGCCGCAGAAAAACCCCCCTGCCCCAGAGAAAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAg                                                                           <  1:736410/151‑1 (MQ=255)                                                         cgcATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGc                                                                          >  2:1157285/1‑151 (MQ=255)                                                         cgcATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGc                                                                          >  1:312765/1‑151 (MQ=255)                                                          gcATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCg                                                                         >  2:549151/1‑151 (MQ=255)                                                          ccccccGACAATGTTGCCAAAATGCAACCCTTTTTATCCTTTCGGGACTCCCCCCGCAACCCACCGCCCCGCAGAAATACCCCCCGCCCCCAGAAAAAACCCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCg                                                                         <  1:1091743/147‑1 (MQ=255)                                                           cATCCGGCAATGGTGTCAAAAGGCAACACGTTTTATCCGTTCTGGACTCCCCCCGCTAACCACCGCGCCGCAGCAATACCCCCCTCCCCCAGAGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCga                                                                        <  1:928363/151‑1 (MQ=255)                                                              ccgaaaaTGGTGTCAAATTGAAACACTTTTTACCCGTTCTGGCCTCCCCCCGCTAACCACCGCGCCGCAGCAATACCCCCCTGCCCCAGAGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCgag                                                                       <  2:830912/144‑1 (MQ=255)                                                                ggCATTGGTGTCAAAATGCAACACTTTTTATCCGTTCTGGCCTCCCCCCGCAAACCACCGCGCCGCAGAAATACCCCCCGGCCCAAGAGAAAACCCCCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATa                                                                    <  2:1136169/150‑1 (MQ=255)                                                                  cAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATg                                                                 >  2:1074859/1‑151 (MQ=255)                                                                  cAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATg                                                                 >  1:645757/1‑151 (MQ=255)                                                                   aagggtgcaaaaaGGCAACCTGTTTTATCCGTTCGGGCCTTCCCCCGTAAACCAACGCGCCCCAGAAAAACCCCCCCGCCCCAGAGAAAAACCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATgc                                                                <  1:500708/141‑1 (MQ=255)                                                                   aaTGGTGCCCAAATCCAACACTTTTTACCCGTTCGGTTCTCCCCCCGCTAACCACCGCCCCGCAGCAATACCCCCCCCCCCAAGAGACAACCCGCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAAt                                                                  <  1:496330/149‑1 (MQ=255)                                                                   aaTGGTGCCAAAAGGAAACCCTTTTTACCCGTTCTGGCCTCCCCCCGCTAACCAACGCGCCGCAGCAATAACCCCCCCCCCAAGAAAAAACCCGCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAAt                                                                  <  2:352899/149‑1 (MQ=255)                                                                      ggTGGCCAAATCAAAAACTTTTATTCCTTTTTGACCTCCCCCCGCTACCCACCGCGCCGAAAAAATACCCCCCTCCCCCATAGACAACCCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATgcg                                                               <  1:871107/149‑1 (MQ=255)                                                                        tgtCAAAATGCAACCCTTTTTATCCGTTCGGGCCTTCCCCCGCAACCCACCGCGCCGCAGAAATACCCCCCTGCCCAAGAGAAAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATgcg                                                               <  1:713340/147‑1 (MQ=255)                                                                           caaaaTGCAACCCGTTTTATCCTTTCGGGACTCCCCCCGCTACCAACCGCCCCGCAGAAATACCCCCCTGCCCCAGAGACAACCCGCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACg                                                        <  1:36063/149‑1 (MQ=255)                                                                             aaaTGCACCACTTTTTTCCCTTTCTGGATTCCCCCCGCTACCCACCGCGCCGCAGAAATACCCCCCGCCCCAAGAGAAAACCCGCCACCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGt                                                       <  2:993646/150‑1 (MQ=255)                                                                             aaaTGCAACCCGTTTTCTCCTTTCTGGCCTCCCCCCGCTACCCAGCGCGCCGCAGAAAAACCCCCCTCCCCAAGAGACAACCCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACg                                                        <  1:695264/149‑1 (MQ=255)                                                                             aaaTGCAACCCGTTTTATCCGTTCTGGACTTCCCCCGCTAACCAACGCGCCGCAGCAATACCCCCCTCCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTg                                                      <  1:1136743/151‑1 (MQ=255)                                                                              atggaaaCACTTTTTACCCTTTCGGGACTCCCCCCGCAACCCAACGCGCCGCAGAAAAACCCCCCTCCCCCAGAGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGc                                                     <  1:238324/146‑1 (MQ=255)                                                                                     ccctttttATCCTTTCTGCTCTCCACCCGCTACCCACCGCGCGGCAGCAATACCCCCCTGCCCCAGAAAAAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACgcag                                               <  1:228897/146‑1 (MQ=255)                                                                                      acGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACgcagc                                              >  1:1176463/1‑150 (MQ=255)                                                                                         ttttATCCGTTCTGTACTCCCCCCGTAAACCAACGGCCCGAAACAATACCCCCCCACCCAAGAGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACgcag                                               <  2:861606/146‑1 (MQ=255)                                                                                         ttttACCCGTTCGGGCCTCCCCCCCCAAACCACCGCCCCGCAGAAACCCCCCCCGCCCCCAGAGAAAACCCCCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCaaa                                           <  2:109580/150‑1 (MQ=255)                                                                                               ccGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTAt                                    >  1:37488/1‑151 (MQ=255)                                                                                               ccGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTAt                                    >  2:665186/1‑151 (MQ=255)                                                                                                    cTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                >  2:1111051/1‑150 (MQ=255)                                                                                                      ggcctccccccGCAACCCACCGCGCCGCAGCAATACCCCCCTCCCCAAGAGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATa                               <  1:22550/141‑1 (MQ=255)                                                                                                      ggACTCCCCCCGCTAACCACCGCGCCGCAGAAATACCCCCCTGCCCCAGAGACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                <  1:266034/148‑1 (MQ=255)                                                                                                          ttCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAAcccc                          >  2:6763/1‑150 (MQ=255)                                                                                                            cACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCccgcc                       >  1:517279/1‑151 (MQ=255)                                                                                                                        aaCGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATAcc             >  2:24018/1‑149 (MQ=255)                                                                                                                                  agcaATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt  >  1:604311/1‑150 (MQ=255)                                                                                                                                  agcaATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt  >  2:206798/1‑150 (MQ=255)                                                                                                                                         |                                                                                                                                               GCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCTAACCAACGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGAT  >  NC_000913/2835083‑2835361

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