Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,840,646	0	T	C	45.5%	68.5 / 33.2	67	D365G (GAT→GGT)	nepI	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (26/10);  new base C (0/30);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 1.57e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGTTAGCCCC  >  NC_000913/3840496‑3840795                                                                                                                                                         |                                                                                                                                                      cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAAtattat                                                                                                                                                       <  2:1008909/151‑1 (MQ=255) cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAAtattat                                                                                                                                                       >  2:350460/1‑151 (MQ=255) cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAAtattat                                                                                                                                                       <  2:914274/151‑1 (MQ=255)   cGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATcc                                                                                                                                                     <  2:678570/151‑1 (MQ=255)    ggATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCa                                                                                                                                                    <  2:752037/151‑1 (MQ=255)       tatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAgcg                                                                                                                                                 >  2:1066593/1‑151 (MQ=255)            aaGAACGGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATaa                                                                                                                                            <  1:101443/151‑1 (MQ=255)              gAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAAcc                                                                                                                                          >  2:328949/1‑151 (MQ=255)                aCAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAAccgc                                                                                                                                        >  1:1030510/1‑151 (MQ=255)                 cAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAAccgcc                                                                                                                                       <  1:217030/151‑1 (MQ=255)                     aaaaacTACTGACCAACTCTTCTTTACCTGATATCAGGATTTCTTCATTTTCCCCTTTGCAGAAACCAAAAATGCAGTCAG‑CAACATCAATGTGCCGGACAACACCAACGGCGAAGTCACAC‑CAATATTATCCAGCGCATAACCGCCGa                                                                                                                                     <  1:5043/146‑1 (MQ=255)                        aacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGcc                                                                                                                                >  2:817573/1‑151 (MQ=255)                              tGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGcc                                                                                                                           >  2:1067270/1‑150 (MQ=255)                                   aaCGTTTCTTTACCTGATATCAGGATTTCTTCTTTTTCACCTTTGCAGTAACAAACAATGCAGTCAC‑CAACATCAATGTGCCGAAAAACACCAACGGCGAATTCAGAC‑CAATTTTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGt                                                                                                                     <  1:1079019/151‑1 (MQ=255)                                      gCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTAtt                                                                                                                  >  1:542743/1‑151 (MQ=255)                                        ttcttTCCCTGTTTTCAGGTTTTCTTCATTTCCCCCTTCGCAGTAACAAAAAATGCAGTCAA‑CAACATCAATGCGCCGAAAAACACCACCGGCAAAGACAGAC‑CAAAATTACCCAGCCCAAAACCGCCGATTGCCGCGCCACAGGTAt                                                                                                                   <  2:1104348/148‑1 (MQ=255)                                            ttACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCa                                                                                                              >  2:18270/1‑149 (MQ=255)                                             tACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAAGGCAGTAAG‑CAACATCAATGGGCCGGACAACACCAACGGCGAAGTCAGAC‑AAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAg                                                                                                            <  2:850544/150‑1 (MQ=255)                                               ccTGATATCAGGTTTTTTTCATTTTCACCTTTGCAGTAACAAACAATGCAGTCAG‑CAACATCAATGTGCCGAAAAACATCAACGGCGAAGTCAGAC‑AAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCt                                                                                                          <  2:542803/150‑1 (MQ=255)                                               ccTAATATAGGGATTTCTTAATTTTCCCCTTTGCAAAACCCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTg                                                                                                         <  2:693974/151‑1 (MQ=255)                                                   atatCAGTTTTTCTTCTTTTTCCCCTTTACAGAAACCAACAATGCAGTCAG‑CAACACCAATGTGCCGGAAAACACCAACGGGGAAGTCAGAC‑CATTATTACCCAGCACAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGc                                                                                                           <  1:817510/145‑1 (MQ=255)                                                   atatCAGGATTTCTTCATTTTCCCCTTTGCAGTAACCAACAATGCAGTCAG‑CACCATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGa                                                                                                        <  2:952405/148‑1 (MQ=255)                                                       cAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCg                                                                                                 >  2:276180/1‑151 (MQ=255)                                                            tttCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAATTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCAcc                                                                                            >  1:547301/1‑151 (MQ=255)                                                             ttctttATTTTCCCCTTTGGAAGAAACAAAAATGGAAGAAGAAAAAATAAATGTGCCGACCAACA‑CAACGGGGAAAGAAGACAAAATTTAACCCAGCCAATACCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCa                                                                                              <  1:179330/149‑1 (MQ=255)                                                                ttcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGa                                                                                         >  1:507189/1‑150 (MQ=255)                                                                ttcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGa                                                                                         >  2:1003804/1‑150 (MQ=255)                                                                 tcATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAAt                                                                                       >  2:813946/1‑151 (MQ=255)                                                                    ttttCCCTTTTCCAGAACCCAAAAGCGCAGTCAC‑CACCATCAATGTGCCGCACAACAGAATCGCCGAAGACAGCC‑CAATTTTCTACGGCCCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTg                                                                                          <  2:1019927/145‑1 (MQ=255)                                                                    ttttCCCCTTTGAATTAACAAAAAATGCAGTCAG‑CAACATCAATGTGCCGGCCAACACCAACGGCGAAGTCAGAC‑CAATATTACCCAGCGCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAg                                                                                     <  2:98546/150‑1 (MQ=255)                                                                    ttttCCCCTTGGCAGTAACCAACAATGCAGTCAG‑CAACATAAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑AAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTg                                                                                          <  1:410947/145‑1 (MQ=255)                                                                        cACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCg                                                                                >  1:9195/1‑151 (MQ=255)                                                                         aCCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGCATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCgg                                                                               >  1:679684/1‑151 (MQ=255)                                                                          ccTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTACCCAGCGCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCgg                                                                               <  2:1024400/150‑1 (MQ=255)                                                                          ccTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAAGGTGCCGGACAACACCACCGGCGAAGCCAGAC‑CAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGAc                                                                                   <  2:660367/146‑1 (MQ=255)                                                                          ccTTTGCAGTAACAAAAAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑AAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCgg                                                                               <  1:401205/150‑1 (MQ=255)                                                                            tttGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCtt                                                                            >  2:290189/1‑151 (MQ=255)                                                                                   aaaCCAACAATGCAGTCAG‑CAACATCAATGTGCCGAACAACACCAACGGCGAAGCCAGAC‑AAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGc                                                                     <  1:1067174/150‑1 (MQ=255)                                                                                    aaCCAAAAAAGGAGTCGC‑AAACATAATTGCGCCGAACAACAAAAAGGGCGAAGCCAGAC‑AAATATCACCCAGCGGATATCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCt                                                                       <  2:452263/148‑1 (MQ=255)                                                                                        aaaaaTGCAGTCAG‑CAACATAAATGTGCGGGACAACATCAACGGCGAAGTCAGAC‑AAATATTATCCGGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTg                                                                  <  2:84762/146‑1 (MQ=255)                                                                                         acaaTGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATc                                                               >  2:142106/1‑151 (MQ=255)                                                                                         acaaTGCAGTCAG‑CAACATCAATGTGCCGGACAACACCAACGGCGAAGCCAGAC‑CAATATTACCCAGCGCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATc                                                               <  1:34298/151‑1 (MQ=255)                                                                                          caaTGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCg                                                              >  2:908220/1‑151 (MQ=255)                                                                                                aGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCag                                                         >  2:650851/1‑150 (MQ=255)                                                                                                 agccaacAACATCAATGTGCCGCACAACACCACCGGCGAAGTAAGAC‑AAATATTACCCGCCCCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGcc                                                           <  1:1055345/142‑1 (MQ=255)                                                                                                       cAACATCAATGTGCCGAACAAAATCAACGGCGAAGTCAGAC‑AAATATTACCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCagcgagc                                                    <  2:426800/149‑1 (MQ=255)                                                                                                        aaaataaaTGTGCCGAAAAACATCACCGGCGAAGTAAAAC‑AAATATTACCCAGCGCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCagcgagcg                                                   <  1:35266/143‑1 (MQ=255)                                                                                                        aaCATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCggg                                                 >  2:552479/1‑151 (MQ=255)                                                                                                        aaCATCAATGTCCCGAACAAAAACAACGGCGAATTCAGAC‑AAATATAACCCAGCAAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCagcgagc                                                    <  1:695047/148‑1 (MQ=255)                                                                                                        aaCATCAATGGCCCGGACAACATCAACGGCGAAGTAAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCagcga                                                      <  2:736453/146‑1 (MQ=255)                                                                                                        aaCACCAATGTGCCGAACAACACCAACGGCAAAGTCAAAC‑AAAAATTACCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCgg                                                  <  1:10061/150‑1 (MQ=255)                                                                                                                  tgCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGt                                       >  1:1041497/1‑151 (MQ=255)                                                                                                                    ccGGACAACAACAAGGGCGAAGTCAGAC‑CAATATTAACAAGCGCATACCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTGCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGtt                                      <  1:618953/150‑1 (MQ=255)                                                                                                                      ggACAACATCAACGGCGAAGTCAGAC‑CAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTg                                     <  2:731292/149‑1 (MQ=255)                                                                                                                          aacaTCAAGGGCGAAGTCAGCA‑AAATATTACCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAg                                <  1:730499/150‑1 (MQ=255)                                                                                                                               cAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGa                            >  1:729033/1‑149 (MQ=255)                                                                                                                                accGGCGAAGTCAGAC‑CAATATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGcc                              <  2:774083/144‑1 (MQ=255)                                                                                                                                aaCGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCACGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCGAGCGAGCGGGTGATCCACGTTGACCAGCCGAc                           >  1:869713/1‑149 (MQ=255)                                                                                                                                aaCGGCAAAGTCAGAC‑CAATATTACCCAGCACAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACg                          <  2:819678/150‑1 (MQ=255)                                                                                                                                   ggggAAGTCAGCC‑CAATATTACCCAGCCAATACCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGa                       <  1:414288/147‑1 (MQ=255)                                                                                                                                   ggCGAAGCCAGAC‑CAATTTTACCCAGCGAATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGaa                      <  2:751944/151‑1 (MQ=255)                                                                                                                                      gAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAAcc                    >  2:424404/1‑150 (MQ=255)                                                                                                                                       aaGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCa                   >  1:221953/1‑150 (MQ=255)                                                                                                                                       aaGTAAGAC‑AAATTTTACCCAGCGCAAAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCaa                  <  2:388531/151‑1 (MQ=255)                                                                                                                                          tCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGc               >  1:261661/1‑151 (MQ=255)                                                                                                                                           cAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGc               >  1:1121380/1‑150 (MQ=255)                                                                                                                                           acaaccaaaTATTATCAGGCTAAAACCCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGc               <  1:76658/144‑1 (MQ=255)                                                                                                                                            aGAA‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATg                <  2:299092/148‑1 (MQ=255)                                                                                                                                                caaaTATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAg          <  1:1066497/149‑1 (MQ=255)                                                                                                                                                  aaTATTACCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGt         <  2:362407/150‑1 (MQ=255)                                                                                                                                                        atCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGTTAGcccc  >  1:1032046/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                      CCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAG‑CAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAC‑CAATATTATCCAGCGCATAACCGCCGATTGCCGCGCCACAGGTATTAGCAAGCTGAATAACCGCCACCTGAATAGACCCGGCTTTTTCTGCCTGATCGGCCAGCGAGCGGGTGATCCACGTTGACCAGCCGACGGGAACCAATGCAAAAGTTAGCCCC  >  NC_000913/3840496‑3840795

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