BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,655,220	T→C	D579G (GAC→GGC)	hscA ←	HscA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,655,220	0	T	C	99.2%	365.3 / ‑4.5	117	D579G (GAC→GGC)	hscA	HscA
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (63/53); total (64/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 126 total aligned reads displayed.

TCAGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGTCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATCA  >  NC_000913/2655131‑2655309

tCAGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCaa                                                                                        >  2:2996597/1‑93 (MQ=255)

  aGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAAca                                                                                      <  1:484369/93‑1 (MQ=255)

    cgcACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAAcatc                                                                                    >  1:4504832/1‑93 (MQ=255)

     gcACGACGAACCGACTTGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAAcatca                                                                                   >  2:2677791/1‑93 (MQ=255)

      cACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAAcatcat                                                                                  <  1:1092564/93‑1 (MQ=255)

         gacgaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGc                                                                               >  1:297647/1‑93 (MQ=255)

             aaacGACTGGTCCATGCGGCGAGCGGCGAAAACCTGGGGTTTTTTTTTTACGTTTTTTATTGCTTTTTTGATGGGGGCAAAATCATCGCCCtt                                                                           >  2:1566064/4‑92 (MQ=255)

             aaCCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTg                                                                           >  2:3277677/1‑93 (MQ=255)

             aaCCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTg                                                                           >  1:862270/1‑93 (MQ=255)

              aCCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTgc                                                                          <  1:483855/93‑1 (MQ=255)

              aCCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTgc                                                                          >  2:4400430/1‑93 (MQ=255)

                cGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTgcgc                                                                        >  1:2443173/1‑93 (MQ=255)

                cGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTgcgc                                                                        >  1:4464767/1‑93 (MQ=255)

                cGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTgcgc                                                                        >  2:3889268/1‑93 (MQ=255)

                 gACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGcc                                                                       >  1:98922/1‑93 (MQ=255)

                   cTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCAc                                                                     <  2:817639/93‑1 (MQ=255)

                    tGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACt                                                                    >  1:3592133/1‑93 (MQ=255)

                     ggTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACtt                                                                   <  1:3564958/93‑1 (MQ=255)

                        ccATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCAc                                                                <  1:4288022/93‑1 (MQ=255)

                         cATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCCACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACt                                                               >  1:986976/1‑93 (MQ=255)

                          aTGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTc                                                              >  1:3620183/1‑93 (MQ=255)

                            gcggcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAg                                                            <  1:4610024/93‑1 (MQ=255)

                            gcggcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAg                                                            >  1:887372/1‑93 (MQ=255)

                             cggcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGTCAACATCATCGCCCTGCGCCACTTCACTCAgg                                                           >  1:1456669/1‑93 (MQ=255)

                             cggcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAgg                                                           <  1:767877/93‑1 (MQ=255)

                              ggcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGt                                                          <  2:555505/93‑1 (MQ=255)

                              ggcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGt                                                          >  2:4373015/1‑93 (MQ=255)

                              ggcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGt                                                          >  2:1712678/1‑93 (MQ=255)

                               gcgagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTg                                                         >  2:789719/1‑93 (MQ=255)

                                  agcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCATCATCATCGCCCTGCGCCACTTCACTCAGGTGAGc                                                      <  1:2654268/93‑1 (MQ=255)

                                   gcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCg                                                     <  1:3986634/93‑1 (MQ=255)

                                    cggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCgg                                                    <  1:3521636/93‑1 (MQ=255)

                                    cggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCgg                                                    >  2:27757/1‑93 (MQ=255)

                                      gcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGtca                                                  >  1:4427892/1‑93 (MQ=255)

                                      gcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGgca                                                  <  1:2849258/93‑1 (MQ=255)

                                        gaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGgcagc                                                <  1:2294536/93‑1 (MQ=255)

                                          aaTCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCAt                                              <  1:2932212/93‑1 (MQ=255)

                                          aaTCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCAt                                              <  2:4373854/93‑1 (MQ=255)

                                           aTCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATc                                             >  1:2474034/1‑93 (MQ=255)

                                           aTCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATc                                             <  2:4535214/93‑1 (MQ=255)

                                           aTCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATc                                             >  2:3416173/1‑93 (MQ=255)

                                            tCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCg                                            >  1:3901077/1‑93 (MQ=255)

                                            tCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCg                                            >  2:3857874/1‑93 (MQ=255)

                                            tCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCg                                            <  1:1582768/93‑1 (MQ=255)

                                              cTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTc                                          >  2:1725371/1‑93 (MQ=255)

                                               cggggTTGTTTGTCTACGCTTTTAATCGCTTGTTCGATGGCGACAACAGCATCGCCCTGCGCCACTTGACTCACGCGAGCGGCAGCATCGGCa                                         >  2:2399819/2‑93 (MQ=255)

                                                gggTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCaa                                        >  1:4569877/1‑93 (MQ=255)

                                                 ggTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAAt                                       <  2:615799/93‑1 (MQ=255)

                                                 ggTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAAt                                       >  1:1127986/1‑93 (MQ=255)

                                                  gtttgtttGTCTACGTTTTTCATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATg                                      <  2:371752/93‑1 (MQ=255)

                                                  gtttgtttGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATg                                      >  2:1146076/1‑93 (MQ=255)

                                                  gtttgtttGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATg                                      >  1:2239069/1‑93 (MQ=255)

                                                   tttgtttgTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGa                                     <  2:1837533/93‑1 (MQ=255)

                                                    ttgtttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGAc                                    >  1:3514968/1‑93 (MQ=255)

                                                     tgtttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGAcc                                   <  1:3391386/93‑1 (MQ=255)

                                                     tgtttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGAcc                                   >  2:4560552/1‑93 (MQ=255)

                                                      ttttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCt                                  <  1:1494760/92‑1 (MQ=255)

                                                       tttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTg                                 <  1:4295775/93‑1 (MQ=255)

                                                       tttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTg                                 >  2:1100294/1‑93 (MQ=255)

                                                       tttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTg                                 >  1:3569871/1‑93 (MQ=255)

                                                        ttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCTGCAGCATCGTCAATGACCTGa                                >  2:3287611/1‑93 (MQ=255)

                                                        ttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGa                                >  1:1292837/1‑93 (MQ=255)

                                                        ttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGa                                <  1:2588921/93‑1 (MQ=255)

                                                        ttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGa                                <  1:2588920/93‑1 (MQ=255)

                                                        ttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGa                                <  2:194628/93‑1 (MQ=255)

                                                         tgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGAc                               <  2:552222/93‑1 (MQ=255)

                                                         tgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGAc                               <  1:2907691/93‑1 (MQ=255)

                                                         tgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGAc                               <  1:517066/93‑1 (MQ=255)

                                                          gtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACg                              >  2:717638/1‑93 (MQ=255)

                                                           tCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGt                             >  2:207088/1‑93 (MQ=255)

                                                             tACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTc                           <  2:4536619/93‑1 (MQ=255)

                                                             tACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTc                           <  2:2725760/93‑1 (MQ=255)

                                                             tACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTc                           >  1:329059/1‑93 (MQ=255)

                                                              tcGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCt                          <  2:3654664/92‑1 (MQ=255)

                                                              aCGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCt                          >  1:4688686/1‑93 (MQ=255)

                                                              aCGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCt                          >  2:2058878/1‑93 (MQ=255)

                                                                 tttttAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTgcg                       <  2:570513/93‑1 (MQ=255)

                                                                  ttttAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTgcgg                      <  1:3280412/93‑1 (MQ=255)

                                                                   tttAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTgcggc                     <  1:2603633/93‑1 (MQ=255)

                                                                    ttAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTgcggcc                    >  2:4090973/1‑92 (MQ=255)

                                                                      aaTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCt                  >  2:2792708/1‑93 (MQ=255)

                                                                       aTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCtt                 >  1:697252/1‑93 (MQ=255)

                                                                       aTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCtt                 >  2:2814356/1‑93 (MQ=255)

                                                                       aTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCtt                 >  2:2740042/1‑93 (MQ=255)

                                                                       aTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCAGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCtt                 >  1:3294189/1‑93 (MQ=255)

                                                                         cGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTaa               >  2:3948247/1‑93 (MQ=255)

                                                                          gCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAAc              <  2:4077084/93‑1 (MQ=255)

                                                                           cTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACa             >  1:4158465/1‑93 (MQ=255)

                                                                            ttGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAg            >  1:3488322/1‑93 (MQ=255)

                                                                               ttgtaTGGCGCCAGCATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAgcgc         <  1:1232444/89‑1 (MQ=255)

                                                                               ttCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAgcgc         <  1:170351/93‑1 (MQ=255)

                                                                               ttCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAgcgc         <  1:54009/93‑1 (MQ=255)

                                                                               ttCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAgcgc         <  1:620038/93‑1 (MQ=255)

                                                                                tCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGcgcg        <  2:1561001/93‑1 (MQ=255)

                                                                                  gATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGc      <  2:2520306/93‑1 (MQ=255)

                                                                                  gATGGCGCAAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGc      <  2:2520305/93‑1 (MQ=255)

                                                                                    tGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCAt    >  1:4320694/1‑93 (MQ=255)

                                                                                      gcgcCACCATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATca  <  2:2529607/93‑1 (MQ=255)

                                                                                      gcgcCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATca  >  1:4144486/1‑93 (MQ=255)

                                                                                      gcgcCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGAGCGGCATca  <  1:4185993/93‑1 (MQ=255)

TCAGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGTCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATCA  >  NC_000913/2655131‑2655309

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