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.4%	548.7 / NA	176	D513G (GAC→GGC)	aceF	pyruvate dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (70/104); total (71/105)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

AGACGGTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAA  >  NC_000913/127143‑127321

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

 gACGGTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTc                                                                                       >  2:3996605/1‑93 (MQ=255)

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

    ggTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaac                                                                                    >  2:700691/1‑93 (MQ=255)

    ggTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACACCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaac                                                                                    >  2:284671/1‑93 (MQ=255)

     gccAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaaca                                                                                   <  2:3177182/91‑1 (MQ=255)

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

     gTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaaca                                                                                   <  1:264481/93‑1 (MQ=255)

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

     gTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTcaaca                                                                                   <  2:4258845/93‑1 (MQ=255)

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

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

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

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

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

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

              tGACCCTGAAGAAATACATTAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGca                                                                          >  2:4562905/1‑93 (MQ=255)

                aCCCTTAAGAAATACATCACCATCGGTGTGGCTGTGGATTCCCCGAACGGTCTGGTTGTACCGGTATTCAAAGGCGTCAACAAGAAAGGcatc                                                                        <  2:2207622/93‑1 (MQ=255)

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

                aCCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGcatc                                                                        <  2:1829563/93‑1 (MQ=255)

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

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

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

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

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

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

                   cTGAAGAAATACATCATCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCatcatc                                                                     <  1:4315418/93‑1 (MQ=255)

                    tGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCg                                                                    <  2:1904438/93‑1 (MQ=255)

                    tGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCg                                                                    >  1:1236723/1‑93 (MQ=255)

                     gaagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGa                                                                   <  2:2017790/93‑1 (MQ=255)

                     gaagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGa                                                                   >  1:2425595/1‑93 (MQ=255)

                     gaagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGa                                                                   <  2:3539396/93‑1 (MQ=255)

                     gaagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGCCATCATCGa                                                                   >  2:2334911/1‑93 (MQ=255)

                      aagaaATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAg                                                                  >  1:524699/1‑93 (MQ=255)

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

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

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

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

                           aTACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGtc                                                             <  1:92140/93‑1 (MQ=255)

                             aCATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGtctc                                                           <  2:3344953/93‑1 (MQ=255)

                             aCATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGtctc                                                           <  2:1552630/93‑1 (MQ=255)

                              cATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGCCGTCAACAAGAAAGGCATCATCGAGCTGTCTcg                                                          <  1:4568797/93‑1 (MQ=255)

                               aTCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTcgc                                                         >  2:2400964/1‑93 (MQ=255)

                                tCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTcgcg                                                        <  2:1420156/93‑1 (MQ=255)

                                tCAACATCGGTGCGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTcgcg                                                        >  2:593128/1‑93 (MQ=255)

                                  aaCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAg                                                      >  1:4490786/1‑93 (MQ=255)

                                   aCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGc                                                     >  1:1372566/1‑93 (MQ=255)

                                   aCATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGc                                                     >  1:3123916/1‑93 (MQ=255)

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

                                    cATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCt                                                    >  1:3240154/1‑93 (MQ=255)

                                     aTCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtg                                                   >  2:1473608/1‑93 (MQ=255)

                                     aTCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCtg                                                   >  2:2683455/1‑93 (MQ=255)

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

                                          tgtgGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGAc                                              >  1:389910/1‑93 (MQ=255)

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

                                           gtgGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACt                                             >  1:2168041/1‑93 (MQ=255)

                                             ggCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTAt                                           <  1:1304164/93‑1 (MQ=255)

                                             ggCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTAt                                           <  2:2962410/93‑1 (MQ=255)

                                             ggCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTAt                                           >  1:3041773/1‑93 (MQ=255)

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

                                               cGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTAttt                                         <  2:60929/93‑1 (MQ=255)

                                                ggTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTc                                        <  2:3635128/93‑1 (MQ=255)

                                                ggTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTc                                        <  1:3449903/93‑1 (MQ=255)

                                                ggTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTc                                        <  1:203243/93‑1 (MQ=255)

                                                 gTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCt                                       <  2:2510022/93‑1 (MQ=255)

                                                  tGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGAATCATCGAGCTGTCTCGCGAGCTGATGACTCTTTCTa                                      >  1:1077723/1‑93 (MQ=255)

                                                   ggATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTaa                                     <  1:4343912/93‑1 (MQ=255)

                                                     aTACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGa                                   <  1:956258/93‑1 (MQ=255)

                                                     aTACCACGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGa                                   <  2:1550155/93‑1 (MQ=255)

                                                         cccGAGCGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAgc                               <  2:3978405/93‑1 (MQ=255)

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

                                                           cGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAgcgc                             >  1:3148811/1‑93 (MQ=255)

                                                           cGAACGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAgcgc                             >  1:1852674/1‑93 (MQ=255)

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

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

                                                             aaCGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGt                           <  1:3380094/93‑1 (MQ=255)

                                                             aaCGGTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGt                           <  1:175371/93‑1 (MQ=255)

                                                                ggTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGAc                        >  2:710289/1‑93 (MQ=255)

                                                                 gTCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACg                       <  1:69440/93‑1 (MQ=255)

                                                                  tCTGGTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACgg                      <  2:4669702/93‑1 (MQ=255)

                                                                     ggTTGTTCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTaa                   >  1:3005755/1‑93 (MQ=255)

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

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

                                                                       ttgttCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATTGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGc                 <  1:4469965/93‑1 (MQ=255)

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

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

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

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

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

                                                                           tCCGGTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGAc             >  1:812180/1‑93 (MQ=255)

                                                                              ggTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGc          <  2:2441173/93‑1 (MQ=255)

                                                                              ggTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGc          >  2:1173419/1‑93 (MQ=255)

                                                                              ggTATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGc          <  1:2457209/93‑1 (MQ=255)

                                                                                tATTCAAAGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTAACTGCgg        >  1:1454823/1‑93 (MQ=255)

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

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

                                                                                     aaaGGCGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGaa   >  2:808056/1‑93 (MQ=255)

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

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

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

AGACGGTCAGCGTCTGACCCTGAAGAAATACATCAACATCGGTGTGGCGGTGGATACCCCGAACGGTCTGGTTGTTCCGGTATTCAAAGACGTCAACAAGAAAGGCATCATCGAGCTGTCTCGCGAGCTGATGACTATTTCTAAGAAAGCGCGTGACGGTAAGCTGACTGCGGGCGAAA  >  NC_000913/127143‑127321

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