BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,249,311	G→A	S216S (TCG→TCA)	yqjG →	putative transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,249,311	0	G	A	84.4%	706.3 / 96.7	262	S216S (TCG→TCA)	yqjG	putative transferase
Reads supporting (aligned to +/- strand): ref base G (22/19); new base A (108/113); total (130/132)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 262 total aligned reads displayed.

TTTATGACACCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGACTACGCT  >  NC_000913/3249220‑3249400

 ttATGACACCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACt                                                                                         >  2:2466524/1‑93 (MQ=255)

   aTGACACCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTgg                                                                                       >  2:1726135/1‑93 (MQ=255)

    tGACACCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGc                                                                                      >  1:4402013/1‑93 (MQ=255)

     gACACCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCa                                                                                     >  1:2209602/1‑93 (MQ=255)

      acacCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCAc                                                                                    >  1:4028029/1‑93 (MQ=255)

      acacCGATAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCAc                                                                                    <  2:2111985/93‑1 (MQ=255)

       cacCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACg                                                                                   >  1:1743653/1‑93 (MQ=255)

       cacCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACg                                                                                   <  1:4317503/93‑1 (MQ=255)

        acCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGa                                                                                  >  2:2252176/1‑93 (MQ=255)

        acCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGa                                                                                  >  2:5350692/1‑93 (MQ=255)

         ccGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGAc                                                                                 <  2:4956427/93‑1 (MQ=255)

          cGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACt                                                                                >  1:4514513/1‑93 (MQ=255)

           gTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTg                                                                               >  1:4214468/1‑93 (MQ=255)

            ttAACAACGGCGTGTATAAAGCTGGTTTTGCCCCCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTgg                                                                              <  2:2009943/93‑1 (MQ=255)

            ttAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTgg                                                                              >  1:4592514/1‑93 (MQ=255)

             tAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGa                                                                             <  1:1472184/93‑1 (MQ=255)

              aacaacGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAg                                                                            >  1:372898/1‑93 (MQ=255)

              aacaacGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAg                                                                            <  1:2815868/93‑1 (MQ=255)

                caacGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCa                                                                          >  1:4103382/1‑93 (MQ=255)

                caacGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCa                                                                          >  2:6200571/1‑93 (MQ=255)

                 aacGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAg                                                                         <  2:6445067/93‑1 (MQ=255)

                  acGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGa                                                                        >  1:5293273/1‑93 (MQ=255)

                  acGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGa                                                                        >  1:5829619/1‑93 (MQ=255)

                  acGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGa                                                                        >  1:1026358/1‑93 (MQ=255)

                   cGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGAt                                                                       >  2:5463938/1‑93 (MQ=255)

                    ggCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGAtt                                                                      >  2:3678270/1‑93 (MQ=255)

                     gCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGAttt                                                                     >  2:2862355/1‑93 (MQ=255)

                       gtgtATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTa                                                                   >  2:2767714/1‑93 (MQ=255)

                        tgtATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAg                                                                  >  2:3769458/1‑93 (MQ=255)

                           ataAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGcc                                                               <  1:5245917/93‑1 (MQ=255)

                           ataAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGcc                                                               <  1:3067853/93‑1 (MQ=255)

                            taAAGCTGTTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCa                                                              <  1:855591/93‑1 (MQ=255)

                              aaGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGc                                                            <  2:5987844/93‑1 (MQ=255)

                               aGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCa                                                           >  1:3588659/1‑93 (MQ=255)

                               aGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCa                                                           >  1:2181923/1‑93 (MQ=255)

                               aGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCa                                                           <  2:839254/93‑1 (MQ=255)

                               aGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCa                                                           >  2:529360/1‑93 (MQ=255)

                               aGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCa                                                           <  1:6488216/93‑1 (MQ=255)

                                gCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCAt                                                          >  1:4662655/1‑93 (MQ=255)

                                gCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCAt                                                          >  1:6345517/1‑93 (MQ=255)

                                 cTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                                         <  2:1271813/93‑1 (MQ=255)

                                 cTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATc                                                         <  2:911948/93‑1 (MQ=255)

                                  tGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCg                                                        <  1:1820149/93‑1 (MQ=255)

                                   ggTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGt                                                       >  2:5617092/1‑93 (MQ=255)

                                    gTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGtt                                                      >  2:1242383/1‑93 (MQ=255)

                                     ttttGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTa                                                     >  1:1266318/1‑93 (MQ=255)

                                        tGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCt                                                  >  1:3015537/1‑93 (MQ=255)

                                        tGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCt                                                  >  1:5390182/1‑93 (MQ=255)

                                         gCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTg                                                 <  1:3369410/93‑1 (MQ=255)

                                           caccaGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGAc                                               >  2:4109016/1‑93 (MQ=255)

                                           caccaGCCAGGAAGCTTACGACGAGGCGGTGGCCAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGAc                                               <  1:6194207/93‑1 (MQ=255)

                                             ccagccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCg                                             <  1:3767977/93‑1 (MQ=255)

                                                gccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCa                                          <  1:2965866/93‑1 (MQ=255)

                                                gccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCa                                          >  1:2050387/1‑93 (MQ=255)

                                                gccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCa                                          <  1:2684402/93‑1 (MQ=255)

                                                gccagGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCa                                          >  1:3395125/1‑93 (MQ=255)

                                                  cagGAAGCTTCCGACGAGGCGGTGGCGCAAGTGTTTGAATCACTGGCACGACTGGAGCAGAGTTTAGGCCAGCATCGTTACCTGACCGGCAAc                                        >  1:3622703/1‑93 (MQ=255)

                                                   agGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAAcc                                       >  2:2993179/1‑93 (MQ=255)

                                                   agGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAAcc                                       <  2:2355295/93‑1 (MQ=255)

                                                   agGAAGCTTACGACGAGGCGGTGGCGAAAGTGGTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAAcc                                       >  1:5508058/1‑93 (MQ=255)

                                                    ggAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCa                                      >  1:5167499/1‑93 (MQ=255)

                                                    ggAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCa                                      >  1:3938636/1‑93 (MQ=255)

                                                    ggAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGAAACCa                                      >  1:3740999/1‑93 (MQ=255)

                                                      aaGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGc                                    >  1:4184664/1‑93 (MQ=255)

                                                       aGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCt                                   >  1:6570359/1‑93 (MQ=255)

                                                        gCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTa                                  <  2:4323299/93‑1 (MQ=255)

                                                        gCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTa                                  <  1:4482486/93‑1 (MQ=255)

                                                          ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAAc                                <  1:2531882/93‑1 (MQ=255)

                                                          ttACGACGAGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAAc                                <  1:6094602/93‑1 (MQ=255)

                                                           tACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAAcc                               <  2:1072106/93‑1 (MQ=255)

                                                            acgacgAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCg                              >  2:3247386/1‑93 (MQ=255)

                                                              gacgaGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGaa                            <  1:1555101/93‑1 (MQ=255)

                                                               acgaGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAg                           <  1:1625754/93‑1 (MQ=255)

                                                               acgaGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAg                           >  1:73050/1‑93 (MQ=255)

                                                                cgaGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGc                          <  2:4656943/93‑1 (MQ=255)

                                                                cgaGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGc                          <  1:6522135/93‑1 (MQ=255)

                                                                 gaGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGcc                         >  2:1564944/1‑93 (MQ=255)

                                                                 gaGGCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGcc                         >  2:3316768/1‑93 (MQ=255)

                                                                  aGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCg                        >  1:3254318/1‑93 (MQ=255)

                                                                   ggCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGa                       >  1:5687637/1‑93 (MQ=255)

                                                                    gCGGTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGat                      <  2:3521920/93‑1 (MQ=255)

                                                                    gCGGTGGCGAAAGTGTTTGAACCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGat                      >  2:656300/1‑93 (MQ=255)

                                                                     cGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGata                     >  2:2321370/1‑93 (MQ=255)

                                                                      ggTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGatat                    >  2:4807575/1‑93 (MQ=255)

                                                                      ggTGGCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGatat                    >  1:5028901/1‑93 (MQ=255)

                                                                        tGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTc                  >  2:5898599/1‑93 (MQ=255)

                                                                         ggCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCg                 <  1:108426/93‑1 (MQ=255)

                                                                         ggCGAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCg                 <  2:4586268/93‑1 (MQ=255)

                                                                           cGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGcc               <  2:4163969/93‑1 (MQ=255)

                                                                            gAAAGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCt              >  2:1011413/1‑93 (MQ=255)

                                                                              aaGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCtgc            >  1:2749467/1‑92 (MQ=255)

                                                                               aGTGTTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCTtgtg           >  2:615312/1‑93 (MQ=255)

                                                                                gtgtTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTgg          >  1:1654820/1‑93 (MQ=255)

                                                                                gtgtTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTgg          >  1:6407545/1‑93 (MQ=255)

                                                                                 tgtTTGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGa         <  2:4085858/93‑1 (MQ=255)

                                                                                   tttGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGACt       >  2:1448383/1‑93 (MQ=255)

                                                                                    ttGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGACTa      >  1:3150336/1‑93 (MQ=255)

                                                                                     tGAATCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGCCTAc     >  1:4992509/1‑93 (MQ=255)

                                                                                       aaTCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGACTACGc   <  2:257231/93‑1 (MQ=255)

                                                                                        aTCACTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGACTACGCt  >  1:5445542/1‑93 (MQ=255)

TTTATGACACCGTTAACAACGGCGTGTATAAAGCTGGTTTTGCCACCAGCCAGGAAGCTTACGACGAGGCGGTGGCGAAAGTGTTTGAATCGCTGGCACGACTGGAGCAGATTTTAGGCCAGCATCGTTACCTGACCGGCAACCAGCTAACCGAAGCCGATATTCGCCTGTGGACTACGCT  >  NC_000913/3249220‑3249400

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