BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,025,492	A→C	5.6%	N341T (AAC→ACC)	helD →	DNA helicase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,025,492	0	A	C	5.6%	338.3 / 8.3	105	N341T (AAC→ACC)	helD	DNA helicase IV
Reads supporting (aligned to +/- strand): ref base A (50/49); new base C (4/2); total (54/51)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 109 total aligned reads displayed.

TAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTAAACGTATCAAG  >  NC_000913/1025345‑1025640

tAGTCAACAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACt                                                                                                                                                    <  2:1013995/150‑1 (MQ=255)

     aGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTgg                                                                                                                                               >  1:479605/1‑150 (MQ=255)

     aGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTgg                                                                                                                                               >  1:454232/1‑150 (MQ=255)

     aGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTgg                                                                                                                                               <  2:1098286/150‑1 (MQ=255)

      gCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTggg                                                                                                                                              >  1:557702/1‑150 (MQ=255)

      ccAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTggg                                                                                                                                              <  2:702141/149‑1 (MQ=255)

       cAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGga                                                                                                                                             >  1:953443/1‑150 (MQ=255)

       cAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGga                                                                                                                                             >  1:1305097/1‑150 (MQ=255)

           cTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGgatgat                                                                                                                                         >  1:1384678/1‑150 (MQ=255)

           cTGGAAAATGATACCGCTGCCAGTCATGAACTCTTTATTCCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGgatgat                                                                                                                                         >  2:138359/1‑150 (MQ=255)

             ggAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGAtgatga                                                                                                                                       >  1:1481150/1‑150 (MQ=255)

                    gATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTCGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGGAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTa                                                                                                                                <  2:1082928/150‑1 (MQ=255)

                              cccGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGcc                                                                                                                      <  2:535810/150‑1 (MQ=255)

                              cccGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGcc                                                                                                                      >  1:358812/1‑150 (MQ=255)

                               ccGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCt                                                                                                                     >  2:809411/1‑150 (MQ=255)

                                  tCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGc                                                                                                                  <  1:797534/150‑1 (MQ=255)

                                  tCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGc                                                                                                                  >  1:1414853/1‑150 (MQ=255)

                                   cATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGcc                                                                                                                 <  1:601589/150‑1 (MQ=255)

                                   cATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGcc                                                                                                                 <  1:560595/150‑1 (MQ=255)

                                   cATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGcc                                                                                                                 <  2:477615/150‑1 (MQ=255)

                                   cATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGcc                                                                                                                 <  2:548238/150‑1 (MQ=255)

                                       aaCTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCtctc                                                                                                             <  1:1218003/150‑1 (MQ=255)

                                        aCTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCg                                                                                                            <  2:401115/150‑1 (MQ=255)

                                        aCTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCg                                                                                                            <  2:521660/150‑1 (MQ=255)

                                          tctTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGcc                                                                                                          <  1:1171378/150‑1 (MQ=255)

                                           ctTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCt                                                                                                         >  1:567629/1‑150 (MQ=255)

                                           ctTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCt                                                                                                         >  2:1318464/1‑150 (MQ=255)

                                           ctTTATTGCTGAGGGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGGGGGCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCt                                                                                                         <  2:335641/150‑1 (MQ=255)

                                            tttATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTc                                                                                                        <  2:1414853/150‑1 (MQ=255)

                                              tATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGa                                                                                                      >  1:1178821/1‑150 (MQ=255)

                                                ttGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATc                                                                                                    >  2:94511/1‑150 (MQ=255)

                                                  gCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGc                                                                                                  >  1:673112/1‑150 (MQ=255)

                                                  gCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGc                                                                                                  >  1:1030621/1‑150 (MQ=255)

                                                  gCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGc                                                                                                  >  2:44876/1‑150 (MQ=255)

                                                   cTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCt                                                                                                 >  1:85173/1‑150 (MQ=255)

                                                    tGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTg                                                                                                <  2:557702/150‑1 (MQ=255)

                                                     gAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCCGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTc                                                                                                                                                  <  2:1481638/99‑1 (MQ=255)

                                                     gAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCCGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTc                                                                                                                                                  >  1:1481638/1‑99 (MQ=255)

                                                      aGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAATGTAACTTCTGTGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTggg                                                                                              <  2:1053933/150‑1 (MQ=255)

                                                      aGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTggg                                                                                              <  2:651743/150‑1 (MQ=255)

                                                       gTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGt                                                                                             <  2:1461707/150‑1 (MQ=255)

                                                        tGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTa                                                                                            <  1:634005/150‑1 (MQ=255)

                                                        gggCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTa                                                                                            <  2:150247/149‑1 (MQ=255)

                                                         ggCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTaa                                                                                           >  1:740950/1‑150 (MQ=255)

                                                           cgcAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGt                                                                                         >  2:593242/1‑150 (MQ=255)

                                                            gcaagcaaTGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTc                                                                                        <  2:1029438/150‑1 (MQ=255)

                                                            gcaagcaaTGCAGCGAAAAGAAAGCGCAGGCGAAGGGCGGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTc                                                                                        <  2:511473/150‑1 (MQ=255)

                                                              aagcaaTGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTg                                                                                      <  2:641186/150‑1 (MQ=255)

                                                                   aTGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCg                                                                                 >  1:640488/1‑150 (MQ=255)

                                                                    tGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCgg                                                                                <  2:1305097/150‑1 (MQ=255)

                                                                         cGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCa                                                                           <  1:809411/150‑1 (MQ=255)

                                                                               gaaagCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGGCTGATGCGGATGCACggtgg                                                                     <  2:1258371/150‑1 (MQ=255)

                                                                                 aagCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGAGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGgtggtg                                                                   <  2:898295/150‑1 (MQ=255)

                                                                                    cgcAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGcac                                                                <  2:1030621/150‑1 (MQ=255)

                                                                                    cgcAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGcac                                                                >  2:343367/1‑150 (MQ=255)

                                                                                    cgcAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGcac                                                                >  1:147437/1‑150 (MQ=255)

                                                                                     gcAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGcaca                                                               >  2:566114/1‑150 (MQ=255)

                                                                                       aGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAgg                                                             <  1:118039/150‑1 (MQ=255)

                                                                                        ggCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGc                                                            <  2:567629/150‑1 (MQ=255)

                                                                                             aGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGaaa                                                       <  2:740950/150‑1 (MQ=255)

                                                                                               ggCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATg                                                     >  1:318550/1‑150 (MQ=255)

                                                                                                gCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACTGGCAGAAATGa                                                    >  1:907947/1‑150 (MQ=255)

                                                                                                   ggcggcAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTg                                                 <  2:640488/150‑1 (MQ=255)

                                                                                                     cggcAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGcc                                               <  2:953443/150‑1 (MQ=255)

                                                                                                       gcAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGc                                                            >  1:189281/1‑135 (MQ=255)

                                                                                                       gcAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGc                                                            <  2:189281/135‑1 (MQ=255)

                                                                                                       gcAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAg                                             <  2:1384678/150‑1 (MQ=255)

                                                                                                       gcAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAg                                             <  2:358812/150‑1 (MQ=255)

                                                                                                            ggCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCAc                                        >  1:106010/1‑150 (MQ=255)

                                                                                                             gCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCAcc                                       >  1:391389/1‑150 (MQ=255)

                                                                                                                gACGGAAGAAATGCAGTGGTCAGTGCCAGAATGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTg                                           <  2:471513/143‑1 (MQ=255)

                                                                                                                gACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTg                                           >  1:471513/1‑143 (MQ=255)

                                                                                                                gACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGACGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGa                                    >  1:1356868/1‑150 (MQ=255)

                                                                                                                  cGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAag                                  <  2:391389/150‑1 (MQ=255)

                                                                                                                    gaagaaATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAagag                                >  2:806758/1‑150 (MQ=255)

                                                                                                                    gaagaaATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAagag                                <  1:208532/150‑1 (MQ=255)

                                                                                                                          aTGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTcgc                          >  1:1529834/1‑150 (MQ=255)

                                                                                                                                 ggTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGt                   >  1:1011703/1‑150 (MQ=255)

                                                                                                                                  gTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCTAGTGCACCCGAAGAGATTCGCGATCTGtt                  >  1:1162761/1‑150 (MQ=255)

                                                                                                                                   tCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTTCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTc                 <  1:1374420/150‑1 (MQ=255)

                                                                                                                                   tCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTc                 <  2:1314488/150‑1 (MQ=255)

                                                                                                                                    tcatcccAGAAGTTACCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCa                >  2:1392792/6‑150 (MQ=255)

                                                                                                                                    cAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCa                >  2:918054/1‑150 (MQ=255)

                                                                                                                                    cAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGACTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCAGAAGAGATTCGCGATCTGTTCa                >  2:1314837/1‑150 (MQ=255)

                                                                                                                                     catcccAGAAGTTACCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAg               >  1:671266/5‑150 (MQ=255)

                                                                                                                                     aGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAg               <  1:1318464/150‑1 (MQ=255)

                                                                                                                                     aGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAg               >  2:512430/1‑150 (MQ=255)

                                                                                                                                      atcccAGAAGTTACCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTc                 <  1:434183/144‑1 (MQ=255)

                                                                                                                                      atcccAGAAGTTACCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTc                 >  2:434183/4‑147 (MQ=255)

                                                                                                                                       tcccAGAAGTTACCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAg                                             <  2:1202238/116‑1 (MQ=255)

                                                                                                                                       tcccAGAAGTTACCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAg                                             >  1:1202238/3‑118 (MQ=255)

                                                                                                                                        gCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTaa            >  1:1071787/1‑150 (MQ=255)

                                                                                                                                           aGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTAAACg         <  2:1236930/150‑1 (MQ=255)

                                                                                                                                            gAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTAAACGt        >  1:499191/1‑150 (MQ=255)

                                                                                                                                              aGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTAAACGTAt      >  2:1032763/1‑150 (MQ=255)

                                                                                                                                              aGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCCCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACACGAAGAGATTCGCGATCTGTTCAGTAAACGTaa      >  2:539452/1‑149 (MQ=255)

                                                                                                                                                  aaCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGACGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTAAACGTATCAAg  >  1:1056347/1‑150 (MQ=255)

                                                                                                                                                  aaCTTCTGGGATGATGAAAAATTACAGCGGCGCCTTGCCTCTCGCCTCGATCGCTGTGTAAGTCTGATGCGGATGCACGGTGGTGCACAGTCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTAAACGTATCAAg  <  2:1481150/150‑1 (MQ=255)

                                                                                                                                                   aCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTc                 >  1:611639/1‑134 (MQ=255)

                                                                                                                                                   aCTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTc                 <  2:611639/134‑1 (MQ=255)

TAGTCAGCAAACTGGAAAATGATACCGCTGCCCGTCATGAACTCTTTATTGCTGAGTGGCGCAAGCAATGCAGCGAAAAGAAAGCGCAGGCGAAGGGCTGGCGGCAATGGCTGACGGAAGAAATGCAGTGGTCAGTGCCAGAAGGTAACTTCTGGGATGATGAAAAATTACAGCGTCGCCTTGCCTCTCGCCTCGATCGCTGGGTAAGTCTGATGCGGATGCACGGTGGTGCACAGGCAGAAATGATTGCCAGTGCACCCGAAGAGATTCGCGATCTGTTCAGTAAACGTATCAAG  >  NC_000913/1025345‑1025640

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