BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	127,232	A→G	D513G (GAC→GGC)	aceF →	pyruvate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	127,232	0	A	G	99.6%	884.0 / ‑5.9	282	D513G (GAC→GGC)	aceF	pyruvate dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (115/165); total (115/167)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 292 total aligned reads displayed.

AAGACGGTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAA  >  NC_000913/127142‑127321

 aGACGGTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGt                                                                                        >  2:5612253/1‑93 (MQ=255)

   aCGGTCAGCGTCTGACCCTGCAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTTCGGTATTCAAAGGCGTca                                                                                      >  2:2169127/1‑93 (MQ=255)

     ggTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaac                                                                                    <  2:1959939/93‑1 (MQ=255)

      gTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaaca                                                                                   >  2:7140974/1‑93 (MQ=255)

       tCAGCGTCTGACCCTGAGGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaacaa                                                                                  <  1:7031498/93‑1 (MQ=255)

       tCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaacaa                                                                                  >  2:6767560/1‑93 (MQ=255)

       tCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaacaa                                                                                  <  1:2445212/93‑1 (MQ=255)

         aGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGa                                                                                <  1:4730333/93‑1 (MQ=255)

         aGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGa                                                                                >  2:3759398/1‑93 (MQ=255)

         aGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGa                                                                                <  1:750361/93‑1 (MQ=255)

         aGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGa                                                                                <  1:7763865/93‑1 (MQ=255)

          gCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGACCGGTCTGGTTGTTCCGGTATTCCAAGGCGTCAACAAGaa                                                                               >  1:266401/1‑93 (MQ=255)

          gCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAAAAAGaa                                                                               >  2:4591185/1‑93 (MQ=255)

           cGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGaaa                                                                              <  1:7178548/93‑1 (MQ=255)

           cGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGaaa                                                                              <  2:3355353/93‑1 (MQ=255)

            gTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAg                                                                             <  1:216044/93‑1 (MQ=255)

            gTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAg                                                                             <  2:6489764/93‑1 (MQ=255)

             tCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTTAAAGGCGTAAACAAGAAAgg                                                                            >  2:2379406/1‑93 (MQ=255)

             tCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAgg                                                                            >  2:6231095/1‑93 (MQ=255)

             tCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAgg                                                                            >  1:6108712/1‑93 (MQ=255)

              cTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGc                                                                           <  1:2326729/93‑1 (MQ=255)

               tGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGca                                                                          <  2:515853/93‑1 (MQ=255)

               tGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGca                                                                          >  2:4459165/1‑93 (MQ=255)

               tGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGca                                                                          >  1:7363693/1‑93 (MQ=255)

               tGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGca                                                                          <  2:277529/93‑1 (MQ=255)

               tGACCCTGAACAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGca                                                                          >  2:1522263/1‑93 (MQ=255)

                gACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcat                                                                         >  2:7250292/1‑93 (MQ=255)

                gACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcat                                                                         <  1:6989179/93‑1 (MQ=255)

                gACCCTGAAGAAATACATCAACATCGGTGGGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcac                                                                         >  1:2751630/1‑92 (MQ=255)

                 aCCCTGAAGAAATACATCAACATCGGTGTGGCGGTGTATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcatc                                                                        <  1:375287/93‑1 (MQ=255)

                 aCCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcatc                                                                        >  2:2640541/1‑93 (MQ=255)

                  cccTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcatca                                                                       <  2:154345/93‑1 (MQ=255)

                   ccTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcatcat                                                                      >  1:963749/1‑93 (MQ=255)

                     tGAAGAAATACATCAACATCGGTGTGGCGGTGGATACACCGAACGGGCTGGATGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCg                                                                    >  1:1036754/1‑93 (MQ=255)

                       aagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGACCGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAg                                                                  <  2:3519646/93‑1 (MQ=255)

                       aagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAg                                                                  <  2:5715282/93‑1 (MQ=255)

                        agaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGc                                                                 >  2:5969197/1‑93 (MQ=255)

                        agaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGc                                                                 <  2:3919473/93‑1 (MQ=255)

                         gaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCt                                                                >  2:5501168/1‑93 (MQ=255)

                          aaaTACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTg                                                               >  2:3092690/1‑93 (MQ=255)

                            aTACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGtc                                                             >  2:5010552/1‑93 (MQ=255)

                             tACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGtct                                                            <  2:3960108/93‑1 (MQ=255)

                             tACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGtct                                                            <  2:4109725/93‑1 (MQ=255)

                                 tCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGGATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTcgcg                                                        >  1:1386167/1‑93 (MQ=255)

                                  cAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGa                                                       <  1:5649838/93‑1 (MQ=255)

                                  cAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGa                                                       >  2:1445111/1‑93 (MQ=255)

                                    aCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGc                                                     <  1:1838910/93‑1 (MQ=255)

                                    aCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGc                                                     <  1:6339768/93‑1 (MQ=255)

                                    aCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGc                                                     <  1:3021355/93‑1 (MQ=255)

                                     cATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCt                                                    <  2:583276/93‑1 (MQ=255)

                                     cATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCt                                                    <  2:704927/93‑1 (MQ=255)

                                     cATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCt                                                    <  1:4311343/93‑1 (MQ=255)

                                       tCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtga                                                  <  1:6461732/93‑1 (MQ=255)

                                       tCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtga                                                  >  2:4846594/1‑93 (MQ=255)

                                        cGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtgat                                                 <  2:900151/93‑1 (MQ=255)

                                        cGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtgat                                                 <  1:1957722/93‑1 (MQ=255)

                                          gtgtGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtgatga                                               >  1:4145157/1‑93 (MQ=255)

                                          gtgtGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtgatga                                               >  2:2477503/1‑93 (MQ=255)

                                          gtgtGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtgatga                                               >  1:286527/1‑93 (MQ=255)

                                          gtgtGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtgatga                                               <  2:6054443/93‑1 (MQ=255)

                                           tgtgGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGAc                                              <  1:6631327/93‑1 (MQ=255)

                                            gtgGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACt                                             <  2:3786276/93‑1 (MQ=255)

                                            gtgGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACt                                             <  1:5313618/93‑1 (MQ=255)

                                             tgGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTa                                            <  2:1115122/93‑1 (MQ=255)

                                               gCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTAtt                                          <  2:1840135/93‑1 (MQ=255)

                                                    ggATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTaa                                     >  2:6357999/1‑93 (MQ=255)

                                                     gATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAg                                    >  1:6909751/1‑93 (MQ=255)

                                                        aCCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGaaa                                 <  2:123807/93‑1 (MQ=255)

                                                        aCCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGaaa                                 <  2:783726/93‑1 (MQ=255)

                                                          cccGACTGGTCTGTTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAgc                               <  2:2811493/93‑1 (MQ=255)

                                                          cccGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAgc                               <  1:4682155/93‑1 (MQ=255)

                                                           ccGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAgcg                              >  1:6224598/1‑93 (MQ=255)

                                                             gAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGcgcg                            <  1:3381560/93‑1 (MQ=255)

                                                             gAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGcgcg                            >  1:508296/1‑93 (MQ=255)

                                                               aCGGTCTGGTTGTTCCGGGATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTg                          >  1:6662866/1‑93 (MQ=255)

                                                                cGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGa                         <  1:5845923/93‑1 (MQ=255)

                                                                 ggTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGAc                        <  1:2250907/93‑1 (MQ=255)

                                                                   tCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACgg                      >  2:2040801/1‑93 (MQ=255)

                                                                     tGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTa                    <  1:1485405/93‑1 (MQ=255)

                                                                      ggTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTaa                   <  1:3584233/93‑1 (MQ=255)

                                                                       gttgttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAg                  >  1:2967136/1‑93 (MQ=255)

                                                                        ttgttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGc                 <  1:1103101/93‑1 (MQ=255)

                                                                        ttgttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGc                 <  2:7531325/93‑1 (MQ=255)

                                                                         tgttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCt                <  1:5923758/93‑1 (MQ=255)

                                                                         tgttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCt                <  1:6444608/93‑1 (MQ=255)

                                                                           ttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGa              <  1:1790909/93‑1 (MQ=255)

                                                                           ttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGa              >  2:1617234/1‑93 (MQ=255)

                                                                            tCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGAc             <  2:7910464/93‑1 (MQ=255)

                                                                             ccGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACt            >  2:157643/1‑93 (MQ=255)

                                                                             ccGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACt            <  1:3232396/93‑1 (MQ=255)

                                                                             ccGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACt            <  2:169803/93‑1 (MQ=255)

                                                                              cGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTg           >  2:655404/1‑93 (MQ=255)

                                                                                 tATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCgg        <  1:4860518/93‑1 (MQ=255)

                                                                                  aTTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCggg       >  2:6320594/1‑93 (MQ=255)

                                                                                  aTTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCggg       <  1:6564968/93‑1 (MQ=255)

                                                                                  aTTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCggg       <  2:7358416/93‑1 (MQ=255)

                                                                                   ttCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGc      <  1:824101/93‑1 (MQ=255)

                                                                                   ttCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGc      >  2:2343416/1‑93 (MQ=255)

                                                                                      aaaGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAGGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGaa   >  1:1174194/1‑93 (MQ=255)

                                                                                       aaGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGaaa  <  1:7853236/93‑1 (MQ=255)

AAGACGGTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAA  >  NC_000913/127142‑127321

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