BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,788,398	A→G	V192A (GTG→GCG)	tdh ←	threonine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,788,398	0	A	G	83.2%	642.0 / 114.5	298	V192A (GTG→GCG)	tdh	threonine dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (24/26); new base G (128/119); total (152/146)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

Only 100 of 315 total aligned reads displayed.

CGTCATTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGA  >  NC_000913/3788306‑3788490

  tCATTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATcacc                                                                                            >  2:608722/1‑93 (MQ=255)

   cATTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATcgcca                                                                                           <  2:4202600/93‑1 (MQ=255)

   cATTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATcacca                                                                                           <  1:5716030/93‑1 (MQ=255)

    aTTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATcgccac                                                                                          >  2:1513301/1‑93 (MQ=255)

    aTTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATcgccac                                                                                          >  1:3783982/1‑93 (MQ=255)

    aTTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATcaccac                                                                                          >  2:1617882/1‑93 (MQ=255)

     ttGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCcaca                                                                                         >  1:7034845/1‑93 (MQ=255)

      tGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAt                                                                                        <  2:5522435/93‑1 (MQ=255)

      tGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAt                                                                                        <  2:6480037/93‑1 (MQ=255)

      tGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAt                                                                                        <  1:1615340/93‑1 (MQ=255)

      tGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAt                                                                                        >  1:1815015/1‑93 (MQ=255)

      tGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAt                                                                                        <  1:5113136/93‑1 (MQ=255)

      tGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAt                                                                                        >  2:6634015/1‑93 (MQ=255)

       nngaTTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAtt                                                                                       <  2:1384876/91‑1 (MQ=255)

       gagaTTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAtt                                                                                       >  1:3546036/1‑93 (MQ=255)

       gagaTTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAtt                                                                                       <  1:847059/93‑1 (MQ=255)

       gagaTTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACAtt                                                                                       <  1:4091358/93‑1 (MQ=255)

         gaTTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGc                                                                                     <  1:6399371/93‑1 (MQ=255)

           ttttCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGTTCGCCACATTGCGt                                                                                   <  2:3182121/93‑1 (MQ=255)

            tttCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTg                                                                                  >  2:1105827/1‑93 (MQ=255)

             ttCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGc                                                                                 >  1:5523595/1‑93 (MQ=255)

               cTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCAc                                                                               >  2:2352453/1‑93 (MQ=255)

               cTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCAc                                                                               <  1:2241878/93‑1 (MQ=255)

                  tGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCaa                                                                            <  2:5300097/93‑1 (MQ=255)

                    gCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACg                                                                          >  1:7482341/1‑93 (MQ=255)

                      gACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACgtg                                                                        <  1:3658639/93‑1 (MQ=255)

                       aCGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACgtgt                                                                       <  1:4212097/93‑1 (MQ=255)

                       aCGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACgtgt                                                                       >  1:5460416/1‑93 (MQ=255)

                        cGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGtt                                                                      >  1:6388971/1‑93 (MQ=255)

                        cGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGtt                                                                      >  1:1958982/1‑93 (MQ=255)

                         gTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGttt                                                                     >  1:6527640/1‑93 (MQ=255)

                          ttAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTc                                                                    >  1:1286931/1‑93 (MQ=255)

                               cGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCAc                                                               <  2:2478064/93‑1 (MQ=255)

                                gCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCAcc                                                              >  2:4440444/1‑93 (MQ=255)

                                 cACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCAccg                                                             >  2:3504066/1‑93 (MQ=255)

                                  aCGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCAccgc                                                            >  1:4797288/1‑93 (MQ=255)

                                    gggTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCAccgccg                                                          >  1:6280017/1‑93 (MQ=255)

                                    gggTGATACCCATTCTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCAccgccg                                                          >  2:6443256/1‑93 (MQ=255)

                                          tACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCa                                                    >  2:5151144/1‑93 (MQ=255)

                                          tACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCa                                                    >  1:4240870/1‑93 (MQ=255)

                                          tACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACTTGTTTCGCCCCCGCCGCTGCCa                                                    <  1:461981/93‑1 (MQ=255)

                                           aCCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCat                                                   <  1:4985528/93‑1 (MQ=255)

                                           aCCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCat                                                   >  1:445837/1‑93 (MQ=255)

                                           aCCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCat                                                   >  2:3972693/1‑93 (MQ=255)

                                            cccATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCata                                                  >  1:6581594/1‑93 (MQ=255)

                                             ccATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCataa                                                 >  1:2116249/1‑93 (MQ=255)

                                             ccATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCataa                                                 <  1:3848989/93‑1 (MQ=255)

                                               aTTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGTACCAACGTGTTTCGCCACCGCCGCTGCCataata                                               <  2:5107036/93‑1 (MQ=255)

                                               aTTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCataata                                               <  2:365659/93‑1 (MQ=255)

                                                ttttACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATAc                                              >  2:740430/1‑93 (MQ=255)

                                                   tACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCaa                                           <  1:4693020/93‑1 (MQ=255)

                                                     cgcgCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATc                                         <  1:3177918/93‑1 (MQ=255)

                                                     cgcgCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATc                                         <  2:6786483/93‑1 (MQ=255)

                                                        gcCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCggg                                      >  2:2213376/1‑93 (MQ=255)

                                                        gcCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCggg                                      >  1:966183/1‑93 (MQ=255)

                                                        gcCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCggg                                      >  1:6311741/1‑93 (MQ=255)

                                                        gcCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCggg                                      >  1:4547449/1‑93 (MQ=255)

                                                         ccAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGc                                     <  2:52250/93‑1 (MQ=255)

                                                         ccAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGc                                     >  1:3278075/1‑93 (MQ=255)

                                                          cAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGcc                                    <  2:7468242/93‑1 (MQ=255)

                                                          cAGCTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGcc                                    >  1:3906173/1‑93 (MQ=255)

                                                             cTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGc                                 <  1:6458027/93‑1 (MQ=255)

                                                             cTCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGc                                 >  1:1132821/1‑93 (MQ=255)

                                                              tCAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCa                                >  1:6909484/1‑93 (MQ=255)

                                                               cAAGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCAc                               <  1:1046956/93‑1 (MQ=255)

                                                                aaGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCAcc                              >  2:5964030/1‑93 (MQ=255)

                                                                aaGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCAcc                              >  1:108150/1‑93 (MQ=255)

                                                                aaGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCAcc                              >  1:6153852/1‑93 (MQ=255)

                                                                 aGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCa                             <  2:65540/93‑1 (MQ=255)

                                                                 aGGCGGTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCa                             >  2:1058307/1‑93 (MQ=255)

                                                                     ggTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGaaa                         >  2:1856547/1‑93 (MQ=255)

                                                                      gTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAAc                        <  2:6155760/93‑1 (MQ=255)

                                                                      gTATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAAc                        >  1:682695/1‑93 (MQ=255)

                                                                       tATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAAcc                       >  2:3750025/1‑93 (MQ=255)

                                                                       tATTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAAcc                       >  1:2779357/1‑93 (MQ=255)

                                                                       tATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAAcc                       >  2:2095156/1‑93 (MQ=255)

                                                                        aTTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCa                      <  1:4839450/93‑1 (MQ=255)

                                                                        aTTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCa                      <  1:624780/93‑1 (MQ=255)

                                                                        aTTCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCa                      >  2:4687703/1‑93 (MQ=255)

                                                                         ttCGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAg                     <  2:7287621/93‑1 (MQ=255)

                                                                         ttCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAg                     <  1:6054510/93‑1 (MQ=255)

                                                                          tCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGc                    <  2:6194881/93‑1 (MQ=255)

                                                                           cGTTAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGca                   <  2:7043480/93‑1 (MQ=255)

                                                                             ttAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGcaca                 >  2:845586/1‑93 (MQ=255)

                                                                             ttAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGcaca                 <  2:711967/93‑1 (MQ=255)

                                                                              tAACATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACAt                >  1:2522814/1‑93 (MQ=255)

                                                                               aaCATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATc               <  1:5954114/93‑1 (MQ=255)

                                                                               aaCATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATc               <  1:4914181/93‑1 (MQ=255)

                                                                                aCATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCt              >  2:4966720/1‑93 (MQ=255)

                                                                                aCATCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCt              <  1:721629/93‑1 (MQ=255)

                                                                                  aTCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTc            <  2:2942017/93‑1 (MQ=255)

                                                                                  aTCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTc            >  2:2015650/1‑93 (MQ=255)

                                                                                  aTCAGTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTc            <  1:6655193/93‑1 (MQ=255)

                                                                                  aTCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTc            <  1:7308905/93‑1 (MQ=255)

                                                                                      gTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGccc        >  2:4305888/1‑93 (MQ=255)

                                                                                      gTGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGccc        >  1:5346866/1‑93 (MQ=255)

                                                                                       tGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCcca       >  2:4167892/1‑93 (MQ=255)

                                                                                       tGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCcca       >  2:329586/1‑93 (MQ=255)

                                                                                       tGATCGCCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCcca       >  2:261509/1‑93 (MQ=255)

                                                                                            accacATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGa  >  1:1444167/1‑93 (MQ=255)

CGTCATTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGA  >  NC_000913/3788306‑3788490

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