BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,676,601	T→C	T74A (ACC→GCC)	yhjH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,676,601	0	T	C	96.6%	537.0 / 5.5	178	T74A (ACC→GCC)	yhjH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/3); new base C (83/89); total (86/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

Only 100 of 184 total aligned reads displayed.

ATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACGGCCATTAACC  >  NC_000913/3676511‑3676690

aTGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTGTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGa                                                                                         <  1:3363000/93‑1 (MQ=255)

  gCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGAtt                                                                                       <  2:1790499/93‑1 (MQ=255)

  gCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGAtt                                                                                       <  1:4439890/93‑1 (MQ=255)

     agcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCa                                                                                    <  2:1942854/93‑1 (MQ=255)

     agcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCa                                                                                    >  1:4213312/1‑93 (MQ=255)

      gcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAg                                                                                   >  1:927564/1‑93 (MQ=255)

        agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTa                                                                                 <  1:2724178/93‑1 (MQ=255)

         ggCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTaa                                                                                >  2:3515756/1‑93 (MQ=255)

          gCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTaaa                                                                               >  2:2035953/1‑93 (MQ=255)

           ccGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTaaaa                                                                              <  1:4177104/93‑1 (MQ=255)

            cgggCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAAt                                                                             <  1:4060192/90‑1 (MQ=255)

            cGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAAt                                                                             <  2:1160965/93‑1 (MQ=255)

              tGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAg                                                                           >  1:4469960/1‑93 (MQ=255)

              tGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAg                                                                           <  1:1796041/93‑1 (MQ=255)

              tGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAg                                                                           >  1:3096662/1‑93 (MQ=255)

              tGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAg                                                                           >  2:4714642/1‑93 (MQ=255)

               gCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGc                                                                          <  2:170847/93‑1 (MQ=255)

               gCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGc                                                                          >  2:4154545/1‑93 (MQ=255)

               gCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGc                                                                          <  1:247132/93‑1 (MQ=255)

                ctctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCg                                                                         <  1:1060351/93‑1 (MQ=255)

                ctctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCg                                                                         <  2:3059494/93‑1 (MQ=255)

                 tctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGa                                                                        <  2:2761023/93‑1 (MQ=255)

                  ctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGAt                                                                       >  1:3830509/1‑93 (MQ=255)

                   tataAAGAGCTCGACTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATc                                                                      <  1:1855984/93‑1 (MQ=255)

                   tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATc                                                                      <  1:4317011/93‑1 (MQ=255)

                   tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATc                                                                      <  1:3213165/93‑1 (MQ=255)

                    ataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATcc                                                                     <  2:4581454/93‑1 (MQ=255)

                     taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCcg                                                                    <  2:3743046/93‑1 (MQ=255)

                     taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCcg                                                                    >  2:2367584/1‑93 (MQ=255)

                     taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCcg                                                                    >  1:2094296/1‑93 (MQ=255)

                      aaaGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCcgg                                                                   <  2:4135853/93‑1 (MQ=255)

                       aagaagTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCcggc                                                                  >  2:442971/1‑93 (MQ=255)

                        agaagTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCcggcg                                                                 >  2:1569273/1‑93 (MQ=255)

                         gaagTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCcggcgg                                                                >  2:2245947/1‑93 (MQ=255)

                            gTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAg                                                             >  1:1465022/1‑93 (MQ=255)

                            gTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAg                                                             >  1:979196/1‑93 (MQ=255)

                             tCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAgg                                                            >  1:275390/1‑93 (MQ=255)

                              cGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGc                                                           <  2:372584/93‑1 (MQ=255)

                              cGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGc                                                           <  1:218088/93‑1 (MQ=255)

                                 ttttttGCGCCAGCAAATCAATCTGCTCTTTCACAACCGCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGtt                                                        <  2:1731127/92‑1 (MQ=255)

                                 cTTTTTGCGCCGGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGtt                                                        <  2:4220949/93‑1 (MQ=255)

                                 cTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGtt                                                        >  1:2278990/1‑93 (MQ=255)

                                  tttttGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGttc                                                       >  1:3448117/1‑92 (MQ=255)

                                  tttttGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTg                                                       <  1:535388/93‑1 (MQ=255)

                                  tttttGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTg                                                       >  2:2161083/1‑93 (MQ=255)

                                  tttttGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTg                                                       >  1:2649912/1‑93 (MQ=255)

                                  tttttGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTg                                                       >  1:2699751/1‑93 (MQ=255)

                                     ttGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGa                                                    <  1:3323560/93‑1 (MQ=255)

                                     tggcgcCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCAATCCGGCGGCAGGCGTTGCGa                                                    <  1:4490781/91‑1 (MQ=255)

                                      tGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGaa                                                   <  1:1295175/93‑1 (MQ=255)

                                      tGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGaa                                                   <  1:2624895/93‑1 (MQ=255)

                                        cgcCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAgg                                                 <  1:3583032/93‑1 (MQ=255)

                                        cgcCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAgg                                                 <  1:3993792/93‑1 (MQ=255)

                                        cgcCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAgg                                                 <  1:2898458/93‑1 (MQ=255)

                                        cgcCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAgg                                                 <  2:3508396/93‑1 (MQ=255)

                                        cgcCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAgg                                                 <  2:2605378/93‑1 (MQ=255)

                                         gcCAGCAAATCAATCTGCGCTTTCACAGCCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAggg                                                <  1:4601245/93‑1 (MQ=255)

                                          ccAGCAAATCAATCTGCTCTTTCACACCCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                                               <  2:3632955/93‑1 (MQ=255)

                                          ccAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                                               <  2:2914429/93‑1 (MQ=255)

                                          ccAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                                               <  2:270462/93‑1 (MQ=255)

                                          ccAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                                               <  1:1733178/93‑1 (MQ=255)

                                             gCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCa                                            >  1:1017663/1‑93 (MQ=255)

                                             acAAATCAATCCGCTCTTTCCCAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCa                                            <  1:3758164/92‑1 (MQ=255)

                                               aaaTCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAg                                          <  1:1584715/93‑1 (MQ=255)

                                                aatcaatcTGCTCTTTCACAGCCACCATCCGATGGCTGACGGCGATTTAAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAgg                                         <  1:3396493/93‑1 (MQ=255)

                                                aatcaatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAgg                                         <  1:3710842/93‑1 (MQ=255)

                                                aatcaatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAgg                                         <  2:1358397/93‑1 (MQ=255)

                                                 atcaatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAggg                                        >  1:4315296/1‑93 (MQ=255)

                                                    aatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATg                                     >  1:802858/1‑93 (MQ=255)

                                                     atcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGc                                    >  1:4117601/1‑93 (MQ=255)

                                                     atcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGc                                    <  1:1153142/93‑1 (MQ=255)

                                                       cTGCTCTTTCACAACCTCCATCCGATGGCTGACGGCGATGTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGt                                  <  1:2824828/93‑1 (MQ=255)

                                                          ctctTTCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAGGGGTTCAAGGGATGCGTGac                               >  1:2370866/1‑93 (MQ=255)

                                                             tttCACAACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGaccac                            <  1:3782690/93‑1 (MQ=255)

                                                                cacaACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGt                         >  1:3699097/1‑93 (MQ=255)

                                                                cacaACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGt                         <  1:110118/93‑1 (MQ=255)

                                                                 acaACCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGtt                        >  1:2363660/1‑93 (MQ=255)

                                                                   aaCCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTaa                      >  1:2771902/1‑93 (MQ=255)

                                                                   aaCCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTaa                      >  1:1931501/1‑93 (MQ=255)

                                                                    aCCTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAAt                     <  2:4271252/93‑1 (MQ=255)

                                                                     ccTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATa                    >  2:144338/1‑93 (MQ=255)

                                                                     ccTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATa                    >  2:4561112/1‑93 (MQ=255)

                                                                     ccTCCATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATa                    >  2:3871470/1‑93 (MQ=255)

                                                                       cccATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGc                  <  2:299353/92‑1 (MQ=255)

                                                                         cATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTc                >  2:1740319/1‑93 (MQ=255)

                                                                         cATCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTc                >  1:3006055/1‑93 (MQ=255)

                                                                          aTCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTcc               >  2:3852742/1‑93 (MQ=255)

                                                                           tCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTccc              >  1:1552145/1‑92 (MQ=255)

                                                                           tCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCa              >  2:4228123/1‑93 (MQ=255)

                                                                           tCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCa              >  2:4613322/1‑93 (MQ=255)

                                                                           tCCGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCa              >  2:2390975/1‑93 (MQ=255)

                                                                             cGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACg            >  1:2842850/1‑93 (MQ=255)

                                                                             cGATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACg            <  2:840741/93‑1 (MQ=255)

                                                                              gATGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACgg           <  2:673840/93‑1 (MQ=255)

                                                                               gtGGCTGACGGCGATTTTAGTAAAATAGCGAGCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCCCGGc          <  2:4077551/92‑1 (MQ=255)

                                                                               aTGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACGGc          >  2:3157460/1‑93 (MQ=255)

                                                                               aTGGCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACGGc          <  1:3726497/93‑1 (MQ=255)

                                                                                tGGCTGAAGGCGATTTCAGTAAAATAGCGATCAGGCGGCAGGCGGTGCGAAGGGTTCAAAGGATGCGTGACCACCGTTAATATCTCCActgtc         >  2:354832/1‑89 (MQ=255)

                                                                                 ggCTGACGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACGGCCa        <  2:4627915/93‑1 (MQ=255)

                                                                                       cGGCGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACGGCCATTAAcc  >  1:1832945/1‑93 (MQ=255)

ATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACCGTTAATAGCTCCACGGCCATTAACC  >  NC_000913/3676511‑3676690

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