BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	922,398	T→C	V88A (GTG→GCG)	yljA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	922,398	0	T	C	99.5%	641.1 / ‑4.5	199	V88A (GTG→GCG)	yljA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (89/109); total (90/109)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

Only 100 of 208 total aligned reads displayed.

AACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTGCCTA  >  NC_000913/922309‑922487

aaCGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGaa                                                                                        >  2:1029960/1‑93 (MQ=255)

   gCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAaaaa                                                                                     >  1:4208479/1‑93 (MQ=255)

   gCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACaa                                                                                     <  2:6536315/93‑1 (MQ=255)

    cAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGTGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAg                                                                                    <  1:5795860/93‑1 (MQ=255)

      aTTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTa                                                                                  <  2:1586606/93‑1 (MQ=255)

      aTTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTa                                                                                  >  1:6585088/1‑93 (MQ=255)

      aTTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTa                                                                                  <  1:201153/93‑1 (MQ=255)

       ttGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGCACAAGGAc                                                                                 >  1:4258482/1‑93 (MQ=255)

       ttGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTAc                                                                                 >  2:381333/1‑93 (MQ=255)

       ttGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTAc                                                                                 <  2:3976625/93‑1 (MQ=255)

          aTGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTAcgcg                                                                              <  1:2269420/93‑1 (MQ=255)

           tGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGa                                                                             <  1:1068986/93‑1 (MQ=255)

            gCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAg                                                                            <  2:140413/93‑1 (MQ=255)

            gCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAg                                                                            <  2:2401400/93‑1 (MQ=255)

            gCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAg                                                                            <  2:3200919/93‑1 (MQ=255)

             cTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAgg                                                                           <  1:3810390/93‑1 (MQ=255)

              tCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCGAAGTGGCGATGGCGAACAAGTACGCGAggg                                                                          >  1:2907148/1‑93 (MQ=255)

              tCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAggg                                                                          <  2:4570268/93‑1 (MQ=255)

              tCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAggg                                                                          >  2:3715216/1‑93 (MQ=255)

              tCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAggg                                                                          >  1:3590580/1‑93 (MQ=255)

               tgCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGga                                                                         <  2:6424974/92‑1 (MQ=255)

               cGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCTATGGCGAACAAGTACGCGAGGga                                                                         <  1:4610861/93‑1 (MQ=255)

               cGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGga                                                                         <  2:5844462/93‑1 (MQ=255)

                gCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGgag                                                                        >  1:141941/1‑93 (MQ=255)

                  tGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGaa                                                                      >  1:6618820/1‑93 (MQ=255)

                  tGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGaa                                                                      >  1:1852362/1‑93 (MQ=255)

                    ttCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATg                                                                    <  1:6307266/93‑1 (MQ=255)

                    ttCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATg                                                                    <  1:3965098/93‑1 (MQ=255)

                     tCACTACCAGGGGAAGGCCATTTGCTGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGa                                                                   <  2:3355695/93‑1 (MQ=255)

                     tCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGa                                                                   <  2:2787637/93‑1 (MQ=255)

                        cTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCa                                                                <  1:875979/93‑1 (MQ=255)

                        cTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCa                                                                >  2:5603629/1‑93 (MQ=255)

                          aCCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATc                                                              <  1:3263971/93‑1 (MQ=255)

                             aGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCAt                                                           <  2:199338/93‑1 (MQ=255)

                             aGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCAt                                                           <  1:4569722/93‑1 (MQ=255)

                              ggggAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCAtt                                                          <  2:6277279/93‑1 (MQ=255)

                              ggggAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCAtt                                                          <  1:90393/93‑1 (MQ=255)

                               gggAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTg                                                         >  2:3907030/1‑93 (MQ=255)

                               gggAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTg                                                         <  2:297419/93‑1 (MQ=255)

                                ggAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGc                                                        >  2:1869417/1‑93 (MQ=255)

                                ggAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGc                                                        <  2:1455537/93‑1 (MQ=255)

                                ggAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGc                                                        <  1:6007802/93‑1 (MQ=255)

                                 gAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCt                                                       <  2:5489630/93‑1 (MQ=255)

                                   aGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCtgt                                                     >  1:6076658/1‑93 (MQ=255)

                                   aGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCtgt                                                     >  1:3114819/1‑93 (MQ=255)

                                       cATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTAc                                                 >  2:822246/1‑93 (MQ=255)

                                        aTTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACg                                                <  2:4273043/93‑1 (MQ=255)

                                        aTTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACg                                                >  2:48698/1‑93 (MQ=255)

                                         tttGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGc                                               >  1:2722157/1‑93 (MQ=255)

                                         tttGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGc                                               >  2:859241/1‑93 (MQ=255)

                                          ttGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCt                                              <  2:4873370/93‑1 (MQ=255)

                                           tGCGTAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTa                                             <  2:1639241/93‑1 (MQ=255)

                                             cGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGa                                           <  2:3130348/93‑1 (MQ=255)

                                              ggAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGaa                                          <  2:4138165/93‑1 (MQ=255)

                                               gAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGaaa                                         >  1:640695/1‑93 (MQ=255)

                                               gAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCCAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGaaa                                         >  2:3932909/1‑93 (MQ=255)

                                                 gTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGaaaaa                                       >  1:5029384/1‑93 (MQ=255)

                                                 gTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGaaaaa                                       <  2:2260034/93‑1 (MQ=255)

                                                  tCTTTACCGCCGAGGTTGCAGAAACCGAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATGCATTGCTGAGTACGCTAGAAAAAg                                      >  1:3143269/1‑93 (MQ=255)

                                                    tttACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGcc                                    <  1:2521740/93‑1 (MQ=255)

                                                     ttACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCt                                   >  2:3889222/1‑93 (MQ=255)

                                                      tACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTg                                  >  2:6683816/1‑93 (MQ=255)

                                                       aCCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGa                                 >  2:4966492/1‑93 (MQ=255)

                                                        ccgccgAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGaa                                <  1:3559203/93‑1 (MQ=255)

                                                        ccgccgAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGaa                                >  1:5345132/1‑93 (MQ=255)

                                                         cgccgAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAAt                               >  2:2560309/1‑93 (MQ=255)

                                                         cgccgAGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAAt                               <  1:171205/93‑1 (MQ=255)

                                                              aGGTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAgg                          >  1:4879065/1‑93 (MQ=255)

                                                               ggTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGc                         >  1:1162464/1‑93 (MQ=255)

                                                               ggTTGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGc                         >  2:2434357/1‑93 (MQ=255)

                                                                 ttGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCAt                       <  1:1250891/93‑1 (MQ=255)

                                                                 ttGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCAt                       <  1:3599703/93‑1 (MQ=255)

                                                                  tGCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATa                      >  2:260717/1‑93 (MQ=255)

                                                                   gCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATaa                     <  2:1147327/93‑1 (MQ=255)

                                                                   gCAGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATaa                     <  1:757493/93‑1 (MQ=255)

                                                                     aGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATaaaa                   >  1:3340694/1‑93 (MQ=255)

                                                                     aGAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATaaaa                   >  2:4753527/1‑93 (MQ=255)

                                                                      gAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATaaaat                  >  2:5202856/1‑92 (MQ=255)

                                                                      gAAACCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATaaaaa                  <  2:638760/93‑1 (MQ=255)

                                                                       aaaCCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAAt                 >  1:309889/1‑93 (MQ=255)

                                                                        aaCCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAAtt                <  1:5682524/93‑1 (MQ=255)

                                                                        aaCCAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAAtt                <  1:544093/93‑1 (MQ=255)

                                                                          ccAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTgg              <  2:4822485/93‑1 (MQ=255)

                                                                          ccAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTgg              <  2:2130600/93‑1 (MQ=255)

                                                                          ccAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTgg              <  1:516232/93‑1 (MQ=255)

                                                                          ccAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTgg              >  1:4463372/1‑93 (MQ=255)

                                                                           cAAAGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTggg             <  2:3628607/93‑1 (MQ=255)

                                                                             aaGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTggggg           >  1:2955091/1‑93 (MQ=255)

                                                                             aaGTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTggggg           >  1:162730/1‑93 (MQ=255)

                                                                               gTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAg         >  2:4501941/1‑93 (MQ=255)

                                                                               gTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAg         <  2:721619/93‑1 (MQ=255)

                                                                               gTGGCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAg         >  2:734587/1‑93 (MQ=255)

                                                                                tggcgatggcgaACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAgg        <  2:6054095/93‑1 (MQ=255)

                                                                                tggcgatggcgaACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAgg        <  2:4387702/93‑1 (MQ=255)

                                                                                  gCGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTg      <  1:1564323/93‑1 (MQ=255)

                                                                                   cGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTGc     <  1:2328011/93‑1 (MQ=255)

                                                                                   cGATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTGc     >  2:1248942/1‑93 (MQ=255)

                                                                                    gATGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTGcc    <  1:5290532/93‑1 (MQ=255)

                                                                                      tGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTGCCTa  <  2:5404668/93‑1 (MQ=255)

                                                                                      tGGCGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTGCCTa  <  1:6538350/93‑1 (MQ=255)

AACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAGGGAGAATGAGCATCCATTGCTGTGTACGCTAGAAAAAGCCTGAATGCAGGCATAAAAATTGGGGGAGGTGCCTA  >  NC_000913/922309‑922487

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