BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	628,816	A→G	V98V (GTA→GTG)	ybdB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	628,816	0	A	G	99.7%	1028.2 / NA	318	V98V (GTA→GTG)	ybdB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (200/117); total (200/118)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

GGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628726‑628905

 gCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTgcgc                                                                                        >  1:2910768/1‑93 (MQ=255)

    ggATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggc                                                                                     >  1:3434132/1‑93 (MQ=255)

    ggATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggc                                                                                     >  1:1869549/1‑93 (MQ=255)

     gATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggcg                                                                                    <  1:7677638/93‑1 (MQ=255)

     gATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggcg                                                                                    >  1:1987296/1‑93 (MQ=255)

     gATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggcg                                                                                    >  2:8514179/1‑93 (MQ=255)

     gATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggcg                                                                                    >  1:9940826/1‑93 (MQ=255)

       tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  <  2:1254647/93‑1 (MQ=255)

          atgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGc                                                                               >  1:2088950/1‑93 (MQ=255)

          atgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGc                                                                               >  2:8288530/1‑93 (MQ=255)

           tgatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGcc                                                                              >  1:4916882/1‑93 (MQ=255)

            gatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCa                                                                             >  2:5976729/1‑93 (MQ=255)

              tgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGc                                                                           <  2:959975/93‑1 (MQ=255)

              tgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGc                                                                           >  1:830404/1‑93 (MQ=255)

                    gcgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGc                                                                     >  2:5518397/1‑93 (MQ=255)

                      gACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCAt                                                                   <  1:2290685/93‑1 (MQ=255)

                      gACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCAt                                                                   >  1:2656628/1‑93 (MQ=255)

                       aCGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCTCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATc                                                                  >  2:3772132/1‑93 (MQ=255)

                        cgaaCAGTTTGTGGTAGGCACAGAACTTAATTCACCCCACCATCGCCCGGTGACTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 <  1:6904376/90‑1 (MQ=255)

                        cGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 >  1:609010/1‑93 (MQ=255)

                        cGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 >  1:8143656/1‑93 (MQ=255)

                        cGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 <  2:911338/93‑1 (MQ=255)

                         ggACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCtt                                                                >  1:309222/1‑93 (MQ=255)

                          gACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTg                                                               <  1:6659609/93‑1 (MQ=255)

                           aCAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              >  2:4641272/1‑93 (MQ=255)

                           aCAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              <  2:9324408/93‑1 (MQ=255)

                           aCAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTgg                                                              >  2:3166820/1‑93 (MQ=255)

                            cAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGt                                                             <  1:2367839/93‑1 (MQ=255)

                            cAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGt                                                             >  2:7320393/1‑93 (MQ=255)

                             aGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTc                                                            >  2:4368295/1‑93 (MQ=255)

                               tgtgtgGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCgg                                                          >  2:5736287/1‑93 (MQ=255)

                                  gtgGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCaa                                                       >  2:389244/1‑93 (MQ=255)

                                   tgGTAGGCACAGAACTTAATGCAACAAACCATCGCCCGGTGATTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCaaa                                                      >  1:1651027/1‑93 (MQ=255)

                                     gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAAt                                                    >  1:6557785/1‑93 (MQ=255)

                                      tAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATc                                                   >  1:5327732/1‑93 (MQ=255)

                                      tAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATc                                                   >  2:1358243/1‑93 (MQ=255)

                                       aGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCa                                                  <  2:507639/93‑1 (MQ=255)

                                       aGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCa                                                  <  1:4693978/93‑1 (MQ=255)

                                        ggCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTGGGCAAAATCag                                                 >  2:8855561/1‑93 (MQ=255)

                                        ggCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCag                                                 >  1:7477324/1‑93 (MQ=255)

                                           acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGGAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              >  1:6255715/1‑93 (MQ=255)

                                           acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              >  1:5740726/1‑93 (MQ=255)

                                           acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              <  2:5919578/93‑1 (MQ=255)

                                           acaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGc                                              >  2:1561578/1‑93 (MQ=255)

                                            caGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCt                                             <  1:5211195/93‑1 (MQ=255)

                                             aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            <  1:4788976/93‑1 (MQ=255)

                                             aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            >  1:9641224/1‑93 (MQ=255)

                                              gAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTgg                                           >  2:1860900/1‑93 (MQ=255)

                                                aCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGa                                         >  2:4724921/1‑93 (MQ=255)

                                                aCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGa                                         >  2:5009041/1‑93 (MQ=255)

                                                aCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGa                                         <  1:5924270/93‑1 (MQ=255)

                                                aCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGa                                         >  1:4306828/1‑93 (MQ=255)

                                                aCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGa                                         <  2:6276029/93‑1 (MQ=255)

                                                  ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       >  2:8315942/1‑93 (MQ=255)

                                                  ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       >  1:4147733/1‑93 (MQ=255)

                                                   tAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAt                                      <  2:9798109/93‑1 (MQ=255)

                                                   tAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAt                                      <  1:3389248/93‑1 (MQ=255)

                                                    aaTGCAAAACACCATCGCCAGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtc                                     <  2:8940099/93‑1 (MQ=255)

                                                      tGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgt                                   >  1:3991465/1‑93 (MQ=255)

                                                       gCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcgtc                                  >  2:2926215/1‑93 (MQ=255)

                                                         aaCACACCGTCGCGCGGTGTCTGAGGGAAAGTTCCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATcgtcgt                                <  2:2443658/93‑1 (MQ=255)

                                                          acacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtt                               >  2:7126806/1‑93 (MQ=255)

                                                          acacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtt                               <  1:5616806/93‑1 (MQ=255)

                                                          acacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtt                               >  2:3422315/1‑93 (MQ=255)

                                                           ctcacaATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              <  2:3146516/91‑1 (MQ=255)

                                                           cacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              >  2:1173051/1‑93 (MQ=255)

                                                           cacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              <  2:8091511/93‑1 (MQ=255)

                                                           cacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              >  1:7169258/1‑93 (MQ=255)

                                                           cacacCATCGCCCGGTGTCGGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              <  2:5302615/93‑1 (MQ=255)

                                                            acacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtttt                             >  2:2125608/1‑93 (MQ=255)

                                                            acacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtttt                             <  2:856054/93‑1 (MQ=255)

                                                             ccccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTc                            <  1:1003473/91‑1 (MQ=255)

                                                             cacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTc                            >  2:456367/1‑93 (MQ=255)

                                                              acCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           >  2:4811567/1‑93 (MQ=255)

                                                              acCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                           <  2:6671227/93‑1 (MQ=255)

                                                               ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  2:7764699/93‑1 (MQ=255)

                                                               ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  2:418527/93‑1 (MQ=255)

                                                               ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  2:6958872/93‑1 (MQ=255)

                                                               ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  1:2509835/93‑1 (MQ=255)

                                                                cATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                         <  2:7496520/93‑1 (MQ=255)

                                                                 aTCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATg                        >  2:1045099/1‑93 (MQ=255)

                                                                  tCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGa                       >  1:6785907/1‑93 (MQ=255)

                                                                  acGCCCGGTGGTGGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGa                       <  1:1088603/92‑1 (MQ=255)

                                                                     cccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACa                    >  1:3301824/1‑93 (MQ=255)

                                                                     cccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACa                    <  1:9186011/93‑1 (MQ=255)

                                                                      ctggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAg                   <  1:793949/91‑1 (MQ=255)

                                                                       cGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgg                  <  1:1208711/93‑1 (MQ=255)

                                                                       cGGTGTCGGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTTGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgg                  <  2:7347492/93‑1 (MQ=255)

                                                                        ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 >  1:1583745/1‑93 (MQ=255)

                                                                         gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                <  2:2157534/93‑1 (MQ=255)

                                                                          tgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggc               <  2:3952606/93‑1 (MQ=255)

                                                                             cTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcggc            <  2:7655504/93‑1 (MQ=255)

                                                                              tGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg           >  1:8641088/1‑93 (MQ=255)

                                                                              tGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg           >  1:9976300/1‑93 (MQ=255)

                                                                               gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          >  2:8622135/1‑93 (MQ=255)

                                                                               gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          >  1:4078818/1‑93 (MQ=255)

                                                                               gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          >  1:2780296/1‑93 (MQ=255)

                                                                                  ggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgc       <  1:7244992/93‑1 (MQ=255)

                                                                                     aaGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctgc    <  1:9693045/93‑1 (MQ=255)

                                                                                       ggTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCat  >  1:8541248/1‑92 (MQ=255)

GGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628726‑628905

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