BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	233,314	T→C	K214K (AAA→AAG)	dniR ←	cytochrome c552

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	233,314	0	T	C	94.1%	1011.5 / 36.5	355	K214K (AAA→AAG)	dniR	cytochrome c552
Reads supporting (aligned to +/- strand): ref base T (14/7); new base C (177/157); total (191/164)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 372 total aligned reads displayed.

GCTGTTTTTGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGATTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAA  >  NC_000913/233224‑233406

 cTGTTTTTGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTc                                                                                           >  1:9702227/1‑93 (MQ=255)

   gTTTTTGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTccc                                                                                         <  1:4645232/93‑1 (MQ=255)

   gTTTTTGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTccc                                                                                         >  1:3564739/1‑93 (MQ=255)

    tttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGATTTCCCa                                                                                        >  1:608150/1‑93 (MQ=255)

    tttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCa                                                                                        <  1:5657180/93‑1 (MQ=255)

    tttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCa                                                                                        >  2:3177719/1‑93 (MQ=255)

    tttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCa                                                                                        <  2:4152658/93‑1 (MQ=255)

    tttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCa                                                                                        <  2:7558873/93‑1 (MQ=255)

    tttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCa                                                                                        >  2:5611433/1‑93 (MQ=255)

    tttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGGAACGACCAGAAGTCCGTGGATTTCCCa                                                                                        <  1:9790119/93‑1 (MQ=255)

     ttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCAc                                                                                       >  2:5997717/1‑93 (MQ=255)

     ttttGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCAc                                                                                       >  1:6998583/1‑93 (MQ=255)

         gagaATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACgcgc                                                                                   >  1:670254/1‑93 (MQ=255)

          agaATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCt                                                                                  >  2:8327777/1‑93 (MQ=255)

           gaATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCtt                                                                                 >  1:1819773/1‑93 (MQ=255)

           gaATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCtt                                                                                 >  1:3832756/1‑93 (MQ=255)

           gaATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGGTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCtt                                                                                 >  2:6882900/1‑93 (MQ=255)

              tatCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGt                                                                              >  2:485840/1‑93 (MQ=255)

               atCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGtt                                                                             <  1:5004130/93‑1 (MQ=255)

                tCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTc                                                                            >  1:2083697/1‑93 (MQ=255)

                 cACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGATTTCCCACGCGCTTTGTTCg                                                                           >  1:6629551/1‑93 (MQ=255)

                   cTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGtt                                                                         >  1:5719178/1‑93 (MQ=255)

                    tCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCATAAGTCCTTGGGCTTCTCACGTGTTTTGATCGttt                                                                        >  2:2075367/1‑93 (MQ=255)

                    tCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGttt                                                                        <  2:2580146/93‑1 (MQ=255)

                     cAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGtttt                                                                       >  1:2232308/1‑93 (MQ=255)

                      aaTGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTa                                                                      >  1:7977830/1‑93 (MQ=255)

                       aTGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTaa                                                                     <  2:1001469/93‑1 (MQ=255)

                       aTGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTaa                                                                     <  2:72335/93‑1 (MQ=255)

                        tGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAAt                                                                    >  1:875611/1‑93 (MQ=255)

                         gCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAAtt                                                                   <  1:8079780/93‑1 (MQ=255)

                         gCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAAtt                                                                   <  1:3387561/93‑1 (MQ=255)

                          ccAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTg                                                                  >  2:9148674/1‑93 (MQ=255)

                           cAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGc                                                                 >  2:1695357/1‑93 (MQ=255)

                            cgCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTAGTTTTAATTGcc                                                                <  1:6393850/92‑1 (MQ=255)

                            aGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGcc                                                                <  1:3689868/93‑1 (MQ=255)

                            aGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGcc                                                                <  2:6559027/93‑1 (MQ=255)

                             gCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCt                                                               <  1:9553721/93‑1 (MQ=255)

                              cATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCtt                                                              <  2:5184696/93‑1 (MQ=255)

                                ttttAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCa                                                            <  1:1044456/93‑1 (MQ=255)

                                ttttAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCa                                                            >  2:2961492/1‑93 (MQ=255)

                                ttttAGGCACGTACTGATTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCa                                                            <  1:8974399/93‑1 (MQ=255)

                                 tttAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCAt                                                           <  2:1193323/93‑1 (MQ=255)

                                 tttAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCAt                                                           <  2:4306826/93‑1 (MQ=255)

                                  ttAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATg                                                          <  2:1903392/93‑1 (MQ=255)

                                  ttAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATg                                                          <  2:285542/93‑1 (MQ=255)

                                    aGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGAc                                                        <  1:8376614/93‑1 (MQ=255)

                                     ggCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGAcc                                                       <  2:584896/93‑1 (MQ=255)

                                     ggCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGAcc                                                       >  1:7146023/1‑93 (MQ=255)

                                       cACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCg                                                     <  1:1451853/93‑1 (MQ=255)

                                        aCGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCgc                                                    >  1:2505416/1‑92 (MQ=255)

                                         cGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGAc                                                   >  1:3135883/1‑93 (MQ=255)

                                          gTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGAcc                                                  >  2:7672659/1‑93 (MQ=255)

                                          gTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGAcc                                                  <  1:4569537/93‑1 (MQ=255)

                                          gTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGAcc                                                  <  1:2006207/93‑1 (MQ=255)

                                           tACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCt                                                 >  2:5401358/1‑93 (MQ=255)

                                           tACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCt                                                 <  2:9400930/93‑1 (MQ=255)

                                           tACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCt                                                 <  1:7310227/93‑1 (MQ=255)

                                           tACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCt                                                 >  2:8234114/1‑93 (MQ=255)

                                            aCTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCtt                                                <  1:5264808/93‑1 (MQ=255)

                                            aCTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCtt                                                >  2:8653166/1‑93 (MQ=255)

                                              tGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTcg                                              >  1:4514461/1‑93 (MQ=255)

                                               gCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTcgc                                             >  2:335037/1‑93 (MQ=255)

                                                 ttCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTcgccg                                           >  1:2438453/1‑93 (MQ=255)

                                                      ttCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGtt                                      >  1:6003081/1‑93 (MQ=255)

                                                       tCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTta                                     >  1:2939013/1‑93 (MQ=255)

                                                        ccTGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTtat                                    >  2:2353197/1‑93 (MQ=255)

                                                          tGCGGCAACGGTAACGACCAGAAGTCCGTGGATTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATaa                                  <  2:666928/93‑1 (MQ=255)

                                                          tGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATaa                                  >  2:968960/1‑93 (MQ=255)

                                                          tGCGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATaa                                  >  1:783547/1‑93 (MQ=255)

                                                            cGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGc                                >  2:1341367/1‑93 (MQ=255)

                                                            cGGCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGc                                <  1:4182087/93‑1 (MQ=255)

                                                              gCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCg                              >  1:6762839/1‑93 (MQ=255)

                                                              gCAACGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCg                              >  2:7977660/1‑93 (MQ=255)

                                                                aaCGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCt                            >  1:2524988/1‑93 (MQ=255)

                                                                 aCGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTa                           >  1:9335652/1‑93 (MQ=255)

                                                                 aCGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTa                           >  2:4934838/1‑93 (MQ=255)

                                                                  cGGTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTAc                          <  2:4050856/93‑1 (MQ=255)

                                                                  cGGAAACGACCAGAAAGCCGTGGACTTCCCACGCGCTTTGTTGGTTTTAATTGCCTTAATGGCCCGACCTTCGCCGCTGTTATAAGCCGCTAc                          <  2:4628267/93‑1 (MQ=255)

                                                                   ggTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACg                         >  2:6342003/1‑93 (MQ=255)

                                                                    gTAACGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACgg                        >  2:8831330/1‑93 (MQ=255)

                                                                      aaCGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGGCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTc                      <  1:3991469/93‑1 (MQ=255)

                                                                      aaCGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCCGTTATAAGCCGCTACGGTc                      >  2:6107015/1‑93 (MQ=255)

                                                                       aCGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCa                     >  2:7850247/1‑93 (MQ=255)

                                                                        cGACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAg                    <  2:3861297/93‑1 (MQ=255)

                                                                         gACCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGa                   <  1:8627392/93‑1 (MQ=255)

                                                                          aCCAGAAGTCCGTGGATTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGaa                  >  2:9874437/1‑93 (MQ=255)

                                                                          aCCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGaa                  >  2:8501055/1‑93 (MQ=255)

                                                                          aCCAGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGaa                  >  2:9204766/1‑93 (MQ=255)

                                                                             aGAAGTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGcc               >  2:2499547/1‑93 (MQ=255)

                                                                              gAAGTCCGTGGATTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCa              <  1:5438841/93‑1 (MQ=255)

                                                                                 gTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTc           <  1:3237785/93‑1 (MQ=255)

                                                                                 gTCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTc           >  2:814238/1‑93 (MQ=255)

                                                                                  tCCGTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCCACGGGCAGGAGCCAGTcc          >  2:145318/1‑92 (MQ=255)

                                                                                     gTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCg       >  2:7458201/1‑93 (MQ=255)

                                                                                     gTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCg       <  1:4925684/93‑1 (MQ=255)

                                                                                     gTGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTCCTATAAACCGCTACGGTCCGAAGCCAGTCGCCg       >  1:7538973/1‑93 (MQ=255)

                                                                                      tGGATTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGt      >  2:4488394/1‑93 (MQ=255)

                                                                                      tGGATTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGt      >  1:5835631/1‑93 (MQ=255)

                                                                                      tGGACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGt      >  2:3610519/1‑93 (MQ=255)

                                                                                       ggACTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGGCGCCGTc     >  2:3950337/1‑93 (MQ=255)

GCTGTTTTTGAGAATATCACTCAATGCCAGCATTTTAGGCACGTACTGCTTCGTTTCCTGCGGCAACGGTAACGACCAGAAGTCCGTGGATTTCCCACGCGCTTTGTTCGTTTTAATTGCCTTCATGACCCGACCTTCGCCGCTGTTATAAGCCGCTACGGTCAGAAGCCAGTCGCCGTCAAA  >  NC_000913/233224‑233406

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