Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,812,792	0	C	A	43.3%	90.3 / 31.7	60	T21N (ACC→AAC)	dfp	fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (9/25);  new base A (26/0);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.04e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGA  >  NC_000913/3812651‑3812939                                                                                                                                               |                                                                                                                                                     aaTCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACt                                                                                                                                             <  1:909062/151‑1 (MQ=255)      cAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTgg                                                                                                                                           <  1:1014632/148‑1 (MQ=255)      cAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTg                                                                                                                                            <  2:619857/147‑1 (MQ=255)       aGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGGCTATAAAACCCCCTAAATGGGGcg                                                                                                                                       >  1:733202/1‑151 (MQ=255)        ggTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGcgt                                                                                                                                      <  1:213012/151‑1 (MQ=255)            aaGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAACTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGATTT‑GTGCCTTAAAAAAACCTGGAccagggggtcg                                                                                                                                    >  1:458931/1‑139 (MQ=255)             aGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGttt                                                                                                                                 <  2:36642/151‑1 (MQ=255)               gCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCGCCCTACAGGAAAAATCATCATGAGCCTGGCCGTTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTgc                                                                                                                               <  1:64632/151‑1 (MQ=255)               gCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGATTG‑CTGCCCATAAAAACCCTGAAATGGTGCGTCGttt                                                                                                                                 >  2:527203/1‑149 (MQ=255)                cACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAAACCCCGGACACGGGCGGCGttt                                                                                                                                 >  1:241088/1‑148 (MQ=255)                 aCGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAAACCCTGAACTGGTGCGTCGgtt                                                                                                                                 >  2:828483/1‑147 (MQ=255)                 aCGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAAACCCCGAAATTGTGGGTCGTTTgcg                                                                                                                              >  2:658288/1‑150 (MQ=255)                    acaacaAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTTCCTATAAAAACCCCGAACTTGTGCGGCGGTTGCGCGa                                                                                                                           >  2:958526/1‑150 (MQ=255)                     caacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAAACCCCCAAATGGTGCGTCGTTTGCGCGATc                                                                                                                         >  2:811210/1‑151 (MQ=255)                         aaaaGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAAACCCCGCACTGGGGCGTCGGTTTCGCGATcgcg                                                                                                                      >  2:1158345/1‑150 (MQ=255)                         aaaaGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCgg                                                                                                                     <  2:660938/151‑1 (MQ=255)                                      tAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAAACCCTGAACTGGGGCGTCGTTTTCGCGGTCGCGGGGCCGACGGcc                                                                                                          >  2:927547/1‑149 (MQ=255)                                         aaTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTa                                                                                                     <  2:628279/151‑1 (MQ=255)                                           tCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCCTTAAAAACCCTGAAACGGTGCGGCCTTTTCGCGAGCGCGGGGGCGACGGCcccgga                                                                                                     >  2:399546/1‑147 (MQ=255)                                             gCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑GCGCCTATAAAAACCCTGAAATGGGGTGTCGTTTGCGCGCTCGCCGGGGCGtcctcc                                                                                                          >  1:12535/1‑142 (MQ=255)                                              ccAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAAACCCTGTACTGGTGCGGCGTTTGGGCGGTCGCGGGGCCGACGTCCGCGTAGcc                                                                                                  >  2:355578/1‑149 (MQ=255)                                              ccAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCAt                                                                                                <  2:461576/151‑1 (MQ=255)                                              ccAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCAGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAAACCCCGAAATGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGcc                                                                                                  >  2:958384/1‑149 (MQ=255)                                               cAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAAACCCTGAACTGGTGGGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATg                                                                                               >  2:145650/1‑151 (MQ=255)                                                aaCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGTTTGCCT‑TAAAAACCCCGAAATGGGGCGGCGTTTTCGGGATCGCGGGGGCGGCGGCCGCCGAGcc                                                                                                  >  1:621547/1‑147 (MQ=255)                                                 aCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑TTGCCTCTAAAAACCCCGAACTGGGGCGGCGGTTTGGCGATCGCGGGGGCGACGGCCCCCGAGCCAtta                                                                                              >  1:903559/1‑148 (MQ=255)                                                  cTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTTCCTATAAAAACCCCGAAATGGGGCGTCGGTTGCGCGAACGCGGGGGCGGCGGCCCCGGAGGCACGAcc                                                                                            >  2:554195/1‑151 (MQ=255)                                                       ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGa                                                                                          <  1:750187/148‑1 (MQ=255)                                                            cACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCa                                                                                   <  2:440070/150‑1 (MQ=255)                                                             aCCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAAACCCCTAACTGGTGCGGCGTTTGTGCGATCGCGGGGCCGACGGCCCCGTAGCCCTGACCGGAGCGGcga                                                                                  >  1:99444/1‑148 (MQ=255)                                                             aCCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCaaa                                                                                 <  1:391615/151‑1 (MQ=255)                                                                  aCAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCt                                                                            <  2:313596/151‑1 (MQ=255)                                                                   cAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCCGAACTGGGGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGGAGCCATGACCGAAGCGGCCAAAGccct                                                                           >  2:1148082/1‑149 (MQ=255)                                                                      gAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTCTAAAAACCCTGAACTGGTGCGGCGTTTGCGCGATCGCGGGGGCGACGGCCGCGTAGACCTGACCGAAGCGGGAAAAACCctttt                                                                        >  2:1075276/1‑149 (MQ=255)                                                                      gAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCAGCACTTGGGGGTCGGTTGCGCGGTTTCGGGGGCGACGTCCCCGTAGCCCTGACCGAAGCGGGCAAAGCCctt                                                                          >  1:853197/1‑149 (MQ=255)                                                                         aaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAAACCCTCAACTGGGGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGCAGCGGCAAAAACCTTTAt                                                                        >  1:870280/1‑148 (MQ=255)                                                                                catGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACtt                                                              <  2:390104/151‑1 (MQ=255)                                                                                       cTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTTCGCGATCGCGGGGGCGACGTCCGCGGAGCCCTGACCGAAGCGGCAAAAACCCTTTTCACCCCACCTAGCtt                                                         >  2:1007472/1‑149 (MQ=255)                                                                                       cTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTT‑CTGCCCATAAAAACCCCGAAATTGGGGGGCGTTTTGGCGAGCGCGGGGCCGACCGCCGCGGAGCCCTGACCGAAGGGGCAAAAAACTTTTTCACCCCCCCTTGCtt                                                         >  2:641778/1‑149 (MQ=255)                                                                                         ggCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑GTGCCTATAAAAACCCCCGACTGGGGGGTGGGTTTCGCGAACGCGGGGCCGCCGCCCGCGGAGCCATGACCGAAGGGGCAAAAACCTTTTTCACCCCACCTAGCtt                                                         >  2:565886/1‑147 (MQ=255)                                                                                         ggCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCa                                                      <  1:598697/150‑1 (MQ=255)                                                                                          gCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTCTAAAAACCCCGAACCGGGGCGTCGTTTGGGCGAACGCGGGGCGGACGGCCGCGCCGCCATGACCGAAGCGGGGGAAGCCCTTTTCAcccccccttgctcttacg                                                    >  1:572104/1‑137 (MQ=255)                                                                                          gCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGGCCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAgg                                                    <  1:134835/151‑1 (MQ=255)                                                                                           ccGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGc                                                   <  1:1199448/151‑1 (MQ=255)                                                                                             ggTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑gg                                                <  1:305641/151‑1 (MQ=255)                                                                                                aaaaaaaTCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAAACCCTGAACCGGTGCGGCGTTTTCGCGAACGCGGGGCCGACGTCCGCGTAGGCATGACCGAAGCGGCAAAAGCCTTTTTTACCCCCCCTAGCTTGCAGGC‑GGgtt                                             >  2:956214/1‑151 (MQ=255)                                                                                                     aaTCGTTCTCGGCGTTAGCGGCGGTATTT‑TTGTCCTTAAAAACCCCGCACCGGGGCGGCGGTTGGGCGGTCGCGGGGGCGCCGTCCCCGGAGCCCCTACCGGAAGGGCAAAAACCCTTATTCCCCCCCCTAGCTTTCAGgg‑ggggttt                                           >  1:62156/1‑143 (MQ=255)                                                                                                      aTCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGtt                                       <  2:982637/151‑1 (MQ=255)                                                                                                        cGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTAt                                     <  1:18061/151‑1 (MQ=255)                                                                                                              cGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAAACCCCGAACTGGGGCGTCGGTTGCGCGGGCGCGGGGGCGACGGCCGCGGAGCCATGACCGAAGCGGCAAAAGACCTTTTTACCCCCCCTTGCCTGCCGGCGGGGTTTTGGTTTTccc                                  >  1:406950/1‑149 (MQ=255)                                                                                                               ggCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGttt                              <  2:682345/151‑1 (MQ=255)                                                                                                                 cGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTcc                            <  2:733227/151‑1 (MQ=255)                                                                                                                     aGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGAc                         <  1:834686/150‑1 (MQ=255)                                                                                                                        ggcggTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTc                     >  2:827558/1‑151 (MQ=255)                                                                                                                              aTTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGGTTCTGGTTATCCCGTTTCCcaca                        >  1:720777/1‑142 (MQ=255)                                                                                                                                    tgcCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGAccc           >  2:803450/1‑150 (MQ=255)                                                                                                                                    tgcCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGAccc           >  1:1028179/1‑150 (MQ=255)                                                                                                                                      ccTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGc        <  2:1028923/151‑1 (MQ=255)                                                                                                                                           aaaaCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGTC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGcc    <  1:827423/150‑1 (MQ=255)                                                                                                                                             aaCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGa  >  2:206836/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                     AATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTG‑CTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGA  >  NC_000913/3812651‑3812939

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