BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	747,952	A→G	V85A (GTA→GCA)	ybgO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	747,952	0	A	G	99.6%	842.1 / ‑6.3	270	V85A (GTA→GCA)	ybgO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (136/132); total (137/133)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 277 total aligned reads displayed.

CGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGATC  >  NC_000913/747862‑748041

  aCTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCTCCCCCGGATAGGGGTTTGCCCa                                                                                       <  2:6440768/93‑1 (MQ=255)

  aCTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCa                                                                                       <  2:6937900/93‑1 (MQ=255)

  aCTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCa                                                                                       <  2:7566447/93‑1 (MQ=255)

    tttttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAgg                                                                                     <  1:5407999/93‑1 (MQ=255)

    tttttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGGTTGCCCAgg                                                                                     >  2:435414/1‑93 (MQ=255)

     ttttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGc                                                                                    <  1:1612213/93‑1 (MQ=255)

        tCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTa                                                                                 <  1:778415/93‑1 (MQ=255)

          ccAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGa                                                                               <  2:1010366/93‑1 (MQ=255)

            aGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACa                                                                             <  1:1346715/93‑1 (MQ=255)

             gACTGGCAACGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACAt                                                                            <  2:4462420/93‑1 (MQ=255)

              gcTGGCATCGTAATCGACACCGTTGTATATCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATg                                                                           <  1:150671/92‑1 (MQ=255)

              aCTGGCATCGTAATCGACACCGTTGTACGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATg                                                                           <  1:532791/93‑1 (MQ=255)

               cTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGc                                                                          >  2:6046296/1‑93 (MQ=255)

                tGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCt                                                                         <  2:2667499/93‑1 (MQ=255)

                tGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCt                                                                         >  2:3636360/1‑93 (MQ=255)

                tGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCt                                                                         >  2:416071/1‑93 (MQ=255)

                 ggCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTc                                                                        <  2:1336234/93‑1 (MQ=255)

                   cATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGc                                                                      <  1:7350720/93‑1 (MQ=255)

                    aTCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGAGTTTGCCCAGGCGTAGACATGCTCGCt                                                                     >  1:3223833/1‑93 (MQ=255)

                     tCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCtt                                                                    >  1:3996525/1‑93 (MQ=255)

                     tCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCtt                                                                    >  2:777906/1‑93 (MQ=255)

                     tCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCtt                                                                    >  2:3528553/1‑93 (MQ=255)

                      cGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCttt                                                                   <  2:785300/93‑1 (MQ=255)

                       gTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCttta                                                                  >  2:5498395/1‑92 (MQ=255)

                        tAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCa                                                                 <  1:7274009/93‑1 (MQ=255)

                        tAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCa                                                                 >  2:1597709/1‑93 (MQ=255)

                           tCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAg                                                              >  1:4085064/1‑93 (MQ=255)

                           tCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAg                                                              >  1:2657759/1‑93 (MQ=255)

                            cGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGt                                                             <  2:3160263/93‑1 (MQ=255)

                             gACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGtt                                                            <  2:2852477/93‑1 (MQ=255)

                              acacCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTg                                                           <  1:6849244/93‑1 (MQ=255)

                              acacCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTg                                                           >  2:2950606/1‑93 (MQ=255)

                               cacCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGc                                                          <  1:7906432/93‑1 (MQ=255)

                               cacCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGc                                                          <  2:744674/93‑1 (MQ=255)

                                cccGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGcc                                                         <  1:6251268/92‑1 (MQ=255)

                                acCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGcc                                                         <  1:6121197/93‑1 (MQ=255)

                                  cGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGt                                                       <  1:1014065/93‑1 (MQ=255)

                                   gTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGtt                                                      >  2:2520533/1‑93 (MQ=255)

                                   gTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGtt                                                      <  1:7570059/93‑1 (MQ=255)

                                    ttGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTg                                                     >  1:4202590/1‑93 (MQ=255)

                                    ttGTATGTAACTCCGCGTTGATAATAACAATCCTGCACCCCCGGACAGGGGCTTGCCCTGTCGTAGACATGCTCGCTATCAAAGTTGGCGTTg                                                     >  2:825265/1‑93 (MQ=255)

                                     tGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTgg                                                    <  1:4863247/93‑1 (MQ=255)

                                     tGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTgg                                                    <  2:4091184/93‑1 (MQ=255)

                                       tATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTa                                                  >  2:2234047/1‑93 (MQ=255)

                                          gTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAttg                                               >  2:916534/1‑93 (MQ=255)

                                           tCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAttgt                                              <  2:7747188/93‑1 (MQ=255)

                                              cGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATtgttgt                                           <  1:7753357/93‑1 (MQ=255)

                                               gCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTc                                          <  1:3969701/93‑1 (MQ=255)

                                                  gAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTcaca                                       >  1:1056009/1‑93 (MQ=255)

                                                  gAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTcaca                                       >  1:4267551/1‑93 (MQ=255)

                                                  gAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTcaca                                       <  1:4343545/93‑1 (MQ=255)

                                                  gAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGATGGTATTGTTGTcaca                                       >  2:3121925/1‑93 (MQ=255)

                                                   aGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACaa                                      >  1:3668519/1‑93 (MQ=255)

                                                   aGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACaa                                      >  1:497769/1‑93 (MQ=255)

                                                    gTTGATAATAACGATCCTGGGCCCCCGGATAGGTGTTTGCCCAGGCGTAGGCATGCTCGCTTTCAAAGTTGCGGTGGGTATTGTTGTCACAAt                                     <  2:5755397/93‑1 (MQ=255)

                                                    gTTGATAATAACGATCATGCACCCCCGAATAGGGAATTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAAt                                     <  1:7023494/93‑1 (MQ=255)

                                                     ttGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATa                                    <  2:6975356/93‑1 (MQ=255)

                                                     ttGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATa                                    <  1:283476/93‑1 (MQ=255)

                                                       gATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGa                                  <  1:3336610/93‑1 (MQ=255)

                                                       gATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGa                                  >  2:7102613/1‑93 (MQ=255)

                                                        ataataACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGGGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGAc                                 >  1:2221056/1‑93 (MQ=255)

                                                         taataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACg                                >  2:6423929/1‑93 (MQ=255)

                                                         taataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACg                                <  1:2987437/93‑1 (MQ=255)

                                                         taataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACg                                <  1:7219222/93‑1 (MQ=255)

                                                          aataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACgg                               <  2:2317952/93‑1 (MQ=255)

                                                           ataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACggg                              >  1:5490553/1‑93 (MQ=255)

                                                           ataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACggg                              <  2:5160744/93‑1 (MQ=255)

                                                            taaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGa                             >  1:1968837/1‑93 (MQ=255)

                                                             aaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGaa                            >  1:469508/1‑93 (MQ=255)

                                                             aaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGaa                            <  2:3366808/93‑1 (MQ=255)

                                                             aaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGaa                            <  2:4437118/93‑1 (MQ=255)

                                                               cGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAAtt                          >  1:2357169/1‑93 (MQ=255)

                                                               cGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAAtt                          >  2:2029914/1‑93 (MQ=255)

                                                               cGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCAGTTGGTATTGTTGTCACAATAGACGGGAAtt                          >  1:3483146/1‑93 (MQ=255)

                                                                   ccTGCGCCCCCGGATAGGGGTTTGCCGCGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTaa                      <  1:6749649/93‑1 (MQ=255)

                                                                   ccTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTaa                      <  2:1823847/93‑1 (MQ=255)

                                                                    cTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAat                     >  2:7251776/1‑93 (MQ=255)

                                                                    cTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAat                     <  2:6059609/93‑1 (MQ=255)

                                                                    cTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAat                     <  1:3783994/93‑1 (MQ=255)

                                                                     tGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAata                    <  2:597989/93‑1 (MQ=255)

                                                                      gCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAatat                   <  1:3153637/93‑1 (MQ=255)

                                                                      gCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAatat                   <  1:847678/93‑1 (MQ=255)

                                                                       cACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATc                  <  1:6510065/93‑1 (MQ=255)

                                                                       cACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATc                  <  1:2449795/93‑1 (MQ=255)

                                                                        aCCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCg                 >  1:1709012/1‑93 (MQ=255)

                                                                           cccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCAAAATAGACGGGAATTTTAATATCGTCa              >  1:2749236/1‑93 (MQ=255)

                                                                            ccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAg             >  1:3467018/1‑93 (MQ=255)

                                                                             cGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGa            <  2:5440955/93‑1 (MQ=255)

                                                                               gATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGANGNGAATTTTAATNTCGTCAGAtt          >  1:7260698/1‑93 (MQ=255)

                                                                               gATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGAtt          >  2:566656/1‑93 (MQ=255)

                                                                               gATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGAtt          <  2:7180349/93‑1 (MQ=255)

                                                                                 tAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTcc        >  2:6576729/1‑93 (MQ=255)

                                                                                 tAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTcc        <  2:4858074/93‑1 (MQ=255)

                                                                                 tAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTcc        >  2:3926060/1‑93 (MQ=255)

                                                                                   ggggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTcccg      >  2:2442799/1‑92 (MQ=255)

                                                                                   ggggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCa      >  2:4528726/1‑93 (MQ=255)

                                                                                     ggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGa    >  1:4585058/1‑93 (MQ=255)

                                                                                     ggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGa    <  1:5151916/93‑1 (MQ=255)

                                                                                     ggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGGATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATGCCCAGa    >  2:7047910/1‑93 (MQ=255)

                                                                                       tttGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGATc  <  1:5896486/93‑1 (MQ=255)

CGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGATC  >  NC_000913/747862‑748041

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: `‑`