Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,531	0	A	C	42.4%	105.3 / 10.3	59	V127G (GTG→GGG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (25/9);  new base C (0/25);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 2.19e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATC  >  NC_000913/3177384‑3177680                                                                                                                                                     |                                                                                                                                                      gCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCg                                                                                                                                                    >  2:882530/1‑151 (MQ=255)  ccTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCgg                                                                                                                                                   <  2:24211/151‑1 (MQ=255)       tCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTgg                                                                                                                                              >  1:399443/1‑151 (MQ=255)            gCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGt                                                                                                                                         <  1:1088772/151‑1 (MQ=255)               tACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTg                                                                                                                                      >  2:206835/1‑151 (MQ=255)                aCCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTg                                                                                                                                      <  1:888671/150‑1 (MQ=255)                    caTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTAt                                                                                                                                 <  2:254208/151‑1 (MQ=255)                        aaCGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                                                                                                             >  1:1103581/1‑151 (MQ=255)                             gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                                                                                                                        <  2:911057/151‑1 (MQ=255)                                 atCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                                                                                                                    >  1:287318/1‑151 (MQ=255)                                         tttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGc                                                                                                            >  1:449103/1‑151 (MQ=255)                                              tcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgc                                                                                                       >  1:691383/1‑151 (MQ=255)                                              ccacccAGCCCGAAAATCCCGCTTGCGGTCGGGGCGCCCATTCCCCCCCCAATAATTGACGAACGCTCCCTGGTGCCCCCCGGGCTTG‑CCAGAAAACTAAACCCGGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCt                                                                                                           <  2:414086/143‑1 (MQ=255)                                               cagccagcCCGTGAATACCGCTTGCGGTCGTGGCGCCCCCTCCCCCCACCATAATTGACGACCGCTCACTGGTGCCCCCCGGGCTTG‑CCAGAAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgcg                                                                                                      <  1:498555/151‑1 (MQ=255)                                               cagccagcCCGTAAAAACCGTTTGGGTGCGTGGCCCCCACTTCGCCAACATTAATTGACGACCGCTCACTGGTGCCCCCGGGGTTTG‑CCAGGAAACTAACCCCCGGTTTGGTCCGTTTGACTATCATTCCCGCCTCTTCAATCGCTTcgc                                                                                                       <  1:888480/150‑1 (MQ=255)                                                  ccAGACGGGGAAACCCGCTTGCGGTCGTGGCGTCCCCTCCCCCCACCATAATGGACGAACGCTCACTGGTCCCCCCCGGGCTGG‑CCAGAAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa                                                                                                     <  2:753885/149‑1 (MQ=255)                                                       ccGTAAATCCCGTTGGGGGTCGTGGGGCCCCCTCCGCCCACCATAATTGCCGAACGCTCCCTGGTGCCCCCCGGGCTTG‑CCAGAAAACTTAACCCCGGTTTGGCCCGTTTGACTATCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCaa                                                                                              <  2:968520/151‑1 (MQ=255)                                                            aaTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATc                                                                                                              >  1:790674/1‑130 (MQ=255)                                                            aaTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCt                                                                                         >  1:1116725/1‑151 (MQ=255)                                                             atcccGCTTGCGGTCGTGGCCTCCACTCCCCCCCCCATAATTGACGAACGCTCACTGGTCCCCCCCGGGCTTG‑CCAGAAAACTAAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCtt                                                                                        <  1:856778/148‑1 (MQ=255)                                                                ccGCTTGCGGTCGTGGCGTCCCCTCCCCCCACCATAATTGACGACCGCTCACTGGTGCCCCCCGGGTTTG‑CCAGGAAACTAAACCCCGGTTTGGTCCGTTTGCCTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTcc                                                                                      <  2:1007255/150‑1 (MQ=255)                                                                    ttGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACt                                                                                 >  1:136675/1‑151 (MQ=255)                                                                    ttGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACt                                                                                 >  2:959870/1‑151 (MQ=255)                                                                    ttGCGGTCGTGGCGCCCCCTCCCCCCACCATAATTGACGACCGCTCACTGGTGCCCCCCGGGCTTG‑CCAGAAAACTAAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCcaca                                                                                   <  2:320678/149‑1 (MQ=255)                                                                     tGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACtt                                                                                >  2:354897/1‑151 (MQ=255)                                                                           cGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCATTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCAccc                                                                          >  1:1120205/1‑151 (MQ=255)                                                                            gTGGCGCCCCCTTCCCCCCCCATAATTGACGACCGCTCACTGGTCCCCCCCGGGTTTG‑CCAGGAACCTTAACCCGGGTTGGGCCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCa                                                                             <  1:830726/147‑1 (MQ=255)                                                                            gTGGCGCCCACTCCCCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGTTTG‑CCAGGAAACTAAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCct                                                                         <  1:779936/151‑1 (MQ=255)                                                                            gTGGCCCCCACTCCCCCCACCATAATTGACGAACGCTCACGGGTCCCCCCCGGGCTTG‑CCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCct                                                                         <  2:493100/151‑1 (MQ=255)                                                                              ggctccccctcccccccccATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAATTTAACCCCGGTTTGGCCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCctc                                                                        <  1:1063259/137‑1 (MQ=255)                                                                                   ccccTTCCCCCCCAATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGAAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCaa                                                                   <  1:1147915/147‑1 (MQ=255)                                                                                       ttCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGATTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAatca                                                                >  1:793730/1‑149 (MQ=255)                                                                                         ccccccccAAAATTGACGACCGCTCACTGGTCCCCCCCGGGCTTG‑CCAGAAAACTTAACCCCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAatcat                                                               <  2:364094/146‑1 (MQ=255)                                                                                                 aTAATTGACAAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGAAAACTTACCCCCGTTTTGTTCCTTTTTACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCa                                                    <  1:460625/151‑1 (MQ=255)                                                                                                   aaTTGACGAACGCTCACGGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCAt                                                   <  1:68130/150‑1 (MQ=255)                                                                                                       gACGAACGCTCACTGGTCCCCCCCGGGCTTG‑CCAGGAATCTTACCCCGGGTTTGGTCGGTTTTACTACCGGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCa                                                        <  1:829856/141‑1 (MQ=255)                                                                                                        aCGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTcc                                               >  1:248279/1‑149 (MQ=255)                                                                                                          gAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCTGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGt                                            >  1:1011103/1‑150 (MQ=255)                                                                                                             cccccACTGGTCCCCCCCGGCCTTG‑CCAGAAAACTTAACCCCGGTTTGTTCCTTTTGTCTACCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGcc                                        <  2:1010557/147‑1 (MQ=255)                                                                                                                tctgctgtccccctcGGCCTTG‑CGATGACACCTACCCCCGGTTTGGACCGTTTAACTACTCTTCCCGCCTCTTCAACCGCTTCGCGACGGGCAACATCTGCCACACTTTCCCCCTCTTCAATCAGCCCGGCAACCATCTCCAGTAGcct                                       <  2:833497/140‑2 (MQ=255)                                                                                                                tCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCa                                       >  2:1064428/1‑149 (MQ=255)                                                                                                                  aCTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAggg                                   >  1:1155102/1‑151 (MQ=255)                                                                                                                     ggTGCCCCCCGGGCTTC‑CCAGAAAACTTAACACCGTTTGGGCCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAgggg                                  <  1:1162127/149‑1 (MQ=255)                                                                                                                     ggTCCCCCCCGGGCTTG‑CCAGGAACCTAAACCCCGGTTTGGCCCTTTTAACTATCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAgggg                                  <  2:1170350/149‑1 (MQ=255)                                                                                                                        gCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCg                             >  1:268067/1‑151 (MQ=255)                                                                                                                         ccccccGGGCTTG‑CCAGAAAACTTAACCCCGTTTTGGCCCGTTGGCCTATAATCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCg                             <  1:816916/150‑1 (MQ=255)                                                                                                                          cccccGGGGTTGCCCAGAAAACTAAACCCCGTTTTGGCCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGc                            <  2:421069/151‑1 (MQ=255)                                                                                                                            cccGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTgg                         >  2:564969/1‑151 (MQ=255)                                                                                                                               gggTTTG‑CCAGGAAACTTAACCCCGTTTGGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGgtgt                      <  1:960188/151‑1 (MQ=255)                                                                                                                               gggCTTG‑CAAGAAACCTTAACCCCGGTTTGGCCCGTTAGCCTTCCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGgtgt                      <  1:442583/151‑1 (MQ=255)                                                                                                                                   ttG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGt                   >  1:1090693/1‑150 (MQ=255)                                                                                                                                       ccAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTAcgc               >  2:729494/1‑151 (MQ=255)                                                                                                                                           gAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGc            >  2:1051710/1‑150 (MQ=255)                                                                                                                                            aaaCTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCa           >  2:525375/1‑150 (MQ=255)                                                                                                                                            aaaCTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCa           >  1:928633/1‑150 (MQ=255)                                                                                                                                            aaaCTTAACACCGGTTTGGTCCGTTTGACTAACAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGc            <  1:212313/149‑1 (MQ=255)                                                                                                                                             aaCTAAACCCCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAAt         <  2:997030/151‑1 (MQ=255)                                                                                                                                                  aaCCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATc        <  1:206878/147‑1 (MQ=255)                                                                                                                                                  aaCCCCGGTTTGGCCCTTTTGATTCTCATTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCg      <  2:918862/149‑1 (MQ=255)                                                                                                                                                    cacCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATc  >  1:636603/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTG‑CCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATC  >  NC_000913/3177384‑3177680

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