BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,108,438	C→A	noncoding (56/76 nt)	pheV →	tRNA‑Phe

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,108,438	0	C	A	99.4%	652.1 / NA	173	noncoding (56/76 nt)	pheV	tRNA‑Phe
Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (72/100); minor base T (1/0); total (73/100)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAA  >  NC_000913/3108348‑3108526

cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGa                                                                                        <  1:392013/93‑1 (MQ=255)

cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGa                                                                                        >  2:320657/1‑93 (MQ=255)

  agacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAtt                                                                                      >  2:3754495/1‑93 (MQ=255)

   gacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTc                                                                                     >  1:545820/1‑93 (MQ=255)

    acgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTcc                                                                                    >  1:685108/1‑93 (MQ=255)

    acgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTcc                                                                                    <  1:2059082/93‑1 (MQ=255)

      gagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGa                                                                                  <  2:2378942/93‑1 (MQ=255)

      gagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGa                                                                                  <  1:2974915/93‑1 (MQ=255)

       agGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAGAATCCCCGTGTCCTTGGTTAGATTCagag                                                                                 >  1:6350642/1‑93 (MQ=255)

          cGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTcc                                                                              <  1:3072874/93‑1 (MQ=255)

           gAATCATGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                             <  1:6086066/93‑1 (MQ=255)

           gAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCg                                                                             >  1:3642670/1‑93 (MQ=37)

            aaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                            >  1:5738721/1‑93 (MQ=25)

            aaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                            >  1:2783791/1‑93 (MQ=25)

            aaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCgg                                                                            >  2:6249119/1‑93 (MQ=25)

             aTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCggg                                                                           <  2:355673/93‑1 (MQ=21)

              tCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGc                                                                          >  2:4499311/1‑93 (MQ=18)

                aGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcac                                                                        >  2:1117591/1‑93 (MQ=14)

                aGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcac                                                                        <  1:1646068/93‑1 (MQ=14)

                   tttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcaccac                                                                     <  1:833514/93‑1 (MQ=14)

                    ttAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACt                                                                    >  2:1709771/1‑93 (MQ=17)

                       aTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAt                                                                 <  1:5060189/93‑1 (MQ=25)

                       aTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAt                                                                 <  2:1314778/93‑1 (MQ=25)

                        tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAAtt                                                                <  2:6183209/93‑1 (MQ=37)

                         gcgcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTc                                                               >  1:315164/1‑93 (MQ=255)

                         gcgcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTc                                                               <  1:5924070/93‑1 (MQ=255)

                          cgcCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCt                                                              <  2:3420881/93‑1 (MQ=255)

                            ccccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                                            >  2:3883702/1‑93 (MQ=255)

                             cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAACTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                           <  2:104714/93‑1 (MQ=255)

                             cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                           <  1:1132531/93‑1 (MQ=255)

                              ccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAg                                                          >  1:4075164/1‑93 (MQ=255)

                              ccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAg                                                          <  2:5995556/93‑1 (MQ=255)

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                                         >  2:1364179/1‑93 (MQ=255)

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                                         <  2:522610/93‑1 (MQ=255)

                               cGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                                         >  1:6458362/1‑93 (MQ=255)

                                gTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                                        >  1:731007/1‑93 (MQ=255)

                                 ttGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAc                                                       >  2:2189342/1‑93 (MQ=255)

                                  tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                                      <  1:2399172/93‑1 (MQ=255)

                                    cccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCg                                                    <  2:1373590/93‑1 (MQ=255)

                                    cccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCg                                                    <  1:1970554/93‑1 (MQ=255)

                                     ccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGc                                                   >  1:4559225/1‑93 (MQ=255)

                                     ccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGc                                                   >  1:4837750/1‑93 (MQ=255)

                                     ccGGATAGCTCAGTCGGGAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCGGc                                                   >  1:852630/1‑93 (MQ=255)

                                       ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                                                 >  1:6239886/1‑93 (MQ=255)

                                       ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                                                 <  2:4785631/93‑1 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                >  2:427078/1‑93 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                >  1:3409496/1‑93 (MQ=255)

                                        gATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                >  2:5970554/1‑93 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                             <  2:720738/93‑1 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                             >  1:5853855/1‑93 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                             <  1:293876/93‑1 (MQ=255)

                                           aGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                                             <  1:2682555/93‑1 (MQ=255)

                                               cAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCg                                         <  1:3420714/93‑1 (MQ=255)

                                               cAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCg                                         >  2:5711274/1‑93 (MQ=255)

                                                aGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCgg                                        <  2:482124/93‑1 (MQ=255)

                                                aGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCgg                                        >  2:3456275/1‑93 (MQ=255)

                                                 gTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCggg                                       <  1:6266747/93‑1 (MQ=255)

                                                 gTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCggg                                       >  2:1902639/1‑93 (MQ=255)

                                                  tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                      <  2:462287/93‑1 (MQ=255)

                                                    ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttt                                    <  2:598643/93‑1 (MQ=255)

                                                     gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCATAAGGCGGGtttt                                   <  1:5152556/93‑1 (MQ=255)

                                                      tAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttt                                  <  1:4571590/93‑1 (MQ=255)

                                                        gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc                                <  2:208955/93‑1 (MQ=255)

                                                        gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc                                <  2:2548235/93‑1 (MQ=255)

                                                         agCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCt                               <  2:5325497/93‑1 (MQ=255)

                                                          gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                              <  2:2906096/93‑1 (MQ=255)

                                                          gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                              >  1:2493979/1‑93 (MQ=255)

                                                          gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                              >  2:6060170/1‑93 (MQ=255)

                                                          gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt                              <  2:430997/93‑1 (MQ=255)

                                                          gCAGGGGATTGAAAAGCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGATCCCGCCCACAAGGCGGGTTTTTGCtt                              <  2:2694632/93‑1 (MQ=255)

                                                            aGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtttt                            <  2:2636444/93‑1 (MQ=255)

                                                             ggggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                           <  1:95406/93‑1 (MQ=255)

                                                             ggggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                           <  1:369792/93‑1 (MQ=255)

                                                              gggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg                          <  1:3201935/93‑1 (MQ=255)

                                                              gggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg                          >  1:799446/1‑93 (MQ=255)

                                                               ggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGa                         <  1:6242444/93‑1 (MQ=255)

                                                                gATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGAt                        >  1:3964296/1‑93 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                       <  1:280525/93‑1 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                       <  2:3589612/93‑1 (MQ=255)

                                                                 aTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                       <  1:19883/93‑1 (MQ=255)

                                                                  ttGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                      >  2:5101066/1‑93 (MQ=255)

                                                                  ttGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                      >  2:1457883/1‑93 (MQ=255)

                                                                     aaaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAc                   <  2:2531234/93‑1 (MQ=255)

                                                                      aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACa                  <  2:4074798/93‑1 (MQ=255)

                                                                        aTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAAt                <  2:1738745/93‑1 (MQ=255)

                                                                          ccccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATaa              <  1:2140977/93‑1 (MQ=255)

                                                                           cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAc             <  2:3228208/93‑1 (MQ=255)

                                                                           cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAc             >  2:5557829/1‑93 (MQ=255)

                                                                            ccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAcc            >  1:2712724/1‑93 (MQ=255)

                                                                            ccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAcc            >  1:1061137/1‑93 (MQ=255)

                                                                             cGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCt           <  1:6173185/93‑1 (MQ=255)

                                                                              gtgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCtt          >  2:883016/1‑93 (MQ=255)

                                                                               tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc         <  1:4254756/93‑1 (MQ=255)

                                                                               tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc         <  1:2415713/93‑1 (MQ=255)

                                                                               tgtCCTGGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc         <  2:797553/93‑1 (MQ=255)

                                                                                gtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCa        <  2:5354200/93‑1 (MQ=255)

                                                                                gtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCa        <  1:3546987/93‑1 (MQ=255)

                                                                                  ccTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACg      <  1:3295394/93‑1 (MQ=255)

                                                                                  ccTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCAGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACg      >  2:3711835/1‑93 (MQ=255)

                                                                                   cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGa     <  1:5251687/93‑1 (MQ=255)

CGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAA  >  NC_000913/3108348‑3108526

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