BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,694,023	(A)_5→4	coding (241/1374 nt)	uidB ←	glucuronide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,694,019	0	A	.	99.2%	537.8 / ‑3.1	120	coding (245/1374 nt)	uidB	glucuronide transporter
Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (49/70); total (50/70)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGTAAAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGTAAC  >  NC_000913/1693921‑1694114

tCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGTaa                                                                                                <  1:1731396/100‑1 (MQ=255)

   aaTATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGTaaaac                                                                                             <  2:168057/100‑2 (MQ=255)

     tatGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGTaaaacgg                                                                                           >  2:2115303/1‑97 (MQ=255)

        gCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCg                                                                                       >  1:1037426/1‑100 (MQ=255)

        gCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCg                                                                                       <  1:2737399/100‑1 (MQ=255)

        gCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCg                                                                                       <  2:937359/100‑1 (MQ=255)

         cATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCgg                                                                                      >  2:2519820/1‑100 (MQ=255)

         cATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCgg                                                                                      >  2:882704/1‑100 (MQ=255)

             caccacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAAtt                                                                                  >  1:654115/1‑100 (MQ=255)

             caccacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAAtt                                                                                  <  2:46687/100‑1 (MQ=255)

              accacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAAttt                                                                                 >  1:2881158/1‑100 (MQ=255)

               ccacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAAtttt                                                                                >  2:2204788/1‑100 (MQ=255)

                cacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTc                                                                               <  1:1894259/100‑1 (MQ=255)

                cacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTc                                                                               <  2:1337910/100‑1 (MQ=255)

                cacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTc                                                                               <  2:1443464/100‑1 (MQ=255)

                cacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTc                                                                               <  2:2478711/100‑1 (MQ=255)

                cacTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTc                                                                               <  2:724388/100‑1 (MQ=255)

                 acTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTcc                                                                              <  1:197728/100‑1 (MQ=255)

                 acTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTcc                                                                              <  2:2962671/100‑1 (MQ=255)

                  cTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTccc                                                                             <  2:400220/100‑1 (MQ=255)

                   tttGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTcccc                                                                            <  1:2426497/100‑1 (MQ=255)

                      gCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGc                                                                         >  2:655591/1‑100 (MQ=255)

                       cTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCg                                                                        <  2:780779/100‑1 (MQ=255)

                        tACCATGGCTCCAGTCGGTCAGCGCCCAGAATACCAGGACGGTGAAGATCACGCACGGCGCCAGCCCGACGAGT‑AAAACGCGCTGAATTTTCACCAGCgg                                                                       <  2:1491694/100‑1 (MQ=255)

                         aCCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCggg                                                                      <  1:714652/100‑1 (MQ=255)

                         aCCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCggg                                                                      <  1:1130671/100‑1 (MQ=255)

                          ccATTGCTCCAGTCGGTCAACACCCAGAATACCAGCACTCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AGAACGGGCGGATTTTTCCCCAGCggg                                                                      <  2:239428/99‑1 (MQ=255)

                            aTGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTAt                                                                   <  1:1152793/100‑1 (MQ=255)

                             tGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTAtt                                                                  <  1:2204807/100‑1 (MQ=255)

                              ggCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTc                                                                 >  2:2389541/1‑100 (MQ=255)

                               gCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTca                                                                <  2:2756698/100‑1 (MQ=255)

                               gCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTca                                                                >  1:2789437/1‑100 (MQ=255)

                                cTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGATATTcac                                                               <  1:1600500/100‑1 (MQ=255)

                                cTCCAGTAGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTcac                                                               >  1:1598898/1‑100 (MQ=255)

                                  ccAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCacac                                                             >  2:2333385/1‑100 (MQ=255)

                                    aGTCGGTCAGCACCCAGAATACCAGCACGCTGAGGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTg                                                           <  2:161156/100‑1 (MQ=255)

                                    aGTCGGTCAGCACCCAGAATACCAGCACGCTGAGGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGAATTCACACTg                                                           >  2:1997717/1‑100 (MQ=255)

                                       cGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTcc                                                        <  1:278856/100‑1 (MQ=255)

                                        ggTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTcca                                                       <  1:307437/100‑1 (MQ=255)

                                        ggTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTcca                                                       >  1:19781/1‑100 (MQ=255)

                                         gTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTccac                                                      >  1:2158593/1‑100 (MQ=255)

                                         gTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTccac                                                      <  1:2625031/100‑1 (MQ=255)

                                          tCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTccacc                                                     <  2:79195/100‑1 (MQ=255)

                                          tCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTccacc                                                     <  1:2943184/100‑1 (MQ=255)

                                           cAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTccacca                                                    <  1:1262849/100‑1 (MQ=255)

                                           cAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTccacca                                                    >  1:388557/1‑100 (MQ=255)

                                             gCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACt                                                  <  1:1820355/100‑1 (MQ=255)

                                              cACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTc                                                 >  1:175186/1‑100 (MQ=255)

                                              cACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTc                                                 >  1:2232514/1‑100 (MQ=255)

                                               aCCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCg                                                <  2:2881132/100‑1 (MQ=255)

                                               aCCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCg                                                >  2:1434103/1‑100 (MQ=255)

                                               aCCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCg                                                <  1:292426/100‑1 (MQ=255)

                                                cccAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGt                                               >  2:892333/1‑100 (MQ=255)

                                                  cAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTcc                                             <  2:2510892/100‑1 (MQ=255)

                                                  cAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTcc                                             <  2:1195031/100‑1 (MQ=255)

                                                  cAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTcc                                             <  2:2987311/100‑1 (MQ=255)

                                                   aGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCg                                            >  1:2032327/1‑100 (MQ=255)

                                                   aGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCg                                            >  1:741613/1‑100 (MQ=255)

                                                    gAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCgg                                           <  2:1024838/100‑1 (MQ=255)

                                                     aagaCCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGc                                          <  2:605952/97‑1 (MQ=255)

                                                     aaTACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCgg                                           <  2:654108/99‑1 (MQ=255)

                                                     aaTACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGc                                          <  2:1037416/100‑1 (MQ=255)

                                                      aTACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCa                                         <  2:2938012/100‑1 (MQ=255)

                                                        aCCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCaaa                                       >  2:2987383/1‑100 (MQ=255)

                                                        aCCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCaaa                                       <  1:231088/100‑1 (MQ=255)

                                                        aCCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCaaa                                       >  1:632408/1‑100 (MQ=255)

                                                         ccAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCGGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAg                                      <  2:1618917/100‑1 (MQ=255)

                                                         ccAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCaaan                                      <  2:1543750/100‑2 (MQ=255)

                                                         ccAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAg                                      <  2:2621228/100‑1 (MQ=255)

                                                         ccAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAg                                      <  1:2923051/100‑1 (MQ=255)

                                                         ccAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAg                                      <  2:3012238/100‑1 (MQ=255)

                                                          cAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGa                                     <  1:657294/100‑1 (MQ=255)

                                                           aGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGAc                                    >  1:1939091/1‑100 (MQ=255)

                                                               cGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCg                                >  1:1701352/1‑100 (MQ=255)

                                                               cGCTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCg                                >  2:1943535/1‑100 (MQ=255)

                                                                cgatctgATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCgg                               <  2:244705/94‑1 (MQ=255)

                                                                 cTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGc                              <  2:19780/100‑1 (MQ=255)

                                                                 cTGAAGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGc                              >  1:2514066/1‑100 (MQ=255)

                                                                     aGATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAg                          >  1:2472625/1‑100 (MQ=255)

                                                                      gATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCggagatcg                         >  1:244708/1‑94 (MQ=255)

                                                                      gATCATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAgg                         <  1:2831012/100‑1 (MQ=255)

                                                                         cATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCGCACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCAt                      <  2:1585413/100‑1 (MQ=255)

                                                                         cATTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCAt                      >  2:751791/1‑100 (MQ=255)

                                                                          aTTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATc                     >  2:2203417/1‑100 (MQ=255)

                                                                          aTTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATc                     >  2:2803532/1‑100 (MQ=255)

                                                                          aTTAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATc                     >  1:1358833/1‑100 (MQ=255)

                                                                           ttAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCg                    <  2:1089669/100‑1 (MQ=255)

                                                                           ttAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCg                    <  1:2950317/100‑1 (MQ=255)

                                                                           ttAACGGCGCAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCg                    >  1:2274401/1‑100 (MQ=255)

                                                                                  cgcAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATAccc             <  1:1741463/100‑1 (MQ=255)

                                                                                   gcAGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCg            >  2:3023215/1‑100 (MQ=255)

                                                                                     aGTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCa          <  1:2620369/100‑1 (MQ=255)

                                                                                      gTACCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCAc         >  2:1348251/1‑100 (MQ=255)

                                                                                        aCCGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCa       >  1:742373/1‑100 (MQ=255)

                                                                                          cGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGt     >  2:1276337/1‑100 (MQ=255)

                                                                                          cGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGt     <  1:1234951/100‑1 (MQ=255)

                                                                                          cGAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACCGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGt     <  2:291957/100‑1 (MQ=255)

                                                                                           gAAGAGT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGTa    <  1:892341/100‑1 (MQ=255)

                                                                                             agagT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGTAAc  <  1:995267/100‑1 (MQ=255)

                                                                                             agagT‑AAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGTAAc  <  2:2032309/100‑1 (MQ=255)

TCAAATATGCATACACCACTTTGCTACCATGGCTCCAGTCGGTCAGCACCCAGAATACCAGCACGCTGAAGATCATTAACGGCGCAGTACCGAAGAGTAAAAACGGGCGGAATTTTCCCCAGCGGGTATTCACACTGTCCACCACTCGTCCGGCAAAGACGTCGGCGAAGGCATCGAATACCCGCACCAGTAAC  >  NC_000913/1693921‑1694114

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