BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,396,204	T→C	Q148R (CAG→CGG)	ydaA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,396,204	0	T	C	88.6%	649.7 / 65.2	254	Q148R (CAG→CGG)	ydaA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (16/13); new base C (112/113); total (128/126)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 254 total aligned reads displayed.

TTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGT  >  NC_000913/1396113‑1396293

tttGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTg                                                                                          <  1:6180357/93‑1 (MQ=255)

 ttGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCgg                                                                                         <  2:4657711/93‑1 (MQ=255)

 ttGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCgg                                                                                         >  2:4409662/1‑93 (MQ=255)

 ttGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCgg                                                                                         >  1:4882093/1‑93 (MQ=255)

  tGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGt                                                                                        <  1:5939307/93‑1 (MQ=255)

    aCCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCt                                                                                      >  2:5017710/1‑93 (MQ=255)

     ccAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCtt                                                                                     >  1:1288685/1‑93 (MQ=255)

     ccAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCtt                                                                                     >  2:7234809/1‑93 (MQ=255)

     ccAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCtt                                                                                     <  1:6519964/93‑1 (MQ=255)

     ccAGTTTTTCATTGAGCGCATCATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCtt                                                                                     >  1:1169051/1‑93 (MQ=255)

      cAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCttt                                                                                    >  1:1857485/1‑93 (MQ=255)

      cAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCAAGGCCGGTCttt                                                                                    >  2:4067661/1‑93 (MQ=255)

         tttttCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCAc                                                                                 >  2:3307258/1‑93 (MQ=255)

          ttgtCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCAcc                                                                                <  2:5869802/90‑1 (MQ=255)

            ttCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCAt                                                                              <  2:6976543/93‑1 (MQ=255)

             tCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATc                                                                             >  2:3876195/1‑93 (MQ=255)

               aTTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCca                                                                           >  2:3744123/1‑93 (MQ=255)

               aTTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCca                                                                           >  1:4653251/1‑93 (MQ=255)

               aTTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCca                                                                           <  2:6594341/93‑1 (MQ=255)

                 tGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCcaca                                                                         >  2:2454947/1‑93 (MQ=255)

                  taGCGCATTATGATACGGCTCTTCACTGGCGAGATTTACCCCCCCCAGCGCCTTGCCCCCTTCCCGCCCCCGCCCGGCTTTTACCATCCccac                                                                        >  1:5114560/2‑93 (MQ=255)

                    gcgcATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCg                                                                      >  1:4319117/1‑93 (MQ=255)

                     cgcATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCgg                                                                     <  1:5026464/93‑1 (MQ=255)

                      gcATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCggg                                                                    <  1:1621052/93‑1 (MQ=255)

                        aTTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCt                                                                  <  1:3479580/93‑1 (MQ=255)

                        aTTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCt                                                                  >  2:105007/1‑93 (MQ=255)

                         ttATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCtt                                                                 >  1:6026845/1‑93 (MQ=255)

                          tATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTg                                                                <  2:1566216/93‑1 (MQ=255)

                          tATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTg                                                                >  1:6933842/1‑93 (MQ=255)

                           aTGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTTACCATCCACACCGGGCTTgg                                                               >  1:2278579/1‑93 (MQ=255)

                            tGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTggg                                                              >  2:1464419/1‑93 (MQ=255)

                            tGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTggg                                                              >  2:2218211/1‑93 (MQ=255)

                            tGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTggg                                                              >  1:7332102/1‑93 (MQ=255)

                              aTACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCa                                                            >  1:4183136/1‑93 (MQ=255)

                              aTACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCa                                                            <  2:6196362/93‑1 (MQ=255)

                                 cGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCAttt                                                         <  2:188098/93‑1 (MQ=255)

                                   gCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGc                                                       <  2:4615068/93‑1 (MQ=255)

                                      ctTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTa                                                    >  2:3990892/1‑93 (MQ=255)

                                       ttCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTaa                                                   <  1:2023373/93‑1 (MQ=255)

                                       ttCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTaa                                                   <  1:1909089/93‑1 (MQ=255)

                                        tCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAAc                                                  >  1:3630706/1‑93 (MQ=255)

                                          aCTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAg                                                >  1:3759637/1‑93 (MQ=255)

                                           cTGGCGAGATTCCCCGCCCCCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGa                                               <  1:1415110/93‑1 (MQ=255)

                                           cTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGa                                               <  2:5456799/93‑1 (MQ=255)

                                           cTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGa                                               >  2:48026/1‑93 (MQ=255)

                                           cTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGa                                               >  2:3021440/1‑93 (MQ=255)

                                           cTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGa                                               <  2:1391096/93‑1 (MQ=255)

                                               cGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGcc                                           <  2:6425298/93‑1 (MQ=255)

                                               cGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGcc                                           <  2:3087980/93‑1 (MQ=255)

                                                gggaTTCCCCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCCCCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCa                                          <  1:534685/91‑1 (MQ=255)

                                                gagaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCa                                          >  2:5950405/1‑93 (MQ=255)

                                                gagaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCa                                          >  2:3504002/1‑93 (MQ=255)

                                                gagaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCa                                          <  1:5912209/93‑1 (MQ=255)

                                                 agaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAg                                         <  2:7122207/93‑1 (MQ=255)

                                                  gaTTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGt                                        >  1:4565327/1‑93 (MQ=255)

                                                    ttCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTcc                                      >  1:7015504/1‑93 (MQ=255)

                                                    ttCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTcc                                      <  2:5382381/93‑1 (MQ=255)

                                                    ttCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTcc                                      >  1:6474425/1‑93 (MQ=255)

                                                     tCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCg                                     >  1:3055940/1‑93 (MQ=255)

                                                       aCCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGtt                                   <  1:6066513/93‑1 (MQ=255)

                                                       aCCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGtt                                   >  2:5072768/1‑93 (MQ=255)

                                                       aCCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGtt                                   <  1:5138394/93‑1 (MQ=255)

                                                       aCCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGtt                                   >  2:4271195/1‑93 (MQ=255)

                                                       aCCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGtt                                   >  1:3239653/1‑93 (MQ=255)

                                                          gCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGc                                >  1:2209906/1‑93 (MQ=255)

                                                          gCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGc                                <  2:2496290/93‑1 (MQ=255)

                                                              ccaGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaa                            <  1:2943552/93‑1 (MQ=255)

                                                              ccaGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaa                            >  2:7115131/1‑93 (MQ=255)

                                                              ccaGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaa                            >  1:3817405/1‑93 (MQ=255)

                                                              ccaGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaa                            >  1:526439/1‑93 (MQ=255)

                                                               caGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaaa                           <  1:745526/93‑1 (MQ=255)

                                                               caGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaaa                           <  1:985080/93‑1 (MQ=255)

                                                                aGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaaaa                          >  1:685435/1‑93 (MQ=255)

                                                                aGCGCCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTaaaa                          >  1:5491378/1‑93 (MQ=255)

                                                                  cgcCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAAt                        <  1:5058914/93‑1 (MQ=255)

                                                                  cgcCTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAAt                        >  2:2767091/1‑93 (MQ=255)

                                                                    ccTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCa                      >  1:5718528/1‑93 (MQ=255)

                                                                    ccTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCa                      <  2:4561745/93‑1 (MQ=255)

                                                                    ccTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCa                      >  2:1801450/1‑93 (MQ=255)

                                                                     cTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCAc                     >  2:2100244/1‑93 (MQ=255)

                                                                     cTTGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCAc                     <  1:5558327/93‑1 (MQ=255)

                                                                      ttGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCAcc                    >  1:1713900/1‑93 (MQ=255)

                                                                       tGCCTCCTTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCg                   >  1:4417038/1‑93 (MQ=255)

                                                                        gCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGc                  <  1:6311504/93‑1 (MQ=255)

                                                                           tcctTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTc               >  2:1714056/1‑93 (MQ=255)

                                                                           tcctTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTc               >  2:5700070/1‑93 (MQ=255)

                                                                             ctTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCa             <  2:6340076/93‑1 (MQ=255)

                                                                             ctTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCa             <  1:1331933/93‑1 (MQ=255)

                                                                             ctTCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCa             <  1:6634072/93‑1 (MQ=255)

                                                                              ttCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAg            >  1:2628455/1‑93 (MQ=255)

                                                                               tCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGa           <  2:2020269/93‑1 (MQ=255)

                                                                               tCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGa           <  1:1558327/93‑1 (MQ=255)

                                                                               tCCGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGa           <  1:3403463/93‑1 (MQ=255)

                                                                                ccGGCCACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGAc          <  1:4125983/93‑1 (MQ=255)

                                                                                    ccACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTc      <  1:2730752/93‑1 (MQ=255)

                                                                                     cACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCa     >  2:6782715/1‑93 (MQ=255)

                                                                                     cACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCa     <  2:4506109/93‑1 (MQ=255)

                                                                                     cACGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCa     <  2:3207680/93‑1 (MQ=255)

                                                                                      aCGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCAt    >  1:7239062/1‑93 (MQ=255)

                                                                                      aCGGCCGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCAt    >  1:1745050/1‑93 (MQ=255)

TTTGACCAGTTTTTCATTGAGCGCATTATGATACGGCTCTTCACTGGCGAGATTCACCGCCACCAGCGCCTTGCCTCCTTCCGGCCACGGCTGGTCTTTCACCATCCACACCGGGCTTGGGCATTTGCGTAACAGATGCCAGTCCGTTGGCGTAAAAATCACCGCTTCCAGACGGTCATGT  >  NC_000913/1396113‑1396293

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