BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,201,729	A→G	V15A (GTG→GCG)	bacA ←	BacA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,201,729	0	A	G	99.3%	895.3 / ‑3.6	287	V15A (GTG→GCG)	bacA	BacA
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (127/158); total (129/158)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 293 total aligned reads displayed.

GCCGTGTCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCACACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGAACTACT  >  NC_000913/3201639‑3201818

  cGTGTCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGTCAGAAATTCTGTCAATCCTTCGACCGCAcc                                                                                       >  2:6455130/1‑93 (MQ=255)

  cGTGTCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCAcc                                                                                       >  1:4867586/1‑93 (MQ=255)

   gtgtCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCAccc                                                                                      >  2:5979777/1‑93 (MQ=255)

    tgtCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCa                                                                                     >  1:3687555/1‑93 (MQ=255)

     gtCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCaa                                                                                    <  1:1421138/93‑1 (MQ=255)

     gtCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCaa                                                                                    >  1:7335922/1‑93 (MQ=255)

     gtCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCaa                                                                                    >  2:6084122/1‑93 (MQ=255)

       cGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAata                                                                                  >  2:2353581/1‑93 (MQ=255)

        gCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAatat                                                                                 <  1:7379474/93‑1 (MQ=255)

        gCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAatat                                                                                 >  1:5598693/1‑93 (MQ=255)

        gCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAatat                                                                                 <  1:5306392/93‑1 (MQ=255)

         cccTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATAtt                                                                                >  2:230478/1‑93 (MQ=255)

         cccTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATAtt                                                                                <  1:7796261/93‑1 (MQ=255)

          ccTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTg                                                                               <  1:6593887/93‑1 (MQ=255)

           cTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGc                                                                              >  1:2730324/1‑93 (MQ=255)

            tCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGCAATTCTGTCAATCCTTCGACCGCACCCAATATTGcc                                                                             >  1:880888/1‑93 (MQ=255)

            tCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGcc                                                                             >  2:3883336/1‑93 (MQ=255)

            tCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGcc                                                                             <  1:7864478/93‑1 (MQ=255)

              aaaCCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGc                                                                           <  2:7685550/93‑1 (MQ=255)

                 ccccAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTAt                                                                        >  2:3400861/1‑93 (MQ=255)

                  cccAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATc                                                                       >  1:6017838/1‑93 (MQ=255)

                   ccAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATca                                                                      <  1:6255697/93‑1 (MQ=255)

                   ccAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATca                                                                      <  1:397391/93‑1 (MQ=255)

                   ccAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATca                                                                      <  2:4096652/93‑1 (MQ=255)

                    caacaaGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATcag                                                                     >  2:1952481/1‑93 (MQ=255)

                     aacaaGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATcagc                                                                    >  2:6906388/1‑93 (MQ=255)

                      acaaGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATcagca                                                                   <  2:4662307/93‑1 (MQ=255)

                       caaGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATcagcag                                                                  >  1:3622709/1‑93 (MQ=255)

                        aaGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCagcagc                                                                 >  1:3367611/1‑93 (MQ=255)

                        aaGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCagcagc                                                                 >  2:7377923/1‑93 (MQ=255)

                        aaGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCagcagc                                                                 <  2:4851405/93‑1 (MQ=255)

                          gTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGa                                                               >  1:955526/1‑93 (MQ=255)

                            gACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGt                                                             <  2:7654025/93‑1 (MQ=255)

                            gACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGt                                                             <  1:4227923/93‑1 (MQ=255)

                             aCCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTg                                                            >  2:1344780/1‑93 (MQ=255)

                             aCCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTg                                                            <  1:2164387/93‑1 (MQ=255)

                               cGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCa                                                          >  1:1688408/1‑93 (MQ=255)

                               cGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCa                                                          <  2:118296/93‑1 (MQ=255)

                                gACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCat                                                         <  1:4100245/93‑1 (MQ=255)

                                gACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCat                                                         <  1:3848906/93‑1 (MQ=255)

                                 aCAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCata                                                        >  2:6379643/1‑93 (MQ=255)

                                 aCAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCata                                                        >  1:4740336/1‑93 (MQ=255)

                                 aCAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCACACCCAATATTGCCGCTATCAGCAGCGAGTGCata                                                        >  1:3897080/1‑93 (MQ=255)

                                  cAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCatat                                                       <  1:3748008/93‑1 (MQ=255)

                                  cAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCatat                                                       <  1:6985007/93‑1 (MQ=255)

                                    ataatCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCg                                                     >  1:4542576/1‑93 (MQ=255)

                                      aatCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCt                                                   <  2:4390132/93‑1 (MQ=255)

                                       atCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCtc                                                  <  1:1963530/93‑1 (MQ=255)

                                        tCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCtca                                                 >  2:7082057/1‑93 (MQ=255)

                                        tCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCtca                                                 <  2:5859958/93‑1 (MQ=255)

                                        tCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCtca                                                 >  1:1198737/1‑93 (MQ=255)

                                        tCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCtca                                                 <  1:3652151/93‑1 (MQ=255)

                                          atatGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCtcatc                                               >  2:94763/1‑93 (MQ=255)

                                            atGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCaa                                             <  1:4390027/93‑1 (MQ=255)

                                             tGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAAt                                            >  1:7178869/1‑93 (MQ=255)

                                             tGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGAATATCGCTCATCAAt                                            >  1:2342504/1‑93 (MQ=255)

                                                cccGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATaaa                                         <  2:4007183/93‑1 (MQ=255)

                                                cccGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATaaa                                         <  2:2795692/93‑1 (MQ=255)

                                                 ccGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAAc                                        >  1:2423591/1‑93 (MQ=255)

                                                   gTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAAccc                                      >  2:1739207/1‑93 (MQ=255)

                                                   gTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAAccc                                      >  1:5566295/1‑93 (MQ=255)

                                                    tGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAAcccc                                     >  1:2837636/1‑93 (MQ=255)

                                                      cTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTa                                   <  2:4731809/93‑1 (MQ=255)

                                                       tGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTaa                                  >  1:2252753/1‑93 (MQ=255)

                                                       tGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTaa                                  >  1:4429230/1‑93 (MQ=255)

                                                       tGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTaa                                  >  1:1019922/1‑93 (MQ=255)

                                                       tGGATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTaa                                  >  2:4536205/1‑93 (MQ=255)

                                                         gATGCCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGATATCAGCAGCGAGTGCATATCGCTCATCAATAACCCCCTAAat                                <  2:2915847/93‑1 (MQ=255)

                                                         gATACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAat                                >  2:5127211/1‑93 (MQ=255)

                                                           tACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAatta                              >  2:7175691/1‑93 (MQ=255)

                                                           tACCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAatta                              <  2:6826923/93‑1 (MQ=255)

                                                            aCCGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAattat                             <  1:5704225/93‑1 (MQ=255)

                                                             ccGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAattatt                            <  1:3110323/93‑1 (MQ=255)

                                                             ccGGCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAattatt                            <  2:6249226/93‑1 (MQ=255)

                                                               ggCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTaa                          <  2:2456786/93‑1 (MQ=255)

                                                               ggCAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTaa                          <  1:4991439/93‑1 (MQ=255)

                                                                 cAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTaaaa                        <  2:1517413/93‑1 (MQ=255)

                                                                 cAGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTaaaa                        <  2:549491/93‑1 (MQ=255)

                                                                  aGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCCATAAACCCCTAAATTATTAAAAt                       >  2:4302140/1‑93 (MQ=255)

                                                                  aGAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAAt                       <  1:3327794/93‑1 (MQ=255)

                                                                   gAAATTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATg                      >  2:3535035/1‑93 (MQ=255)

                                                                    aaaTTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGt                     <  1:1881454/93‑1 (MQ=255)

                                                                    aaaTTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGt                     >  1:5115072/1‑93 (MQ=255)

                                                                     aaTTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTa                    <  2:7514667/93‑1 (MQ=255)

                                                                     aaTTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTa                    >  1:2673135/1‑93 (MQ=255)

                                                                     aaTTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTa                    <  1:6350825/93‑1 (MQ=255)

                                                                      aTTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTAc                   >  1:7632242/1‑93 (MQ=255)

                                                                      aTTCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTAc                   <  1:354407/93‑1 (MQ=255)

                                                                       ttCTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTAcc                  >  1:5755845/1‑93 (MQ=255)

                                                                         cTGTCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGc                <  2:3310481/93‑1 (MQ=255)

                                                                            tCAATCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTg             >  1:3015389/1‑93 (MQ=255)

                                                                               aTCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTcc          <  2:2084084/93‑1 (MQ=255)

                                                                                tCCTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCg         >  2:268042/1‑93 (MQ=255)

                                                                                accTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCg         <  1:2474589/92‑1 (MQ=255)

                                                                                  cTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGaa       >  2:4997198/1‑93 (MQ=255)

                                                                                  cTTCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGaa       >  2:311080/1‑93 (MQ=255)

                                                                                  cTTAGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGaa       <  2:5344643/93‑1 (MQ=255)

                                                                                   ttCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGAac      >  2:6542475/1‑93 (MQ=255)

                                                                                    tCGACCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGAact     >  2:1794985/1‑93 (MQ=255)

                                                                                       aCCGCACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGAactact  <  2:3688897/93‑1 (MQ=255)

GCCGTGTCGCCCTCAAACCCCAACAAGTGACCGACAATAATCATATGGCCCGTGCTGGATACCGGCAGAAATTCTGTCAATCCTTCGACCACACCCAATATTGCCGCTATCAGCAGCGAGTGCATATCGCTCATCAATAAACCCCTAAATTATTAAAATGTACCGCTTGTCCGAACTACT  >  NC_000913/3201639‑3201818

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