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.5%	634.4 / ‑4.7	198	D579G (GAC→GGC)	hscA	HscA
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (100/97); total (100/98)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 204 total aligned reads displayed.

TTCAGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGTCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATCA  >  NC_000913/2655130‑2655309

ttCAGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCa                                                                                         >  2:1883785/1‑93 (MQ=255)

   aGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAAca                                                                                      >  2:6123885/1‑93 (MQ=255)

    gcgcACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAAcat                                                                                     <  2:1370817/93‑1 (MQ=255)

      gcACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAAcatca                                                                                   >  1:2155366/1‑93 (MQ=255)

        acgacgAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACatcatc                                                                                 >  2:4138800/1‑93 (MQ=255)

         cgacgaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCg                                                                                >  2:2185105/1‑93 (MQ=255)

         cgacgaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCg                                                                                <  2:1486079/93‑1 (MQ=255)

           acgaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGcc                                                                              >  1:2369927/1‑93 (MQ=255)

            cgaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGccc                                                                             <  2:4625942/93‑1 (MQ=255)

            cgaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGccc                                                                             <  1:4357487/93‑1 (MQ=255)

             gaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCt                                                                            >  1:3054543/1‑93 (MQ=255)

             gaACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCt                                                                            >  1:1050906/1‑93 (MQ=255)

             gaACCGACTGGTCCATGCGGCGAGCGGAGAAATCCTGGGTTTGTTTTTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCt                                                                            >  1:2889690/1‑93 (MQ=255)

              aaCCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCtc                                                                           >  2:5668904/1‑92 (MQ=255)

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

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

               aCCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGGTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCCACATCATCGTCCTgc                                                                          >  1:3574769/1‑93 (MQ=255)

                ccGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTgcc                                                                         >  2:5874281/1‑92 (MQ=255)

                 cGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTgcgc                                                                        <  2:2725424/93‑1 (MQ=255)

                  gACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGcc                                                                       <  2:4115902/93‑1 (MQ=255)

                  gACTGGTCCATGCGACGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGcc                                                                       <  1:4974377/93‑1 (MQ=255)

                   aCTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCa                                                                      <  2:3301938/93‑1 (MQ=255)

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

                     tGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACt                                                                    <  1:1034261/93‑1 (MQ=255)

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

                        tCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCa                                                                 <  2:4369060/93‑1 (MQ=255)

                        tCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCa                                                                 <  1:154999/93‑1 (MQ=255)

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

                          cATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACt                                                               <  1:2962024/93‑1 (MQ=255)

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

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

                           aTGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTc                                                              <  1:19757/93‑1 (MQ=255)

                           aTGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTc                                                              <  1:980186/93‑1 (MQ=255)

                            tGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCa                                                             >  1:1442982/1‑93 (MQ=255)

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

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

                                  gagcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAg                                                       <  1:3324445/93‑1 (MQ=255)

                                    gcggcgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCg                                                     >  1:2289951/1‑93 (MQ=255)

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

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

                                        cgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCACCATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGgcag                                                 <  2:4817949/93‑1 (MQ=255)

                                        cgaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGgcag                                                 >  2:1041618/1‑93 (MQ=255)

                                         gaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGgcagc                                                >  2:3862718/1‑93 (MQ=255)

                                         gaAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGgcagc                                                >  1:954670/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                              ccTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGt                                           <  1:4832620/93‑1 (MQ=255)

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

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

                                                 gggTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCaa                                        <  2:1407228/93‑1 (MQ=255)

                                                 gggTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCaa                                        <  2:3066065/93‑1 (MQ=255)

                                                 gggTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCCTCaa                                        <  1:4715251/93‑1 (MQ=255)

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

                                                    tttgtttgTCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGa                                     >  2:5263851/1‑93 (MQ=255)

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

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

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

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

                                                       gtttgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCt                                  <  1:4766381/93‑1 (MQ=255)

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

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

                                                          tgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGAc                               >  1:5057672/1‑93 (MQ=255)

                                                          tgtCTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGAc                               >  1:2497361/1‑93 (MQ=255)

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

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

                                                             cTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGtt                            <  1:4485841/93‑1 (MQ=255)

                                                             cTACGTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGtt                            <  2:6195270/93‑1 (MQ=255)

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

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

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

                                                                 gTTTTTAATCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTgc                        >  1:6093132/1‑93 (MQ=255)

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

                                                                       aaTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCt                  <  2:3628283/93‑1 (MQ=255)

                                                                        aTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCtt                 <  1:5954923/93‑1 (MQ=255)

                                                                        aTCGCTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCtt                 <  1:5385312/93‑1 (MQ=255)

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

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

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

                                                                            cTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACa             <  2:3682663/93‑1 (MQ=255)

                                                                            cTTGTTCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACa             <  1:2546472/93‑1 (MQ=255)

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

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

                                                                               tttCGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAgcg          <  1:6443228/92‑1 (MQ=255)

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

                                                                                  cGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCgg       >  2:2315830/1‑93 (MQ=255)

                                                                                  cGATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCgg       >  2:1181112/1‑93 (MQ=255)

                                                                                   gATGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGc      >  1:3954302/1‑93 (MQ=255)

                                                                                    aTGGCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCa     >  1:5600968/1‑93 (MQ=255)

                                                                                      ggCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATc   >  2:5164861/1‑93 (MQ=255)

                                                                                      ggCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATc   <  2:4288053/93‑1 (MQ=255)

                                                                                      ggCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATc   >  2:2209066/1‑93 (MQ=255)

                                                                                      ggCGCCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATc   >  1:1374893/1‑93 (MQ=255)

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

                                                                                       gcgcCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATca  <  2:3163417/93‑1 (MQ=255)

                                                                                       gcgcCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATca  <  2:6704371/93‑1 (MQ=255)

TTCAGCGCACGACGAACCGACTGGTCCATGCGGCGAGCGGCGAAATCCTGGGTTTGTTTGTCTACGTTTTTAATCGCTTGTTCGATGGCGTCAACATCATCGCCCTGCGCCACTTCACTCAGGTGAGCGGCAGCATCGTCAATGACCTGACGTTCTGCGGCGCTTAACAGCGCGGCATCA  >  NC_000913/2655130‑2655309

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