Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,653	0	T	G	46.7%	72.2 / 32.0	60	P46P (CCA→CCC)	dppB	dipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (23/9);  new base G (0/28);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 1.72e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACA  >  NC_000913/3705504‑3705789                                                                                                                                                       |                                                                                                                                          aaGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTg                                                                                                                                          <  2:1139913/151‑1 (MQ=255)  aGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTgg                                                                                                                                         <  2:490787/151‑1 (MQ=255)     gcgGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGga                                                                                                                                      <  1:1017183/151‑1 (MQ=255)           cGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATccc                                                                                                                                <  1:620518/151‑1 (MQ=255)                   tCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGtt                                                                                                                           <  2:1176368/148‑1 (MQ=255)                    cccAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑cgcc                                                                                                                      >  2:349486/1‑151 (MQ=255)                     ccAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑cgccc                                                                                                                     <  1:976249/151‑1 (MQ=255)                       aaaCCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑cgcccgc                                                                                                                   >  2:498309/1‑151 (MQ=255)                           cGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATg                                                                                                               >  1:380156/1‑151 (MQ=255)                                 gCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGatcatc                                                                                                         >  2:724885/1‑151 (MQ=255)                                  cGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGatcatc                                                                                                         >  2:98414/1‑150 (MQ=255)                                     cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCg                                                                                                     >  2:530298/1‑151 (MQ=255)                                     cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCg                                                                                                     >  1:1155850/1‑151 (MQ=255)                                     cTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCg                                                                                                     >  1:308811/1‑151 (MQ=255)                                      ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCgg                                                                                                    >  1:1132139/1‑151 (MQ=255)                                           aTTAAAATCCCTACACCCCCATACATAACCCCCCAAATGTAATGGAGATACTGCTCCCACCTCTTTTTATCTAACCCGATTC‑CACCCAGCCCCTGCCCGTCACCCTCGGGGGAGTCCCCACTTC‑CCCCCGCCATGATCATCACCGGATc                                                                                                 <  1:125500/149‑1 (MQ=255)                                            ttaaaaTGCCTAGATCCCCATGCAAACCCCCCAAAATGTAAGGGGGATCCTGCTGCCACATGGTTTTATCTAAGCCGGTTC‑CAGCAAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTC‑CGCCCGCCATGATCATCACCGGATc                                                                                                 <  1:860923/145‑1 (MQ=255)                                              gAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCg                                                                                            >  1:373491/1‑151 (MQ=255)                                                aaTGCCGAGGCCGCCTGGCAAAACGCCCAAAATGTAAGGGATTTAGCGCTCCAACATGGTTTTATAT‑AGCCGGCTTCCAGCCGCCGGCTGGTCGTGGCGCGGGGGGGAGCCCCCACGTC‑CGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                          <  2:771307/151‑1 (MQ=255)                                                    cctatacccccATGCATACCCCCCAAAATGAAAGGGAGATCCTCCTCCCACATCGTTTTATCAAACCCGATTC‑CAGCCAGCAGCTGCGCGGGACGCTCTGGGGAGACCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                       <  1:887084/141‑1 (MQ=255)                                                    cctagtcccccATGCATACCCCCCAAAATTTAAGGGAGATCCTGCTCCCACTTCGTTTTATCAAACCCGAGTT‑CAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGATCCCACTTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGa                                                                                         <  1:920294/139‑1 (MQ=255)                                                    ccTAGACCGCCATGTATACCCCCCAAAATAAAAGGGGGATACGGCCGCCACCGTGTTTTATCAGAGCCGATTC‑CACCCCGCCGCTGCGGGTGCCGCTGGGGGGAGA‑CCCCCGTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                       <  2:79178/150‑1 (MQ=255)                                                         acccccATGCAAAACCCCCAAAATGTAGTGGAGATCCTGCTGCCCCATCTGTTTACCAAGCCCGAGTT‑CAGCCAGCCGCGGCGCGGGCCGCTCGGGGGAGCCCCCGCGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                     <  2:184294/143‑1 (MQ=255)                                                         aTCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTgg                                                                                 >  1:1115107/1‑151 (MQ=255)                                                           ggCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                               >  2:465755/2‑151 (MQ=255)                                                           cGCCTTGCATACCCCCCAAAATGTAAGGGAAATACTGCTCCCACATCGTTTTATCTAAGCCGAGTT‑CAGCCAGCAGCGGCGGGTGCCGCTCGGGGGAGATCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGa                                                                                <  1:158358/150‑1 (MQ=255)                                                           cGCCATGCAAACCCCCCAAAATTTAATGAAGATCCTGCTCCCACATCGTTTTATCTAAGCGGAGTT‑CAGCCAGCGGCGGCGCGTGACGCTCGGGGGAGACCCCTCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                               <  1:1003042/151‑1 (MQ=255)                                                           cGCCATCCATAACCCCCAAAATGTAGGGGAGATACTGCTCCCCCATGGTTTTATCTAAGCCGATTT‑CGGCCAGCGGGTGCGGGTGCCGCTCGGGGGAGCCCCCCCTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                               <  1:346322/151‑1 (MQ=255)                                                               ataaaaacccccaaaaaTTTAAGGGAGATCCGGCGGCCACATCGGTTTATCTACGCCGAGTT‑CAGCCCCCACCGGCGCGTGCCGCTGGGGGAAGCTCCCACGTC‑GGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                   <  2:149375/129‑1 (MQ=255)                                                                tGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAgg                                                                          >  2:523401/1‑151 (MQ=255)                                                                   ataccccccaAAATGTAATGGAGATACTGCTGCCCCGTCGTTTTATCTAAGCCGAGTT‑CAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaaa                                                                            <  1:779857/140‑1 (MQ=255)                                                                   aTACCGCCCAAAATGTAAGGGAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTT‑CAGCCGGCAGCTGCGCGTGCCGCTGGGGGGAGATCCCACTTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaa                                                                       <  1:553147/151‑1 (MQ=255)                                                                       cccccaAAATTTAAGGGAGATCCTGCTGCCACATGGTTTTATCAAAGCCGAGTT‑CGGCCAGCAGCTGCGCGTGCCGCTGGGGGGACACCCCGCGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaaa                                                                      <  1:807537/146‑1 (MQ=255)                                                                       cccccaAAATGAAAGGGATATACTGCCCCCACATCGTTTTATCTAAGCCGAGTC‑CAGCCGGCGGCGGCGCGGGACGCTCGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgt                                                                   <  1:259307/149‑1 (MQ=255)                                                                             aaaTGTAAGGGAGATTCTGCTGCCCCATCTTTTTATCTAACCCGAGTC‑CGGCCGGAGGGTGGGCGTGCCGCTGGGGGGAGATCCCAGGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                               <  2:868494/149‑1 (MQ=255)                                                                                 gaaagggggatcgtgccgcCACATGGGTTTATAAAAGCGGATTC‑CGCCCGGCGGCGGCGGGGGCCGTTGGGGGGAGACCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg                                                             <  2:863099/134‑1 (MQ=255)                                                                                     gggagatctgcctcccACGTCGTTTTACCTACCCGGAGTC‑CGCCCACCCGCTGCGCGTCACGCTGGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                     <  2:590021/137‑1 (MQ=255)                                                                                       gagttagcgctgcCACTTCGTTTTATCAAAGCCGAGTT‑CCGCCAGCAGCGGCGGGGGGCGCTGGGGGGAGACCCCACGTG‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                          <  2:1057387/136‑1 (MQ=255)                                                                                         gaTCCTGCGCCCACATCGTTTTATCAAAGCCGAGTT‑CGCCCGGCAGCTGCGGGTGTCGCTGGGGGGAGACCCCTTTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                   <  1:851578/149‑1 (MQ=255)                                                                                         gaTACTGCTGCCACATCGGTTTATCTAACCCGATTT‑CAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCCCGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGa                                                  <  1:513531/150‑1 (MQ=255)                                                                                           tccgcctcccACGTCGTTTTATCAAAGCCGAGTC‑CAGCCAGCAGCGGCGCGTGCCGCTCGGGGGAGACCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATa                                                <  1:966462/142‑1 (MQ=255)                                                                                             ctgctgCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAAc                                             >  2:942382/1‑151 (MQ=255)                                                                                               gctgcCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACg                                            >  2:366000/1‑150 (MQ=255)                                                                                                ctgcCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                          >  2:785567/1‑151 (MQ=255)                                                                                                      caTGGGTTTATCTAAGCGGTGTT‑CACCCGGCGGCTGCGGGTCCCGCTGGGGGGAGA‑CCCCCGTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATg                                    <  1:980530/151‑1 (MQ=255)                                                                                                         cGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                 >  1:190390/1‑151 (MQ=255)                                                                                                          ggTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                 >  2:874318/1‑150 (MQ=255)                                                                                                            tttATCAAAGCGGTGTT‑CAGCCAGCCGCGGCGCGTGCCGCTCGGGGGAGCCCCCACTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAg                               <  1:863053/150‑1 (MQ=255)                                                                                                               aTCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTccc                           >  1:673885/1‑151 (MQ=255)                                                                                                                         aGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                  >  1:461655/1‑150 (MQ=255)                                                                                                                         aGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                  >  2:259118/1‑150 (MQ=255)                                                                                                                           tt‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgc                  >  2:896219/1‑147 (MQ=255)                                                                                                                            tccagccagcGGCTGGGCGGGCCGCTCGGGGGAGACCCCCCTTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgag                <  1:221499/148‑1 (MQ=255)                                                                                                                            tccagccagcCGCGGGGCGTGCCGCTGGGGGGATCCCCCACGTG‑CCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                 <  2:959227/147‑1 (MQ=255)                                                                                                                              cagccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt             >  2:810365/1‑151 (MQ=255)                                                                                                                              cacccgccggcTGCGCGTGACGCGCGGGACATCCCCCACTGG‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                 <  2:324173/138‑1 (MQ=255)                                                                                                                                 ccagcCGCGGGGCGTGCCGCTGGGGGGAGCCCCCCCTTC‑GGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATa            <  2:866491/149‑1 (MQ=255)                                                                                                                                  cagcggCTGCGCGTCACGCTGTGGGGACACCCAAGGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt             <  2:1167721/147‑1 (MQ=255)                                                                                                                                          gcggGTGCCGCTGTGGGGAGCCCCCCCGTC‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAAc   <  1:610333/149‑1 (MQ=255)                                                                                                                                            gggTGCCGCTGGGGGGAGCCCCCCCCTG‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACa  <  2:997221/146‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                          AAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTT‑CAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTT‑CGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACA  >  NC_000913/3705504‑3705789

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