Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,580,634	0	C	T	39.3%	99.2 / 20.0	62	E277K (GAA→AAA)	hsdS	specificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base C (5/32);  new base T (24/0);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 2.00e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAAT  >  NC_000913/4580489‑4580783                                                                                                                                                    |                                                                                                                                                           aaTCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑ttt                                                                                                                                                      <  2:632367/151‑1 (MQ=255)   tCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTAAt‑ttttt                                                                                                                                                    >  2:935431/1‑149 (MQ=255)    cGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCt                                                                                                                                                    <  1:560402/150‑1 (MQ=255)    cGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTg                                                                                                                                                   <  2:415799/151‑1 (MQ=255)      aGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTT‑TTTAC‑ttttt                                                                                                                                                    >  1:47307/1‑146 (MQ=255)          cGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACa                                                                                                                                              <  2:802293/149‑1 (MQ=255)           gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCCC‑TTTCTGA‑ACAtt                                                                                                                                            <  2:451071/150‑1 (MQ=255)             aTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTT‑TTTTC‑TTTTTTT‑ACATTCt                                                                                                                                          >  1:333208/1‑150 (MQ=255)                 gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAC‑TTTCTGA‑ACAATTTAG‑aa                                                                                                                                     >  2:882053/1‑150 (MQ=255)                 gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCAACTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑aaa                                                                                                                                    <  2:426956/151‑1 (MQ=255)                    tATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAT‑TTTAC‑TTTTTTA‑ACAATCTAG‑aaa                                                                                                                                    >  1:238586/1‑148 (MQ=255)                           atatAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCg                                                                                                                          <  2:585785/151‑1 (MQ=255)                           atatAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCg                                                                                                                          <  2:1030990/151‑1 (MQ=255)                             ataGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGTCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGt                                                                                                                         <  2:221634/150‑1 (MQ=255)                               aGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTA‑TTTAC‑TTTTTGA‑ACAATCTAA‑AAAACGAATATTGtttt                                                                                                                      >  1:102419/1‑151 (MQ=255)                               aGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGtttt                                                                                                                      <  2:105207/151‑1 (MQ=255)                                    aTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTAC‑TTTCTTT‑ACATTCTAG‑AAAACGGATATCGttttttttt                                                                                                                 >  1:299868/1‑146 (MQ=255)                                         gATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAC‑TTTTTTT‑ACATTCTAG‑AAACCGAATATTGGTTTTAttt                                                                                                                 >  1:179609/1‑144 (MQ=255)                                            gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGATTTTG‑TTC‑C‑TTTTTTA‑AAATTCTAGAAAACCGAAATTCGttttatattcattgcc                                                                                                          >  1:938480/1‑139 (MQ=255)                                            gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAC‑TTTTTTA‑ACATTTTTG‑AAAACGAATATCGTTTTTAttt                                                                                                                 >  1:394906/1‑141 (MQ=255)                                             ttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTAC‑TTTTTGT‑ACATTCTAT‑AAAACGAATATTGTTTTTATCT‑ACATGGcc                                                                                                        >  1:414386/1‑150 (MQ=255)                                               gtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGc                                                                                                     <  1:625472/151‑1 (MQ=255)                                                 aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAC‑TTTCTTA‑ACATTCTAG‑AAACCGAATATTGTTTTTATTT‑ACATTGcc                                                                                                        >  1:104649/1‑146 (MQ=255)                                                        tCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAA‑TTTCTTA‑ACATTTTAT‑AAAACCAATATCGTTTTTATTTAACATGGcagc                                                                                                      >  2:1058958/1‑139 (MQ=255)                                                          aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAAc                                                                                          <  2:474326/151‑1 (MQ=255)                                                          aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAAc                                                                                          <  2:747677/151‑1 (MQ=255)                                                                  acaAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCgt                                                                                  <  1:427109/151‑1 (MQ=255)                                                                   caAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTGTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCgta                                                                                 <  2:518015/151‑1 (MQ=255)                                                                     aaCACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTT‑TTTAC‑TTTTTTA‑ACATTTTAT‑AAAAAGAATATTGTTTTTTTTA‑ACATTGCCAGCACCGACAAAAACAATGCgt                                                                                  >  2:789766/1‑148 (MQ=255)                                                                      acacCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTAC‑TTTCTTA‑ACATTCTAG‑AAACCGAATATCGTTTTTTTTT‑ACATGGGCAGCACGGACAGAACTAATGCgta                                                                                 >  1:281090/1‑148 (MQ=255)                                                                        acCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTAtt                                                                            <  1:845035/151‑1 (MQ=255)                                                                           aacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGa                                                                         <  1:530397/151‑1 (MQ=255)                                                                            acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTAC‑TTTCTGT‑AAAATCTTT‑AAAACGAATAACGTTTTTATCT‑ACATGGCCCGCACGTACAAAAATCAAGCGTAGTATTTGAt                                                                        >  2:155958/1‑151 (MQ=255)                                                                             caaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTTC‑TTTTTGA‑ACATTCTAG‑AAAACGAATATTGTTTTGATCT‑ACATGGCCAGCACGGACAGAAATAAAGCGTAGTAtttg                                                                          >  1:1015484/1‑146 (MQ=255)                                                                              aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGa                                                                      <  2:333061/151‑1 (MQ=255)                                                                                     aaaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAaca                                                                <  1:568128/150‑1 (MQ=255)                                                                                          tCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAAcaccac                                                            <  1:812828/149‑1 (MQ=255)                                                                                           ccGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑GTTTG‑TTTCA‑TTTTTTA‑ACATTTTTA‑AAAACGAATATCGGTTTGATTT‑ACATTGGCAGCACGCCCAGAAAAAAAGCGTGGTATTTGGTGGACCACACCAC‑ttt                                                        >  2:944518/1‑151 (MQ=255)                                                                                            cGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTA‑TTTAC‑TTTTTTG‑ACATTTTAG‑AAACCGCATAACGTTTTTATTT‑ACATGGCCAGCACGTACAGAAATAATACGTAGGATTTGGTGAACAACAACCC‑ttt                                                        >  1:1059508/1‑150 (MQ=255)                                                                                              ttATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTAA‑TTTTTTA‑ACATTTTTG‑AAAACGAATATCGTTTTGATCT‑TCATTGCCACCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑att                                                        >  2:670192/1‑148 (MQ=255)                                                                                                ataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAC‑TTTTTTA‑ACATTCTTA‑AAACCCAATATCGTTTTTATCT‑ACATTGCCCGCACGTACAGAAAAAATGCGTAGGATTTGGTGACCAACACCAC‑TTTCCttt                                                   >  1:115432/1‑151 (MQ=255)                                                                                                  aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCCC‑CTTCTTA‑ACATTTTAG‑AAACCGAATATCCTTTTTATCT‑TAATGGGCAGCCAGGACAAGAATAATGGCGAGTATTGGATGGACAAAACCAC‑CTTTAttt                                                   >  1:386245/1‑149 (MQ=255)                                                                                                    cGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCt                                               <  1:681809/151‑1 (MQ=255)                                                                                                       gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAC‑TTTTTTA‑ACATTTTTA‑AAAACGAATATTGTTTTGATCT‑ACATTGCCAGCACGTACAGAACAAAAGCGTAGTATTGGATGACCAACACCAC‑CTTCATTTTGCttt                                             >  1:172186/1‑150 (MQ=255)                                                                                                       gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTtg                                            <  1:1039420/151‑1 (MQ=255)                                                                                                         aaataaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTtga                                           <  2:192846/150‑1 (MQ=255)                                                                                                            taaaAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAGTTTCAT‑TTTTTTACACATTATA‑‑AAACC‑AATATTGTTTTTTTTT‑ACATGGCCAGCACGTACAAAAATAATGAGGAGTATTGGATGATCAAAACCACATTTCTTTTGG‑TTtgat                                          >  2:1010173/1‑148 (MQ=255)                                                                                                             aaaaGATCTCCATCTTGTAATTTGTGGCGG‑GTTAG‑TTTCC‑TTTTTGA‑ACAATTTTA‑AAAACGAAAATTGTTTTTATCT‑ACATTGCCAGCACGTACAGGAATAAATCGTCGGATTTGGTTAACAACACCAC‑CTTTATTTGTCttt                                             >  1:172680/1‑144 (MQ=255)                                                                                                               aaGATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTAC‑TTTTTTA‑ACATTTTAG‑AAAACGAAAATCGTTTTTATCT‑ACATTGCCAGCACGGACAGAAATAAAGCGTAGTATTTGATGAACAACACCAC‑CTTTATTTTGCTTTTATGaaa                                      >  1:846204/1‑149 (MQ=255)                                                                                                                 gATCTCCATCTTGTAATTTGTGGCGG‑TTTTG‑TTTCC‑CTTTTTA‑CCATTTTTG‑AAAACCAATATCGGTTTTTTTT‑ACATTGCCCGCACGTACAAAAATAATACGTAGTATTGGATGGCCAACACCCC‑ATTCATTTTGCTTTGTTGaaaaa                                    >  2:1034671/1‑147 (MQ=255)                                                                                                                  aTCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTTAC‑TTTTTGA‑AAAATCTTG‑AAAACGAAAATCGTTTTGATTT‑ATATTGCCAGCACGGACAAAAATAAATCGGAGTATTTGATGAACAACACCAC‑CTTTATTTTGGTTTGATGAAAGaac                                  >  1:833276/1‑148 (MQ=255)                                                                                                                   tctcCATCTTGTAATTTGTGGCGG‑TGTAT‑TTTACTTTTTTTA‑CCCTTCTTA‑AAACC‑AATATCGGTTTTATTT‑ACATTTCCAGCACATCAAAAAATAATGAGTAGGATTTTGTGATCAACAACCA‑TTTTTTTTTGCTTTGTTGAAAGaa                                   >  2:18892/1‑147 (MQ=255)                                                                                                                      ccATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTAc                             <  2:500003/151‑1 (MQ=255)                                                                                                                      ccATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTAc                             <  1:234540/151‑1 (MQ=255)                                                                                                                            tGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGt                           <  1:601400/147‑1 (MQ=255)                                                                                                                                        cGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagatag           <  1:1059316/151‑1 (MQ=255)                                                                                                                                          g‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagatagat         <  1:230/151‑1 (MQ=255)                                                                                                                                              tAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTc       <  2:938171/150‑1 (MQ=255)                                                                                                                                               aG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGTACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCTC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCa      <  2:397938/150‑1 (MQ=255)                                                                                                                                                g‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCaa     >  1:111722/1‑150 (MQ=255)                                                                                                                                                   tCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCaaaa   <  2:278036/150‑1 (MQ=255)                                                                                                                                                   tCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCaaaa   <  1:694203/150‑1 (MQ=255)                                                                                                                                                    cAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAAt  <  2:545096/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                           AATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑TTTAG‑TTCAC‑TTTCTGA‑ACATTCTAG‑AAACCGAATATCGTTTTGATCT‑ACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCAC‑TTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAAT  >  NC_000913/4580489‑4580783

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